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Zur nichtlinearen Torsion eines Kreiszylinders * 


Von E. Koppe 


1. Einleitung. Die nichtlineare Torsion eines Kreiszylinders ist bereits mehrfach behandelt 
worden, z. B. von R. S. Rivlin!, A. E. Green and R. T. Shield?. Die einfachste Form der Behand- 
lung scheint jedoch bisher nicht erreicht zu sein. Im folgenden soll auf theoretische und experi- 
mentelle Methoden aufmerksam gemacht werden, die dem Problem besonders gut angepaft sind. 
| Der als inkompressibel vorausgesetzte gerade Kreiszylinder soll einer ,,reinen‘* Torsion unter- 

worfen werden, d.h., die Querschnitte sollen eben bleiben. WVerhindert man eine Streckung des 
gesamten Zylinders, so bleiben infolge der Inkompressibilitat auch die radialen Abstande erhalten. 
Gefragt wird nach der Krafteverteilung, die man an den Endflachen des Zylinders anbringen muB, 
um diesen Verfolgungszustand zu erreichen. 


2. MaBtensoren und Christoffelsymbole bei der Torsion. In Zylinderkoordinaten (siehe Abb. 1) 
erhalt man fiir den unverformten Zustand als Quadrat des Linienelementes 


ds = g,, dq* dq? = dr? + r2 dq? + dz. (1) 


(Es wird von der Summationskonvention Gebrauch gemacht; iiber gleiche Indizes, die einmal 

oben (kontravariant) und einmal unten (kovariant) auftreten, wird summiert.) Das Zeichen ° 

bezeichnet GréBen, die auf den unverformten Kérper bezogen sind, z. B. 
-8xp = MaBtensor des unverformten Kérpers. 

Nach einer Drehung im positiven Sinn um den Winkel w z erhalt man J Le» 

fiir den verformten Zustand als Quadrat des Linienelementes 


ds? = g,, dq* dq? = dr? + 1? dy? + 21? w dp dz + (1 + 1? w?) dz*. (2) 


Im folgenden werden Glieder mit m und wm? mitgenommen, dagegen 
héhere Potenzen von ~ vernachlassigt. 


2 
Spater benétigt man bei der Differentiation nach den Koordinaten des * 
verformten Kérpers die Christoffelsymbole zweiter Art ['?,, die aus den 
Christoffelsymbolen erster Art 
_ 1 (ac ap | Agpo\ 
Dae, she ee dqo | Age (3) 
a 


durch Uberschieben mit dem reziproken MaBtensor g” hervorgehen: 


faa 1 rs GLY (4) 


Abb. 1. Die benutzten 


Von Null verschieden sind hier allein die GréBen Koordinaten, 
i= — FL=—ro, Th=—ro’, I? ee (5) 
op ee? 23 ? ae = oeeech 12 5) 13 = 


Sind g = |g,,| und g = |Z,g| die Determinanten der Maftensoren, so wird ganz allgemein 
der Quotient der Volumenelemente des verformten und unverformten Kérpers durch 


ave 2 ae 
ee aired (6) 


ave 
gegeben. Aus dem in der Literatur fast ausschlieBlich benutzten Greenschen Verzerrungstensor 
2 ap = 2.8 ~ Bap (7) 


* Zusammenfassung eines Vortrages, der im Mechanik-Koiloquium der Niedersichsischen Hochschulen, 
Gottingen am 25. Januar 1955 gehalten wurde. 

1R.S. Rivlin, J. Appl. Phys. 18 (1947) S. 444 u. 837. 

2 4, E. Green und R. T. Shield, Proc. Roy. Soc. London A 202 (1950) S. 407. 
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erhalt man durch Uberschieben mit dem reziproken MaBtensor des unverformten Kérpers g* 
2 je = 9 g,,— 08. (8) 


Die Volumendehnung wird hier durch den unhandlichen Ausdruck 


dy Rae 
itm /\65 +275 | 


= : (9) 
+ ey of Op Oy i) Oy Ox 4 oy 2a © 4 Ke og 3 yom 2 
= \1 | 2 Vy 2 Vy og Tw: 2 Vn Y» =F ey toy Ne Ye 123% 4 Vy Yo Vx a com V6 Vo ys 


beschrieben. Im Falle der reinen Torsion ist nun 


6 dV 
Ce Gly Iie, Pe iat (10) 


Es ware deshalb von Vorteil, wenn man die Torsion durch einen Verzerrungstensor beschreiben 
kénnte, bei dem ganz allgemein eine einfache Invariante die Volumendehnung kennzeichnet, wie 
es z. B. im linearen Fall sehr schwacher Formanderung beim Greenschen Verzerrungstensor 


geschieht : 

dV oe 
rag ahr (11) 
dv 


3. Der Henckysche (logarithmische) Verzerrungstensor. Definiert man den Logarithmus einer 


Matrix 9{ mit den Elementen as durch 


eine — OI ; 
und bezeichnet man die Komponenten von In Y mit L az , so gilt allgemein 
Lat =In |as| ry (12) 


d.h. die erste Invariante, die ,,Spur“, des Logarithmus einer Matrix ist gleich dem Logarithmus 


ihrer Determinante. : 
Sucht man nach einem Verzerrungstensor, beidem auchim nichtlinearen Fall die Volumen- 
dehnung allein durch seine erste Invariante beschrieben wird, so wird man durch die Beziehung (12) 


dazu veranlaft, einfach den Logarithmus von (6) zu bilden: 


dV 1 
i => > in B™ 8 up| 
und den Tensor 
2 hs = Lg aye (13) 
als neuen Verzerrungstensor einzufiihren. Fiir diesen Verzerrungstensor ist 
1 ° 1 ° dV 
h’ a oy Lg Sae 3 in greg = gs 


Mit dem Greenschen Verzerrungstensor haingt er durch die Beziehung 


he = L (62 + 2 93) (15) 


pe2 
zusammen. Seine erste Invariante kennzeichnet die Volumendehnung. Der Tensor h®, der als 


Henckyscher Verzerrungstensor! bezeichnet werden soll, 1a8t sich auch auf andere als die gerade 
gegebene Weise einfiihren. So wird man z. B. auf ihn gefiihrt, wenn man ein Superpositionsprinzip 
fiir koaxiale Streckungen fordert oder wenn man die Verformung multiplikativ aus Drehung und 
Streckung zusammensetzen will. Wird hy auf Hauptachsen gebracht, so stellen seine Komponenten 


gerade die Logarithmen der Quotienten von Endlange durch Anfangslange der Strecken in Haupt- 
achsenrichtung dar. Es sei im folgenden eine andere, mehr formale Einfithrung gegeben, die die 
Gleichberechtigung von Greenschem und Henckyschem Verzerrungstensor besonders gut hervor- 
treten laBt?. 

Bezeichnet man den auf den verformten Kérper bezogenen Spannungstensor bei Zerlegung der 
Spannungen nach den Koordinatenrichtungen mit 1% und die Dichte des verformten Kérpers mit 9, 


Vel. C. Truesdell, J. Rat. Mech. and Analysis 1 (1952) S. 125; 2 (1953) S. 593. 
* Man vergleiche auch E. Koppe, Z. angew. Math. Mech. 36 (1956) S. 455. 


3 


Koppe: Zur nichtlinearen Torsion eines Kreiszylinders 


XXYV. Band 1957 
so ist die Anderung der spezifischen Formanderungsenergie w bei Benutzung des Greenschen Ver- 


zerrungstensors gegeben durch 
e dw = tv dys . 


Geht man zum gemischtvarianten Spannungstensor iiber, d. h. zerlegt man die Spannungen nach 
(16) 


Richtungen senkrecht zu den Koordinatenflachen, so erhalt man 
e dw = 7? g dyog . 


Da man Differentialbildung und Uberschiebung mit dem Ma8tensor des verformten Korpers nicht 
vertauschen darf, kann man hier nicht in einfacher Weise einen aus dem Greenschen Verzerrungs- 


tensor abgeleiteten gemischtvarianten Tensor einfiihren. 


Mit (7) la&t sich jedoch fiir (16) schréiben 
e dw = 78 ge is 


und aus’ 


folgt rein formal 
Man erhalt so 

o dw =t! dhj. (17) 
Der Henckysche Verzerrungstensor erscheint hier in natiirlicher Weise dem gemischtvarianten 
Spannungstensor zugeordnet, wie der Greensche Verzerrungstensor dem rein kontravarianten Span- 


nungstensor zugeordnet ist. 
4, Das Elastizitatsgesetz bei inkompressiblem Material. Sind die dh; unabhangig vom Verfor- 


Ow 


mungsweg, so ist dw ein vollstandiges Differential der Verzerrungen, d. h. 
E 
o dw =o ape oP 
B 


Ein Vergleich mit (17) liefert 
9 =e ahs =r dhs. (18) 


oh’, 
B 
Ist der Kérper inkompressibel, so sind die dh* nicht mehr vom Verformungswege unabhangig, da 


stets die Volumendehnung ungeandert bleiben muB 
dh; = 56 dhg = 0. 


Man erhalt so anstelle von (18) mit einem beliebigen Lagrangeschen Multiplikator p 


ve dh; = (¢ - 
oder mit 9 = 0 
pie a OF (19) 
% Oh‘, 


Die Formanderungsenergie mége allein von den Invarianten der Verzerrungen abhangen. Die ein- 
. & é x €7a 7,0 . 
fachsten Invarianten, aus denen sich alle anderen zusammensetzen lassen, sind h, , hy he hte hes 


. . & 
da Inkompressibilitat vorausgesetzt ist, verschwindet h, 


Mit 
é(0 w) VD a(0 w) tape 
SUING See et te 
ane he 8 arene he 8 
Bi pdh 4 Ohe 4 Pre he. (2 


erhalt man so aus (19) 
Te 
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5. Verzerrungstensor und Spannungstensor bei der Torsion. Die Gleichung (13) gibt zusammen 
mit den Gleichungen (1) und (2) fiir die MaBtensoren im Fall der reinen Torsion den Henckyschen 


Verzerrungstensor 
i @ @\ Oh OQ) 
Dee elie OFS Saris (21) 
ORM 0 ra r2m2 


Bevor durch Reihenentwicklung bis zu Gliedern mit w* der Tensor (h®) berechnet wird, sollen die 
Eigenwerte dieses Tensors ermittelt werden. Die Eigenwerte A von (g g,-) folgen aus der Gleichung 


(1 —A) ?—A24+7?o") + 1] =0. 


(he) = eS In 


Ks ist 
yirees dl. poy lk 


Die Eigenwerte 7 von (he) sind die Logarithmen der Eigenwerte von (g g,,); d. h. es ist 


my =9, No + 43 = 9. 


In der Gleichung 


—7 + J, 7? — Jon + Jz = 9 
fiir die Kigenwerte 7 ist 
é€ 1 Ea é€ Xx 1 Ej,u Oo E€E 7K 7O E 7a 1o 
Jy =he, Jp =x (tte My — hry he) » Jy = G (he he hg — 3 he Hg hy + 2 h,, he h.) - 


Man erkennt, daB im Falle der reinen Torsion, wo h, = 0 ist, auch h, hj h? = 0 wird; d. h. allein 


die Invariante h;, h; ist von Null verschieden. Diese iibersichtlichen Verhaltnisse treten bei Ver- 
wendung des Henckyschen Verzerrungstensors stets ein, wenn inkompressibles Material so verformt 
wird, daB die Dehnung in einer Hauptdehnungsrichtung — hier in r-Richtung — verschwindet. 


In (20) sind also ® und Y Funktionen der einzigen nichtverschwindenden Invarianten hehe 
Entwickelt man in (21) bis zu Gliedern mit jw? so ergibt sich 


0 0 0 \ 
1 1 
0 ——roa — w 
(he) = : Z | (22) 
2 2 2 
0 zo ZrO ) 
. . . . . E X 1 
Die einzige Invariante ist h, hg => rw? . 


Mit Hilfe von (22) und (20) kann man jetzt die Spannungen bis zu Gliedern mit w* anschreiben. 
Die Funktionen ® und Y’ reduzieren sich in dieser Naherung auf Konstanten. 


il © © 0 0 0 0 0 0 | 
1 1 1 2 
(cf) =p|° ! he Oe Oa | rege (23) 
1 iL 
0 Roe 2 2S aD ao 
OM 0 YO @ 0 09 [re 
6. Die Gleichgewichts- und Oberflachenbedingungen. Die Gleichgewichtsbedingungen lauten | 
ax? B 
€ 1B v Bey pe ¢ 
aie + iS oe sae se (), (24) 


Setat Senin ae und (5) ein, so verschwinden fiir ¢ = gy und ¢ = 2 alle Glieder, wihrend man fiir 
€ =r erhalt 


ope. t iba 


1 1 t 0 m 
pile Pe Ip 4 OO eo —“.>Po@=0, 5 = O+¥). (25) 


Der Zylindermantel muB spannungsfrei sein. Diese Forderung vereinfacht sich wegen le =t,=0 | 


zu tT, =p = 0 fir r=a und es folgt 


9 
wm 


pe sgl a) Omar: (26) 
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Man kann nun die an den Endflachen angreifenden Krafte bestimmen. Dazu geht man von den 
bisher benutzten reziproken Basisvektoren wieder zu den Basisvektoren in Richtung der Koordi- 
naten iiber. (Von der Summationskonvention soll im folgenden kein Gebrauch mehr gemacht wer- 
den; es werden Summationszeichen gesetzt.) Die in azimutaler Richtung an den Endflachen an- 


greifende Kraft ist gegeben durch 
9° ~ Van) os & 
Pe => Ti 8’? n, + héhere Glieder = 1”? n, , 


die Kraft senkrecht zur Endflaiche durch 


- ee 1 nS . om 
pS vt; g”*n, + héhere Glieder = t* n, . 


Mit den aus (23) und (26) resultierenden Spaunungen 


Fe ol hanes 
eras a (27) 
und 
9 2 
On Te. 


Bog Moe Oe a Ong) (28) 


erhalt man folgendes Ergebnis: Um bei einem geraden Kreiszylinder aus inkompressiblem Material 
reine Torsion, d. h. Ebenbleiben der Querschnitte zu erzwingen, mu® an den Endflachen in azimu- 
taler Richtung die Kraft 


pe = - w+ @ + Glieder der Ordnung w?, (29) 
in normaler Richtung die Kraft! 

Pe = © (Fa) (3 ¥— 9) — Tota (@— P) (30) 
angebracht werden. Dabei ist zu beachten, daB P® auf ey mit | ey | = r bezogen ist. 


= @ ist gleich G, dem Schubmodul der linearen Theorie. Die Konstante Y tritt in der linearen 


2 


Theorie tiberhaupt nicht auf und kann nur durch ein Experiment bestimmt werden. Aus der oben 
- : : es 1 
angegebenen Krafteverteilung folgt, da ® = 2G positiv ist: wenn auch VY positiv und ¥Y > zy @ 


ist, so miissen an der Endflache von innen nach auBen abnehmende Druckkrafte angebracht wer- 
den, sofern man reine Torsion er- 
halten will. Solange Y < @ ist, 
miissen auf der gesamten End- 
flache einschlieBlich des Randes 
Druckkrafte angebracht werden. 
Man ist also bereits mit wenigen Abstanus- 
qualitativen Aussagen in der 
Lage, eine Abschatzung fiir die 
neu auftretende Elastizitatskon- 
stante Y zu geben. 


aN 


A) Oftiaung fir Halteschraube | 
4 Torsionshebe/ 


7. Beschreibung einer Ver- 
suchseinrichtung zur Torsion 
eines geraden Kreiszylinders. Das 
Versuchsstiick ist ein Gummi- 
zylinder yon 110mm _ Durch- 
messer und 25mm Héhe, an —~—_ ee oe pe 
dessen Endflachen Messingplat- Abb. 2. Der Torsionsapparat (schematisch). 
ten aufvulkanisiert sind. In der ae 
oberen Messingplatte befinden sich drei MeBlécher von 12 mm Durchmesser, durch die die obere 
Endflache des Gummizylinders zuginglich ist. 

Dieses Versuchsstiick wird in einem Torsionsapparat, von dem Abb. 2 eine schematische Dar- 
stellung zeigt, auf einem drehbar gelagerten Sockel befestigt. Die obere Messingplatte wird bei der 
Torsion von einem Deckflansch festgehalten. Die Abstandshiilsen, die den Deckflansch tragen, 
sind gerade so lang, daf im unverformten Zustand kein Druck auf das Versuchsstiick ausgetibt 


Gs Truesdell schlug die Bezeichnung ,, Poyntingeffekt* fir diese Erscheinung vor. 
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wird. Dies wurde durch Parallelversuche mit stark und schwach angezogenen Deckflanschschrau- 
ben kontroilliert. Das Verdrillen geschieht mit Hilfe eines Spannschlosses; der Verdrillungswinkel 
kann an einer Skala abgelesen werden. 

Die Verhinderung der Verschiebung der Endflachen in achsialer Richtung bewirkt, das an den 
Endflichen bei der Torsion nicht nur Drill- und Einspannmoment, sondern auch Krafte in senk- 
rechter Richtung auf den sich auszudehnen trachtenden Gummizylinder ausgetibt werden. An den 
Stellen, an denen sich in der oberen Messingscheibe die MeBlécher befinden, fehlen diese Krafte und 
man erwartet, daB sich der Gummi bei Verdrillung in diese MeBlécher hinein auswolbt. Diese Aus- 
wilbung, die als MaB fiir die dort angreifenden Krafte angesehen wird, soll ermittelt werden. 


aoe 8. Die Mefeinrichtung. Bei den 
von R. S. Rivlin’ beschriebenen Ver- 
suchen — wie auch bei Versuchen 
On anderer Autoren — wurden zur Mes- 


Lutt vom Kompressor fesie Diise Dp Mehdise 


sung solcher Auswélbungen Meb- 

/Tiohe, deren uhren benutzt. Mit MeBuhren lassen 

ee es sich relative Anderungen von 10mm 

werden soll noch gut feststellen, die Kraft zur 

Betatigung der Uhr mu jedoch von 

Abb. 3. Prinzipskizze des MeByerfahrens. der sich verwélbenden Gummiflache 


aufgebracht werden, so dafi die Mes- 
sung die zu messende GréBe beeinfluBt. Ein weiterer Nachteil der Verwendung von MeBuhren 
besteht in der unvermeidlichen Reibung der Teile der Uhr gegeneinander. Bei der Messung elasti- 
scher Verschiebungen entsteht dadurch ein Spiel zwischen elastischer Spannung und Haft- bzw. 
Gleitreibung in der Uhr, das ein kontinuierliches Messen unméglich macht. 

Die angefiihrten Nachteile vermeidet man, 
wenn man zur relativen Abstandmessung die 
Anderung des Druckabfalls an einer von Luft 
Schlauchtiivle durchstrémten Diise bei Annaherung der Gummi- 
flache an diese Diise benutzt. In Abb. 3 ist das 
Prinzip des MeBverfahrens skizziert. Die Disen D ¢ 
und D,, sind so ausgelegt, daB der gesamte Druck 
praktisch bei D; abfallt, wenn die zu vermessende 
Flache F weit genug von D,, entfernt ist. Wird 
F der Diise D,, genahert, so fallt ein wachsender 
Teil des Druckes beim Durchstrémen des Spaltes 
zwischen D,, und F ab. Die an der Skala ab- 


gelesene Druckdifferenz ist also ein MaB fiir die 


Teilkreis 
7Teilstrich=10 “mm SOQ 


Entfernung D,,F. Mit diesem MeBprinzip ist es 
méglich, kontinuierlich mit groBer Ge- 
nauigkeit Abstandsinderungen zu messen. Auf 
die zu vermessende Flache wird dabei praktisch 
keine Kraft ausgeibt. 
Abb. 4. Eichanordnung der MeBdiise (schematisch). Bei den Versuchen wurde eine Diise D,, be- 
nutzt, deren Ende aus einem unten rechtwinklig 
abgeschnittenen Schaft von etwa 20 mm Lange und 5 mm @ besteht, der eine Bohrung von 2 mm & 
hat. Eine Prinzipskizze dieser Diise ist in Abb. 4 zu sehen. Die Diise wurde empirisch geeicht. Im 
benutzten MeBbereich entsprach einer Abstandsinderung von 10-3 mm eine Druckdifferenz von 
5mm W.-S. Da die Messungen spater in den Mefléchern des Versuchsstiicks ausgefiihrt werden 
sollten, wurden Paralleleichungen unter entsprechenden Verhaltnissen durchgefiihrt. Bei dem be- 
nutzten MeSlochdurchmesser von 12 mm konnten keinerlei Abweichungen festgestellt werden. 
Die gesamte Versuchsanordnung zeigt Abb. 5. In der Mitte steht der auf eine starke Tischler- 
platte montierte Versuchsapparat mit aufgesetzter MeBdiise, die in das mittlere MeBloch hinein- 
taucht. (Die Diise wurde in den MeBléchern durch eine nicht gezeigte Vorrichtung zentriert.) Der 
Torsionshebel wird mit Hilfe eines auf derselben Holzplatte befestigten Spannschlosses bewegt, so 
da neben einer gleichmaBigen Drillbewegung auch das Festhalten des tordierten Versuchsstiicks 
in jeder Lage gegeben ist. Rechts in Abb. 5 ist der von einem Elektromotor getriebene Kompressor, 


1 R. S. Rivlin, a.a. O. 
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links die MeBsaule zu sehen, an welcher der Differenzdruck, d. h. die Abstandsanderung der sich in 
dem MeBloch auswélbenden Gummiflache, abgelesen werden kann. 


a pee oe MeBergebnisse sind in Abb. 6 bis 8 zusammengefaBt. In den Diagrammen 
ist fur die einzelnen MeBléc jeweils di 6 uber i i 
eee ee ee eils die Ausw élbung d tiber dem Quadrat des Verdrillungswinkels 

. getragen. (4ur Vereinfachung ist unter @ hier der Gesamtverdrillungswinkel des Gummi- 
e inders Food in Grad verstanden.) Zunachst ist das streng lineare Anwachsen der Verwil- 

ung mit dem Quadrat des Verdrillungswinkels fiir nicl i i i 

Is cht Jerdr f 3 

Tc ger mee 5 mu groBe Verdrillungswinkel (bei 25 mm 

des ylinders Gesamtdrillwinkel bis zu etwa 5°) festzustellen. Abweichungen von der 
Linearitat bei gréBeren Drillwinkeln treten zuerst am Rande des Zylinders auf, in der Mitte des 
Zylinders ist auch bei Drillwinkeln von etwa 9° der lineare Zusammenhang noch gegeben. 


Abb. 5. Gesamtansicht der Versuchsanordnung. 


In Abb. 8 unterscheiden sich die Auswélbungen fiir die beiden Drillrichtungen deutlich. Der 
Unterschied ist wohl darauf zuriickzufiithren, daB die beiden Messingplatten an dieser Stelle nicht 
streng parallel stehen. 

Nimmt man an, daB die Auswélbung proportional der zwecks Auswélbungsverhinderung anzu- 
bringenden Kraft ist, so wird zunachst das Ergebnis der Theorie bestatigt: Man mu, um reine 
Torsion zu erzwingen, an den Endflachen des Zylinders Druckkrafte anbringen, die mit dem Qua- 
drat des Verdrillungswinkels ansteigen.. Die Abweichungen vom linearen Ansteigen bei gréBeren 
Verdrillungswinkeln kénnen einmal darauf zuriickgefiihrt werden, daf sich bei gréBeren Winkeln 
bereits die vernachlassigten Glieder héherer Ordnung bemerkbar machen, zum anderen daf} die 
Voraussetzung, die Verwélbung sei den dort herrschenden Kraften proportional, sicher nur fir 
kleine Verdrillungswinkel — d. h. kleine Auswélbungen — giiltig ist. 

Bei den von R. S. Rivlin’ mitgeteilten Ergebnissen wurden bei etwa gleicher Héhe des Gummi- 
zylinders Gesamtverdrillungswinkel bis zu 25° benutzt. Die mit MeBuhren gewonnenen Werte fiir 
die Auswolbung streuen jedoch zu stark, als daB mit Sicherheit eine systematische Abweichung 


festzustellen ware. 


Dee Sa wins ana 
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Abb. 8. Auswélbung d als Funktion des Quadrates des 


Verdrillungswinkels. 


MeB stelle bei r = 45 mm. 


Abb. 9, Verteilung der Druckkriafte iiber die Endflache 


des Zylinders. 
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Fir die Verteilung der Krafte auf den Endflachen des Zylinders hat sich qualitativ ergeben, daB 
diese Krafte in der Zylinderachse am gréBten sind und nach dem Rande hin abnehmen. Auch in 
der Nahe des Randes wurden noch geringe Druckkrafte festgestellt. Die Gleichung (30) erlaubt es, 
daraus die Folgerung zu ziehen, daf fiir die neue Konstante Y die Ungleichung @/3 < VY < @ gilt. 

Macht man die weitere Annahme, daf der Proportionalititsfaktor zwischen Auswélbung d und 
Kraft unabhangig von der Lage der MeBstelle auf der Endflache des Zylinders ist, so miiBte sich, 
wenn man das Verhiltnis d/q? iiber (a? — r?) auftragt, nach (30) eine Gerade ergeben. (Dabei ist a 
der Zylinderradius, r der Abstand der Mef stelle von der Zylinderachse.) Beschrankt man sich auf 
Gesamtdrillwinkel von w < 5°, so erhalt man die in Abb. 9 gezeichneten Punkte. Es sieht so aus, 
als liefe die Kurve d/q@? als Funktion von (a? — r®), die gestrichelte Gerade, genau durch den Null- 
punkt. Dies wiirde bedeuten, daf® fiir das verwendete Material Y = @ ist, d. h. die durch diesen 
Versuch bestimmte in der linearen Theorie nicht enthaltene neue Elastizitatskonstante Y ware 
gleich dem doppelten Schubmodul in der linearen Theorie. Dieses eigenartige Ergebnis lat sich 
theoretisch weder beweisen noch widerlegen, solange man mit einem elastischen Kontinuum 
rechnet; denn die in der Kontinuumstheorie auftretende Konstante Y wird ja gerade erst durch 
diesen Versuch bestimmt. Da die Versuchseinrichtung lediglich zu Demonstrationszwecken ge- 
dacht war, wurden die Versuche in dieser Richtung nicht fortgesetzt. Ausgehend von einer mole- 
kularen Theorie hat L. R. G. Treloar! eine Form der Formanderungsenergie abgeleitet, die gerade 
zu dem Ergebnis Y = 26 fihren wiirde. Bei den in der Literatur mitgeteilten Messungen an 
Gummi verschiedener Zusammensetzung und verschiedenen Vulkanisationsgrades wurden jedoch 
unterschiedliche Verhaltnisse Y/® gefunden, die im allgemeinen nicht gleich Eins waren, 


10. Zusammenfassung. Am Beispiel der Torsion eines geraden inkompressiblen Kreiszylinders 
wird gezeigt, daB der Henckysche (logarithmische) Verzerrungstensor zur Behandlung von Kérpern 
aus inkompressiblem Material besonders gut geeignet ist. Es wird auf ein MeBverfahren aufmerk- 
sam gemacht, das es gestattet, die Verschiebungen freier elastischer Oberflachen praktisch krafte- 
frei, kontinuierlich und mit groSer Genauigkeit zu messen. Einige mit diesem Verfahren gewonnene 
Ergebnisse werden mitgeteilt. 


(Eingegangen am 22. Marz 1956.) 


Anschrift des Verfassers: Dr. Eberhard Koppe, Gottingen, Bottingerstr. 6-—8 
Max-Planck-Institut fiir Stromungsforschung, Gottingen. 


11, R.G. Treloar, Trans. Faraday Soc. 39 (1943) S. 36 u. 241. 
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Turbulente Vermischung ebener Heifluftstrahlen 
Von W. Szablewski 


1. Einleitung. Auf der Grundlage eines von mir aufgestellten Gleichungssystems zur Berechnung 
turbulenter Strémungen von Gasen stark veranderlicher Dichte habe ich in einer Arbeit’ des Jahres 
1952 das turbulente Strémungsfeld ebener HeiGluftstrahlen, die zunachst durch eine Scheidewand 
getrennt sind und sich dann mischen, unter der Einschrankung berechnet, da® die Strahlen fast 
gleiche Geschwindigkeit haben. In der vorliegenden Arbeit lassen wir diese Einschrankung fallen 
und behandeln das Problem fiir den Unterschallbereich in voller Allgemeinheit. 

Die theoretische Beherrschung des vorliegenden Stroémungsmodells ist in folgender Hinsicht 
von technischer und theoretischer Bedeutung: Betrachten wir die turbulente Ausbreitung eines 
ebenen oder runden heiBen Diisenstrahls im umgebenden Medium, so wird das von uns studierte 
Modell realisiert durch das in unmittelbarer Nahe der Diisenkante befindliche Mischungsfeld. Mit 
der Kenntnis der Modellstrémung wird uns damit Lage und Breite des Winkelraumes der Mischung 
(mit Scheitel in der Diisenkante) im infinitesimalen Ausmaf fiir die technisch wichtigen Diisen- 
strémungen bekannt sein, woraus man gewisse Schliisse auf den weiteren Verlauf dieser Strémungen 
ziehen kann. Theoretisch ist die Kenntnis der Modellstrémung insofern von Belang, als sie uns fiir 
das mathematische Problem der Lésung der partiellen Differentialgleichungen, die das Mischungs- 
feld eines Diisenstrahls beschreiben, die Anfangsprofile (Anfangsbedingungen fiir Abstand Null 
von der Diisenkante) liefert. Diese letztere Aufgabe wird in einer folgenden Arbeit behandelt 
werden. 

Wir méchten noch bemerken, da theoretische Untersuchungen iiber Phanomene der freien 
Turbulenz von Gasen stark veranderlicher Dichte in neuerer Zeit auch von S. I. Pai? durchgefihrt 
wurden. Wir werden darauf in Abschn. 6 zuriickkommen. 


2. Bewegungsgleichungen. Mit Bezug auf meine Arbeit I kénnen wir uns kurz fassen. Ein 
turbulentes ebenes, im zeitlichen Mittel stationares, Mischungsfeld wird unter der Voraussetzung 
konstanten Drucks durch folgende Gleichungen beschrieben: 


00 uu 0ouv O=/ =n OU, =@ 
Impuls: g . = z : 
puls ee a a &(x) ay \o ay | Eu ay)? (1) 
: Gore, Tit. oe. Tv aoc, T 
. it eC Q c 
Energie: i ay == Ei e(x) Raye ’ (2) 
00 u 00 D 00 
Masse: = —— = ae 
Aca ay E (x) Bye (3) 


Hierbei bedeuten x die Hauptstrémungsrichtung, y die Querachse, u bzw. v die Geschwindigkeits- 
komponenten in der x- bzw. y-Richtung, @ die Dichte, T die absolute Temperatur, Cp die spezi- 
fische Warme bei konstantem Druck. E ist eine empirische 
Zahl, die das Ubertragungsverhaltnis von substantiellen 
Eigenschaften und der (Druckschwankungen wahrend des 
Transports ausgesetzten) Geschwindigkeit beim turbulenten 
Austausch angibt. Das Uberstreichen bedeutet zeitliche 
Mittelwertbildung; im folgenden werden wir die Quer- 
striche fortlassen. 
Die scheinbare kinematische Zahigkeit wird durch 

e(x) = x, b,(x) (uy — uy) (uy > uy) (4) 
beschrieben, wo b,(x) die Mischbreite des Geschwindigkeits- 
feldes, x, ein empirischer Koeffizient, uy und u, die Ge- 
schwindigkeiten der ungestérten Strahlen 0 und 1 bedeuten 


(Abb. 1). 
1 W. Szablewski, Ing.-Arch. 20 (1952) S. 67 und S. 73. Die letztgenannte Arbeit wird im folgenden als I 


zitiert. 


®S. I. Pai, J. Aeron. Sci. 16 (1949), S. 463; Quarterly Appl. Mathem. 10 (1952) S. 141. 
’ Grenzschichteinfliisse der Wand werden in dieser Untersuchung nicht beriicksichtigt. 
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Die GroBe c, wird als Konstante angesehen. Es wird die Ubertemperatur 
'— T— T, 
eingefiihrt. Mit der Zustandsgleichung fiir ideale Gase bei konstantem Druck 


o T = konst. 5 
ersetzen wir 0 durch ?: (5) 
__ konst. __ konst. ] 
T. SET aon 


Unter der Annahme von Affinitat der Profile iiber der Mischbreite ergibt sich aus dimensions- 
analytischen Griinden 


b~x, b~x, (6) 


wo 6, die Mischbreite des Temperaturfeldes bedeutet. Wir transformieren auf die Ahnlichkeits- 
koordinate 


My 
und machen mit uo, Jy (= T) — T,) dimensionslos. Fiihren wir noch mit c, = 6,(x)/x die Variable 
1 
(8) 


y2 m1 Cy (Uy — Uy)/Up 


ein, so erhalten wir schlieBlich das Gleichungssystem 


Gaye mut) C= 


du/u du/u u ov (E+ 1) dd/T. 

i . Oey 0 i. 

pews d(a 7)” ldo n P (z oa a 1+ #/T, don | 0, (9) 
ae PHF, , db/O,[2 fu Ov 2 d0/T,| _ 

Energie: d(a n)* ls do 7 E & ol a 1+ 0/T, don | K (10) 

Masse: on aa +- ae = (11) 


Gegeniiber dem Feld konstanter Dichte treten dabei die unterstrichenen Glieder zusatzlich auf. 
Es werde darauf hingewiesen, da Gleichung (11) identisch ist mit 


und damit zum Ausdruck bringt, daB gemaB der Annahme konstanten Drucks und der Vernach- 

lassigung molekularer Wirkungen die Mischung isochor (unter Erhaltung der Volumina) erfolgt. 
Es sind folgende Randbedingungen zu erfillen: 

Geschwindigkeitsverteilung 


u— Uy jl oi =A1ED ie 
ara all fiir ee (12) 
Temperaturverteilung 
o 1 — 
Bae 13 
rtd 0 fiir mean (13) 
SchlieBlich fordern wir noch in Hinsicht auf Anwendungen der Theorie auf Diisenstrémungen 
(2° _,0 fir on>—oo. (14) 
Ug 


(In dem durch die Vermischung noch nicht berihrten Teil des Diisenstrahls, dem Strahlkern, ist 


Verschwinden der Querkomponente der Geschwindigkeit anzunehmen.) es. 
Die letztere Randbedingung bedeutet keine Einschrankung der Allgemeinheit. Denn, wie be- 
reits H. Gértler! fir den Fall konstanter Dichte bemerkt hat, ist mit 


u(y) > Py)» v(m) 
u,=u(n+a), 0,=9+4), v, =v (yn +a)—au(y + a) 


(a beliebige Konstante) eine Lésung des Gleichungssystems. 


auch 


1 H. Gértler, Z. angew. Math. Mech. 22 (1942) S. 244. 
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Die Lésung unseres Gleichungssystems (9) bis (11) ist nun leider nicht explizit beschreibbar 
und erfordert die Anwendung numerischer Methoden. Nur im Spezialfall (up — U,)/Up = O ist die 
Lésung explizit angebbar (vgl. I). ; ; 

In Abschn. 4 werden wir die Lésung mittels eines Iterationsformalismus berechnen. Dieser er- 
fordert sowohl hinsichtlich der Beschrankung des Arbeitsumfanges wie hier auch der Konvergenz 
iiberhaupt die Bereitstellung angemessener Kingangsfunktionen. Diese werden wir in Abschn. 3 
gewinnen, wobei die dort erzielten Ergebnisse tiber den genannten Zweck hinaus eine allgemeine 
Bedeutung besitzen. 


3. Beschreibung des Mischungsfeldes mittels der Fehlerintegralfunktion. Nach vorliegenden 
experimentellen! und theoretischen (vgl. I) Untersuchungen werden die Verteilungen der Geschwin- 
digkeit und Temperatur bei der turbulenten Vermischung ebener Strahlen gestaltlich recht brauch- 
bar durch die Fehlerintegralfunktion beschrieben. Es ist daher naheliegend, fiir die Eingangs- 
funktionen in den Iterationsformalismus des Abschn. 4 die Fehlerintegralfunktion heranzuziehen. 

c 


Dee [ eta. 


as 


Das fiihrt zunachst zu dem die Randbedingungen (12) erfiillenden Ansatz 


on 
hy 1 1 —(on+oa)? copter ' ba 
Uy — Uy oS Ve [e doy = 5} 9 D(C) ae © = Oo” - oa) : (15) 


Die zunachst unbestimmte Konstante oa (Parallelverschiebung der Funktion langs der o 7-Achse 
um den Betrag o a) legt dabei die Lage des Winkelraumes der Mischung zur x-Achse (vgl. Abb. 1) 
fest bzw. die jeweilige, von Geschwindigkeits- und Temperaturdifferenz der Ausgangsstrahlen ab- 
hangige, Drehung des Winkelraumes im (x, ¢ y)-Koordinatensystem nach der einen oder anderen 
Seite hin (vgl. Abb. 4 und 5). Uber den Streckungsfaktor o selbst treffen wir im Rahmen der 
Untersuchung dieses Abschnitts keine Annahme. 
Entsprechend setzen wir fiir die Temperaturverteilung an 

: on Bee ; : 7 1 

—ob (16) 
Der Faktor JE ergibt dabei das Mischbreitenverhaltnis },/b, = jE. (Im Fall (uy — u,)/uUp & 0, 
0)/T, ~ 0 ergibt unser Gleichungssystem (9) bis (11) die Verteilungen 
Lane! 1 oT aan 1 
=—50on)+z5 wd 5 Bat = 


was b,/b, = jE entspricht, siehe I; vgl. ferner Abschn. 4.) 
Fiir die Querkomponente der Geschwindigkeit erhalt man dann aus (11) unter Beachtung der 
Randbedingung (14) 


phage pele 
2 


u— Uy 


Up — Uy 


Up Up 


v Uy — Uy a Ls eee 
a" = (5 +o? +0). (17) 


Zur Bestimmung der Konstanten o a und o b ziehen wir nun die Integralsatze des Impulses und 
der Energie heran. 

Bezeichnen wir fiir den Augenblick die untere geradlinige Begrenzung des Mischfeldes der Ge- 
schwindigkeit (vgl. Abb. 1) mit dj(x) (< 0), die obere mit d,(x) (> 0), so gewinnt man durch Inte- 
gration von (9) unter Beachtung von (11) den Impulsintegralsatz 

d,(x) 
F A Qu (u— uy) dy = Qo Uy (Uy — Uy) a d,(x)). 
Transformiert man auf 7 = y/x, so erhalt man unter Benutzung von o/ey9 = (1 + 9,/T,)/(l + 9/T,) 
die unserem Problem angepaBte Form 
0 ee) 
[ 1+ %6/T, vu u—u Vane ‘1L+%0/T, uw u—uy 
J \it OT, uy ty — ty | 1+ 8/T, uy uy —u, 


dy = 0. (18) 


1 H. Reichardt, VDI-Forschungsheft 414 (1942). 
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Entsprechend gewinnt man den Energieintegralsatz 


0 [o.) 
1+%/T, u & “14+ %,/T, u @ 
1G + O/T, uy q sil | 1+ aT, Up vy Lap (19) 
eee é 


Der Impuls- und Energieintegralsatz 
system und sind hypothesenfrei. 


Gehen wir nun mit den Ansiitzen (15) und (16) in die Integralsitze ein, so erhalten wir fiir oa 
und ob das Gleichungssystem 


gelten natiirlich auch unabhangig von unserem Gleichungs- 


; Ll; 3/2; u u—u , Leer — 
Impuls: — Oe is ee ‘ 15 af 1 Th iD Uy a 
4 | 1+ (Go/Ty) 8/99 Uy uy — uy : | ( 1+ (8/73) 8/9) uy Up — Uy a, ly) 


0a 


Abb. 2, 
es, BY / oa 
A Ufo! Te u il ae o ie 5, #9 ; 
E : Sin OLA a “as a, he ae 
pers hi ai (8o/T,) B/D Ug Do | (1 1 + (9o/T;) 8/9, Up A dC 9 (21) 
5 des ae o 1 eng pean 
mit Tae = 2 2 D(C) > % = yi — a oD) an ; 


Dieses nicht einfache Gleichungssystem haben wir fiir vorgegebene Geschwindigkeits- bzw. Tem- 
peraturdifferenzen (uy — u,)/uy bzw. %,/T, in folgender Weise graphisch gelést: 

Fiir ein vorgegebenes festes o a wurden fiir eine Reihe von o b-Werten die’ Integrale auf der 
linken und rechten Seite z. B. der Gleichung (20) graphisch ermittelt, so daB wir ,,links- bzw. rechts- 
seitige‘‘ Integralkurven tiber der Variablen o b erhalten. Thr Schnittpunkt bestimmt das dem vor- 
gegebenen o a gemaB Gleichung (20) zugehérige o b. Dieselbe Rechnung wurde dann fiir weitere 
vorgegebene o a-Werte durchgefiihrt, so daB wir auf diese Art und Weise die Kurve der zu vorge- 
gebenen go a-Werten gemaB Gleichung (20) zugehérigen o b-Werte erhalten. Der geschilderte 
Rechengang, angewandt auf Gleichung (21), liefert eine zweite Kurve o b iiber o a. Der Schnittpunkt 
beider Kurven ergibt dann die den vorgegebenen Parametern (uy — u,)/up, Po/T, entsprechenden 
Werte o a,o b. 

Unserer Rechnung haben wir die empirische Zahl 


R=? (22) 
zugrundegelegt. Dieser Wert ergibt sich nach den Messungen von 0. Pabst! an runden Heifluft- 


strahlen, auf die wir in Abschnitt 5 zuriickkommen werden. L. Prandtl? gibt fir Phinomene der 
freien Turbulenz bei konstanter Dichte ebenfalls den experimentell ermittelten Wert E = 2 an. 


1Q, Pabst, U. u. M. 8004 (1944) und U. u. M. 8007 (1944). 
2 L. Prandtl, Fiihrer durch die Strémungslehre, S. 114, Braunschweig 1949. 
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da® sich aus von O. Pabst! an einem ebenen Leuchtgasstrahl durch- 
— 1,4 ergibt. Die theoretische Untersuchung des Modells 


Wir weisen jedoch darauf hin, 
gefiihrten Messungen auch der Wert E 
eines ebenen Leuchtgasstrahls soll in einer spateren Arbeit erfolgen. 

Das Ergebnis unserer Rechnungen haben wir in den Abb. 2 und 3 dargestellt. Es zeigt sich, daB 
der EinfluB der Geschwindigkeitsdifferenzen (ug — u,)/Uo und Temperaturdifferenzen Yo/T; sich 


im x, oy = Koordinatensystem in Drehungen des Winkelraumes der Mischung aufsert. 


{oy 


Temperaturfeld 


GS 
——S=—— 


Oy 


Geschwindigkertsteld 


7 L 
Abb. 4. 
4O0Y 

ab Geschwindigkeitsteld 

‘ 

0 & 

10 
aiff 
=7 e 
Abb. 5. 


Man kann folgendes feststellen: Halten wir die Geschwindigkeitsdifferenz (uy — u,)/Up fest und 
lassen die Temperaturdifferenz %)/T, anwachsen, so dreht sich der Winkelraum nach dor aene des 
heiBeren Strahls hin. Dabei ist fir ?)/T, = 0 baw. konstante Dichte die Lage des Winkelraumes im 
allgemeinen nicht symmetrisch zur x-Achse (bzw. der Diisenkante) gelegen. Man entnimmt Abb. 2 
bzw. der hier erhaltlichen expliziten Lésung o a = (1/2 2) (uy — u,)/Up, daB fiir irgendein vorge- 
gebenes (uy — U)/Uy (> 0) die Lage des Winkelraumes durch Drehung aus der cane iis e 
des Falles (wy — u,)/Up ~ 0 nach der Seite des langsameren Strahls hin gewonnen wird. i 

Wir haben diese Verhaltnisse noch geometrisch an Hand der beiden Beispiele (uy — u,)/u, = 1,0 
,/T, variabel und 3/T, = 0, (up — u,)/Up variabel in den Abb. 4 und 5 dareestelle ah ag 


Unsere Ergebnisse sind noch mit dem Streckungsfaktor o behaftet, durch den in 7-Koordi- 
naten die Mischungsbreite, in x, y-Koordinaten der Offnungswinkel des Winkelraumes festgelegt | 
5 || 


10, Pabst, FW-Bericht 09 021 T (1942). 


XXYV. Band 1957 Szablewski: Turbulente Vermischung ebener HeiGluftstrahlen 15 


wird. Das Verhaltnis von AuBen- zu Innenwinkel bzw. X,/x%; (Abb. 6), diese von der 
x-Achse aus zu den Grenzen 0,95 und 0,05 gemessen, bleibt jedoch durch den Streckungsfaktor o 
unberiihrt. Der Streckungsfaktor o selbst kann mit den hier zur Anwendung gekommenen Mitteln 
natirlich nicht bestimmt werden, sondern ist nur mittels einer Hypothese iiber den turbulenten 
Austausch berechenbar. 

Wenn man bedenkt, daB die Grundlage unserer Rechnung die hypothesenfreien Integralsatze 
des Impulses und der Energie bilden; da8 ferner der angenommene, durch die Fehlerintegralfunk- 
tion beschreibbare Verlauf nicht allzuweit von der Wirklichkeit entfernt sein diirfte, ist den hier 
gewonnenen Ergebnissen eine gewisse Allgemeingiiltigkeit nicht abzusprechen. 


Xa 4 Oa 


Oy Xe; 


Geschwindigkeitsfeld Ze lemperaturteld 
=e 


ere 16+ 
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4, Iterative Lésung der Differentialgleichungent. Unser System von Differentialgleichungen 
(9) bis (11) ist explizit nicht lésbar und erfordert zu seiner Behandlung numerische Methoden. 

Naheliegend ist eine Behandlung mittels Stérungsrechnung baw. ein Lésungsansatz mit einer nach 
Potenzen des Stérungsparameters (uy — u,)/Uy fortschreitenden Reihe, wie das H. Gértler? bei der 
Behandlung des Falles konstanter Dichte getan hat. Diesbeziigliche orientierende Rechnungen 
lassen diese Methode jedoch fiir den hier vorliegenden allgemeinen Fall als ungeeignet erscheinen, 
da sie einen sehr groBen Rechenaufwand erfordert und, wenn iiberhaupt, dann jedenfalls nur sehr 
schlecht konvergiert. Wir sind deshalb einen anderen Weg gegangen. 

Wir fiihren der Kiirze halber die Bezeichnungen gy = u/uy, y = 3/0), yp =o v/uy und & = on 
ein. Durch formale Integration erhalt man dann aus (9) und (10) die Gleichungen 


Fe 
F 
Ej1 —2§ (p&—y)aé 
5 0 


ten l ac * 1) : Hea en en (23) 
0 


und 


1 
eres De ites e dé + D,. (24) 
Lop T, iG a 
Die Integrationskonstanten ergeben sich dabei aus den Randbedingungen (12), (13) zu 
00 
C Uy —u, I C ie mtone ts i P 
cena | Uy < ane Ug cons 
J j 
—oo —oo 
oO 
{ 
Bo/T, 1 1) Gee Bo/ Ty : a 
a cae aie Bo/T; i / ‘ 14: Fo/ Ty i 
—oo —0o 


1 Die Durchfiihrung der numerischen Rechnungen lag in den Handen meines bewahrten Mitarbeiters Herrn 


Herbert Moch. 
2 H. Gortler, a. a. O. 
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Fir y erhalt man, eingedenk der Randbedingung (14), 


(25) 


2 we 


jee ee 


Lym, 
én-7. ’ ém-n 
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Dieses System von Gleichungen erméglicht einen Iterationsformalismus, dessen Erfolg (Arbeits- 
aufwand und Konvergenz) natiirlich weitgehend von den Eingangsfunktionen abhangen wird. Als 
Funktionen @o, 7%, Yo benutzen wir die in Abschn. 3 bereitgestellten Funktionen (15) bis (17), wobei 
nunmehr o durch (8) definiert ist. Mit ihnen gehen wir in die rechten Seiten der Gleichungen (23), 
(24) ein. E wird (vgl. Abschn. 3) gleich 2 gesetzt. Auf diese Weise erhalten wir P1> %1- Die Funk- 
tion y, wird dann mittels y, durch (25) bestimmt. Die auftretenden Integrale wurden dabei auf 
graphischem Wege ermittelt. Und so weiter. 
Wir haben nach der geschilderten Methode folgende Falle der Parameter behandelt: 


pes 


1,00 


Up 


0,25 | 0,50 | 0,75 


O/T | 0, | i | 2, | a, 


In den Fallen (uy — u,)/up = 0,25 und 0,50 fiihrte der Rechengang nach zwei bis drei Iterationen 
zum Ziele, d. h. es ergab sich in der graphischen Auftragung eines Din A 4-Blattes ein Stillstand der 
Rechnung. In den Fallen (uy — u,)/uy = 0,75 und besonders 1,0 wurde die Konvergenz jedoch 
wesentlich schlechter. Hier wurden die Eingangsfunktionen durch Extrapolation aus den bereits 
iterativ gewonnenen Verteilungen der Falle (uy — u,)/uy = 0,25 und 0,50 konstruiert. Des weiteren 
erméglichte der alternierende Charakter der Iteration eine Beschleunigung des Rechenganges durch 
Mittelwertbildung. Auf diese Weise kam man auch in diesen Fallen mit vier bis fiinf Iterationen aus. 

In Abb. 7 bis 9 haben wir die berechneten Verteilungen der Geschwindigkeit, der Temperatur 
und der Querkomponente der Geschwindigkeit dargestellt. Man stellt zunichst — in Ubereinstim- 
mung mit Abschn. 3 — aus Abb. 7 und 8 fest, da8 der TemperatureinfluB sich in einer ,,Parallelver- 
schiebung” nach negativen o 7 hin auBert. Die Abweichung von der exakten Parallelverschiebung 
ist nur gering (vgl. hierzu Abb. 11); da der Streckungsfaktor o nach (8) unabhangig von @,/T, ist, 
heiBt das, daB die Mischbreite nur geringfiigig durch die Temperatur beeinfluBt wird. Damit kénnen 
wir fiir physikalische x, y-Koordinaten unser Ergebnis auch so aussprechen: Der Temperatur- 
einfluB auBert sich in einer mit der Temperaturdifferenz wachsenden Drehung 
des Winkelraumes der Mischung zum heif®eren Strahl hin. 

Fiir konstante Dichte bzw. @)/T, = 0 ist dabei der Winkelraum im allgemeinen nicht symme- 
trisch zur x-Achse gelegen. Die sich hier ergebenden Verhaltnisse haben wir in Abb. 10 besonders 
dargestellt, der man entnehmen kann, daf mit wachsender Geschwindigkeitsdifferenz (uy — u,)/Up 
der AuBenwinkel gegeniiber dem Innenwinkel anwichst. 

Auffallend beziiglich der Querkomponente der Geschwindigkeit (vgl. Abb. 9) ist, da® bei Uber- 
schreitung einer gewissen Temperatur- bzw. Dichtedifferenz sich der Zustrom zum Mischungsraum 
in einen Abstrom wandelt. Interpoliert man aus Abb. 9 auf diese Schwelle, so erhalt man als Be- 
dingungsgleichung 

Oo Mora Or Ui 
eine Bedingung, die aus Symmetriegriinden zu erwarten war. (Vgl. hierzu noch Abschn. 7.) 
Die Mischbreiten werden nach (8) durch den Streckungsfaktor 


a 1 


2 x © (Up — 1y)/ Uo 


bzw. 


Z : (26) 


~ 2 ty (© €) (Uo — Uy)/ Uo 


festgelegt. Hier ist x, ein empirischer Koeffizient der Theorie, auf den wir in Abschn. 5 zu sprechen 
kommen werden. Die Breite o c, im o 7-MaBstab ist der Abb. 7 zu entnehmen. Die Festlegung der 
Grenzen ist daher hier mit einer gewissen Willkiir behaftet. Das mit den gut ablesbaren Grenzen 
0,95 und 0,05 erhaltene oc, ist in Abb. 11 wiedergegeben. Eine benachbarte andere Wahl der 
Grenzen ist in der Sache belanglos und bedeutet praktisch nach (26) nur einen anderen Zahlenwert 
des Koeffizienten x,. Man sieht, da8 die Mischbreiten des Geschwindigkeitsfeldes sich nur wenig 
mit der Temperatur andern. Das gleiche trifft fiir die Temperaturfelder zu. 

Im Vergleich mit den Annahmen des Abschn. 3 kann man zeigen, da die hier erhaltenen Ver- 
teilungen der Geschwindigkeit und Temperatur mit wachsendem (uy — U)/Uy gréBer werdende 
Abweichungen von der Feblerintegralfunktion aufweisen, wobei diese mit (9 — u,)/u) marschieren- 
den Abweichungen fiir gréRere 3}/T, gréBere Werte annehmen. Vgl. hierzu Abb. 6, aus der ersicht- 
| lich ist, daB hier relativ gréBere AuBenwinkel auftreten; siche ferner Abb. 10. Fiir das Mischbreiten- 
2 
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verhaltnis hatten wir in Abschn.3 b,/b, = 2 angenommen; hier ergeben sich fiir (tg — U4)/Uo - tee 
0.50: 0,75 mit diesem Parameter wachsende Mischbreitenverhaltnisse in den Grenzen 1,38 bis 1,47, 


wihrend sich fiir (wy — U,)/U) = 1,0 ein mittleres b,/b, = 1,72 ergibt. 


Verteilung der Querkomponente der Geschwinaigkelt 
4 Ov/iug 


Ce, 4 ony) yy 
Ug Uy 
up 


yy 


Ug—-Uy _ L 
eee 


Abb. 9. 


2 0 Z q Rome 7 2 5 
Abb, 10. Abb, 11, 


Die so insbesondere fiir kleinere AuBengeschwindigkeiten (u,/uy < 1) in Erscheinung tretenden 
stirkeren Abweichungen vom durch die Feblerintegralfunktion gelieferten Bild ditrften vermutlich 


bei einer erweiterten Konzeption der Feldgleichungen (1) bis (3) gemildert werden. Nach den expe- 


| 
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rimentellen Untersuchungen, die von S. Corrsin! an Diisenstrahlen in ruhender Luft durchgefiihrt 
wurden, sind am auBeren Rand des Mischungsfeldes starke Schwankungen u! /u zu konstatieren, die 
in der Randzone Anla8 zu nicht mehr — wie hier geschehen — zu vernachlassigenden Druckspan- 
nungen o, geben diirften. Wir gehen jedoch hier nicht weiter darauf ein. : 
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5. Vergleich mit dem Experiment. Zum Vergleich mit der Theorie ziehen wir die Messungen 
yon O. Pabst? heran, die an HeiBluftstrahlen in bewegter Luft durchgefiihrt wurden. Wir beziehen 


1 §. Corrsin, ACR 3 L 23 (1943); S. Corrsin u. M.S. Uberoi, NACA Rep. 998 (1950). 
2 Q. Pabst, a. a. O. 
Dn 
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uns auf die Messungen im sogenannten Kernbereich des Strahls, der durch den Strahlkern gekenn- 
zeichnet ist. Der aus einer Diise mit 10 mm Durchmesser tretende Strahl hatte eine Geschwindig- 
keit von ~ 400 m/sec und eine Temperatur von ~ 400° C. Dieser HeiGluftstrahl trat in einen Frei- 
strahl von 80 mm Durchmesser mit Geschwindigkeiten von ~ 180, 100 und 26 m/sec sowie nor- 
maler Temperatur. Die dimensionslosen Daten betrugen 


Uy — Uy 
Uy 


0,94 | 0,74 053 


O/T, | 1,29 | 1,26 “1,33 1a 1,3 


Abb. 12 und 13 zeigen die von 0. Pabst gemessenen Geschwindigkeits- und Temperaturfelder. 
Es ist zunachst zu sagen, daB die diisennahen Mefschnitte noch deutlich den Einflu8 der Grenz- 
schichtstromung an der Innen- und AuBenwand der Diise einschlieBlich der Warmeabgabe an der 
DiisenauBenwand aufweisen. Es zeigt sich jedoch, daB diese Wandeinfliisse bereits in einem Ab- 
stand von wenigen Diisendurchmessern weitgehend durch den Mechanismus der freien Turbulenz 
beseitigt sind, wie man an der Inzidenz der nachfolgenden MeBschnitte erkennt. 

Da die theoretische Untersuchung unter Vernachlassigung der Grenzschichteinfliisse an der 
Wand durchgefiihrt wurde, kénnen zum Vergleich nur die nicht ganz diisennahen experimentellen 
Profile herangezogen werden, die den Typus der freien Turbulenz zeigen. Da diese letzteren fast 
inzidieren und man annehmen kann, dafs bei Ausschaltung der Grenzschichteinfliisse (bzw. Reali- 
sierung rechteckiger Profile der Geschwindigkeit und Temperatur an der Diisenmiindung) auch die 
sehr diisennahen Schnitte sich diesem Bild einordnen wiirden, erscheint ein Vergleich mit den oben 
erwahnten experimentellen Profilen sinnvoll. 

Es ist jedoch noch folgendes zu bemerken: Auch diese letzteren Profile zeigen noch ein der 
theoretischen Vorstellung nicht gemaBes Verhalten. Entsprechend dem Umstand, dafs die Iso- 
tachen und Isothermen der Diisenstrémung zur x-Achse hin gekriimmt verlaufen miissen, miiBten 
die Profile eine, wenn auch geringe, mit dem Abstand von der Diise fortschreitende Staffelung in 
Richtung der negativen 7-Achse hin aufweisen. Diese Linksstaffelung tritt in den vorliegenden 
Messungen nicht in Erscheinung, was mit den genannten EKinfliissen in Zusammenhang steht, wie 
wir noch auseinandersetzen werden. 

In der Theorie treten zwei empirische Koeffizienten auf: Das Ubertragungsverhaltnis E ergibt 
sich aus den Messungen von O. Pabst! zu EF = 2. Mit diesem Wert von E wird jedenfalls durch die 
Theorie der gréferen Mischbreite des experimentellen Temperaturfeldes im Kernbereich gut Rech- 
nung getragen (vgl. Abb. 13 und 14). Wir wollen noch bemerken, da die Messungen von O. Pabst 
in dem an den Kernbereich anschlieBenden Bereich das Verhaltnis F = 2 mit groBer Scharfe her- 
ausstellen, wie in einer nachfolgenden Arbeit gezeigt werden wird. 

Es verbleibt der Koeffizient ~,, den man als MaB fiir den Mischungsweg ansehen kann. Leider 
ist kaum anzunehmen, daB es sich hier um eine Konstante universellen Charakters handelt. Wahr- 
scheinlicher ist es schon, daB x, den Charakter einer Modellkonstanten hat, ebenso, wie es fiir den 
Koeffizienten x des logarithmischen Geschwindigkeitsgesetzes turbulenter Grenzschichten der Fall 
zu sein scheint. Doch lat sich zur Zeit wenig dazu sagen. Aus dem Vergleich der theoretischen mit 
den vorliegenden experimentellen Verteilungen haben wir uns hier fiir einen mittleren Wert von 
#, entschlossen: 


x, = 0,0082 , (27) 


wobei in den einzelnen Fallen mit modifizierten Werten von x, bessere Ubereinstimmung zu er- 
zielen gewesen wiare®. Mit dem gewahlten Wert von x, ergibt sich jedenfalls eine brauchbare Uber- 
einstimmung der theoretischen Verteilungen mit den MeBreihen von O. Pabst. 

In Abb. 13 haben wir die sich demgema® ergebenden theoretischen Temperaturverteilungen 
eingetragen. Die Streckungsfaktoren ergaben sich dabei nach Formel (26), wobei die Mischbreiten 
o ¢, zwischen den Grenzen 0,95 und 0,05 abgelesen wurden (vgl. Abb. 11), zu o = 18, 26,6 und 41,5 
fiir (uj — u,)/Uy = 0,94, 0,74 und 0,53. Die Ubereinstimmung mit dem Experiment ist sowohl hin- 
sichtlich des mittleren Gradienten bzw. der Mischbreite wie auch hinsichtlich des Verhaltnisses von 
AuBen- zu Innenwinkel als gut zu bezeichnen. 


TO mPabstrasiac®: 


2 Dabei wurde — vgl. Abschn. 4 — das in Formel (26) fur den Streckun c i i 
1 gsfaktor o eingehende o h 
den Grenzen 0,95 und 0,05 abgelesen. Wirde man die Ablesung zwischen den Grenzen 0,99 und 0,01 Behe 
so wurde man praktisch dieselben Verteilungen erhalten, wenn man %, = 0,0053 setzt. : 
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Anders sieht es mit den Geschwindigkeitsverteilungen aus. Ein unmittelbarer Vergleich, wie 
oben, wiirde zwar in bezug auf die mittleren Gradienten noch brauchbare Ubereinstimmung er- 
geben. Doch wiirden die experimentellen gegeniiber den theoretischen Verteilungen um eine mehr 
oder minder groBe Distanz nach negativen 7 hin parallel verschoben erscheinen. Dieser Umstand 
hangt mit den Grenzschichteinfliissen zusammen. Gegeniiber dem theoretischen Modell mit recht- 
eckig angenommenen Profilen an der Diisenmiindung weist die experimentelle Diisenstrémung in- 


folge der Grenzschichteinfliisse an der Diisenmiindung einen Impulsverlustgraben auf. Die experi- 
mentellen Profile vom Typus der freien Tur- 


bulenz sind demgemaB als Profile der effek- 
tiven Impulsstrémung anzusehen. Vom geo- 
metrischen Beobachtungspunkt in der Diisen- 
kante (Abb. 1) aus miissen diese Profile gegen- 
iiber der Theorie zu negativen 7 hin verscho- 
ben erscheinen, wie man sich an Hand eines 
Schemas (Abb. 14) klar machen kann. Damit 
findet gleichzeitig die in Erscheinung tretende 


anfangliche Rechtsstaffelung statt der theoretisch zu erwartenden Linksstaffelung der Profile 
eine Erklarung. Angemessener fiir den Vergleich mit der Theorie erscheint die Wahl eines anderen 
geometrischen Beobachtungspunktes, der dem erwahnten Impulsverlust Rechnung tragt. Das 
kann nach einem Vorschlag von O. Pabst', dem wir uns hier anschlieBen, durch Einfiihrung 
eines sog. effektiven Diisenhalbmessers rj geschehen mittels der Gleichung 


iS trahlachse x 


Abb. 14, 


io) *2 
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Gegeniiber der alten Auftragung der MeBpunkte ttber 77 = (r — ro)/x ergibt sich dann mit ie = 
(r —r%s)/x eine Korrektur in Gestalt einer Parallelverschiebung der einzelnen MeBschnitte, die 
also die Gradienten unberihrt laBt, nach positiven 7 hin um den Betrag? [1 — (r9/r»)|/(«/ro). Die 
Auswertung der experimentellen Verteilungen fiir die Schnitte x/ry = 2 und 6 ergab folgende 
effektive Halbmesser: 


UU, 


| 0,94 | 0,74 | 0,53 


Up 


! : 
Fors | 0,93 | 0,88 | 0,88 


In Abb. 15 haben wir unter Weglassung der diisennachsten Schnitte sowie der nicht mehr dem 
Kernbereich angehérigen Schnitte die auf den effektiven Radius bezogenen experimentellen Ge- 


10. Pabst, a. a. O. . Be 
2 a aae fs technische Anwendungen der Theorie: Die theoretischen Geschwindigkeitsverteilungen 
beschreiben das Feld bei Auftragung iiber 1* = (r — rf)/x. 
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schwindigkeitsverteilungen mit den theoretischen verglichen. Auch hier ist die Ubereinstimmung 
sowohl hinsichtlich des mittleren Gradienten, der durch die obige Koordinatentransformation ja | 
nicht beriihrt wird, wie nunmehr auch hinsichtlich des Verhaltnisses von AuBen- zu Innenwinkel i 
als befriedigend anzusehen. 

Wiinschenswert wire noch ein Vergleich mit Messungen, die den Effekt der Temperaturdifferenz 
herausstellen. Solche Messungen scheinen fiir das hier untersuchte Problem nicht vorzuliegen. Wir 
haben jedoch in Abb. 16 die Theorie mit der Messung (ug — u,)/U = 9,94, O/T, = 1,3 und der von 
H. Reichardt! durchgefiihrten Messung (uy — u,)/U9 = 1,0, Dy/T, = 0 (o = 16,9) verglichen. Die 
Reichardtsche Messung wurde dabei einem Diagramm der gemessenen Staudruckverteilung ent- 
nommen; diese Messung ebenen Charakters erscheint mit nur einem geringen Impulsverlust 
behaftet, der auBer Acht gelassen werden kann. Die geringe Geschwindigkeitsdifferenz der beiden 
Messungen zueinander hat theoretisch gegeniiber der gegenseitigen Temperaturdifferenz nur einen 
geringen HinfluB. Der Vergleich in Abb. 16 zeigt, daB die Theorie den Temperatureffekt an- 
scheinend gut wiedergibt. Eingetragen wurde ferner noch in die Abb. die altere Tollmiensche? 
Verteilungsfunktion®. 


4 W-U, 
Ug 
Fe ties on Messung 1. Reichardt 
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Abb. 16. 


6. Die Theorie von I. Pai. Phanomene der freien Turbulenz von Gasen stark veranderlicher 
Dichte sind in jiingster Zeit auch von S. I. Pai‘ theoretisch behandelt worden. Die von Pai auf- 
gestellten Feldgleichungen sind fundamental von unseren Gleichungen (1) bis (3) verschieden. Wir 
weisen hier nur auf die Unterschiede hin: 


1. Die Impulsgleichung lautet bei Pai 


A(guu) , Aguv) a @ du 
Diese Gleichung gewinnt Pai aus dem Ansatz 
ea ee 
m2 (pute); 


wahrend wir (vgl. I) unsere Gleichung (1) aus dem Ansatz 


or =z (-8w Sa v’') 


dy 
erhielten. 


2. Fir die Temperaturverteilung setzt Pai in Analogie zu laminaren Strémungen mit der 
Prandtlzahl Eins 
ae 


2 Cp 


T=A-+ Bu— (A, B Konstanten). (IT) 


DH wINeichanat.amaaO): 
* W. Tollmien, Z. angew. Math. Mech. 6 (1926) S. 468. 
8 Ohne Anbringung einer Korrektur fiir den Impulsverlust im Falle (ug — u)/uy = 0,94, 3)/T, = 1,3 


wirde sich ein noch groBerer Temperatureffekt ergeben. 
Se Hh Raises A, O), 
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Diese Verteilung entspricht der Energiegleichung 


A(0 cp T u) , Aoep Tv CalmrO(Gn el: du / 
ee ey 


Diese Gleichung enthalt also ein ,,dissipatives“ Glied. Fiir den turbulenten Energieaustausch wird 
derselbe Transport wie fiir den turbulenten Impulsaustausch angenommen, vel. (I). (In unserer 


Notierung also E = 1.) Demnach sind gleiche Mischbreiten des Geschwindigkeits- und Temperatur- 
feldes baw. b,/b,; = 1 zu erwarten.1 


3. Die Kontinuitatsgleichung lautet bei Pai 


IG u) AQ v) ape (IIL) 


Ox dy 


Die Strémungen verlaufen im Feld konstanten Drucks hier nicht isochor: Andernfalls miiBte mit 
o@ T = konst. die linke Seite von (II’) identisch Null zu setzen sein, was nicht sein kann. 


7. Nachpriifung der Annahme p = konst. Unserer Rechnung lag die Annahme konstanten 
Drucks im Mischungsfeld zugrunde. Wir wollen diese Annahme in bekannter Weise? an Hand der 
zweiten Bewegungsgleichung nachpriifen. 


Aus der Bewegungsgleichung in der Querrichtung 


Aovu) , Aqvv) _ _— op , at , Aoy (29) 
Ox dy dy | ax | ay’ $ 
du 00 
wo t= es) (05 + Bus) 


und o, die scheinbare Druckspannung in der y-Richtung ist, folgt (vgl. 1) mit den Abkiirzungen 
y =ulUuy, y = 8/9), p=ovrvjuy und £=o7 die Gleichung 


o 1 , 1 dp 
A (09/2) ug 1 + G/T, dé 
< 1 dp ea o/T, Dieses os eel 1 dog 
“T+ @lton & °F —Y— 3 TE eITomR eP—Y) +e @p ya IF OT, a 
(30) 
Integration iiber € zwischen den Grenzen — ov, oo ergibt 
o* 1 1 ( ) 
4 (0/2) ug 1+ 8,/T, Pe Ps 
= ie E dp/d& 1% a Po: gE. 31 
He -@JTox ?§ —YE—s x, | THe RE —Y)e mo 
Explizite Auswertung ist in zwei Grenzfallen méglich: 
1. Fir 9,/T, = 0 erhalt man die Formel 
(Poo — P—oo) 9 | i at (32) § 
(20/2) us Oki 
2. Fiir (uy — u,)/uy 0 erhalt man (vel. I) 
(Poo — P—co) ou =2 3, : (33) 


(0/2) up T, 


1 Fiir das hier diskutierte Modell ergibt sich aus (I1) 


a (Up — U4)? u — Uy 


BS Pr 2 cp Oo g(l—), wo deme 5 ae 


2 W. Tollmien, a. a. O. 
3 Vel. H. Gortler, a. a. O. 
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Fiir die numerische Auswertung der Integrale im allgemeinen Fall mit den in Abschn. 4 berechneten 
Funktionen 9, 7, y formt man die Integrale zweckmafigerweise um+ 

o 1 1 

4 (0/2) ug 1 + Bo/Ty 


0 co 
gem - 1 a y [th 9 , ai 
ae | [rede trey # | peerme (a2) | Bas (34) 


—oo 0 


(Poo — P—ce) 


Das Ergebnis der numerischen Auswertung von Formel (34) mit den Funktionen des Abschn. 4 
zeigt Abb. 17. 


ee aeores 


(0/2) up 


Abb. 17. 


Wenn man beachtet, da o? eine groBe Zahl ist (vgl. Abschn. 5), so ergibt sich eine sehr kleine 
Druckdifferenz im Mischungsfeld. Thre Vernachlassigung erscheint also gerechtfertigt. 

Bemerkenswert ist noch folgender Sachverhalt: Fir die Temperaturdifferenz Null bzw. kon- 
stante Dichte haben wir nach (32) bekanntlich Uberdruck im Strahl gréBerer Geschwindigkeit 


1 Unter Beachtung der Kontinuitatsgleichung (25) erhalt man zunachst 


co 


‘| : | 
fo) 4 2 
af 1 + (o/ 71) 41 € ee 


d& 2 Z 


il yr 
We Ps) 9 Tt Of TOF 


—oo 


das Integral formen wir um in 
0 


s d | (a 1 | 
| 1 + (8/T1) % 1 + (8/T;) ee 
d§ 2 


dé + 


ee 
| 1 + (8/T) x Uy) | &? 
—co 0 d& 2 


0 


ge 1 ye: P il g2 
=| |; + (Po/T) % 1 + (8/T) eau 1 ot (Pol Fa) yin Lea aT 2 


—0o. 


—0o 


= eae ea ee 
: + (8o/T1) % Up Ve Cal Lae Uo 2 Io 
6 
Beachtet man, daB lim e = : und lim i a be : sowie da} 
poe cc Ley Lay 1 + (8/T;) Peg a ota) % Uo)’ 
dy 0 (e— 5°) fur E<—l 
= zaieil Et kW 
BN Maer r=) fir £S 1, 


so folgt (34). 
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(Diisenstrahl). Bei Anwachsen der Temperatur dieses Strahls gegentiber dem anderen Strahl fallt 
dieser Uberdruck, um schlieBlich in Unterdruck tiberzugehen. Der Ubergang ist nach Abb. 17 
bestimmt durch die Bedingungsgleichung 


Qo Uy = 01 Uy. 


Das ist dieselbe Grenze, bei welcher der Zustrom zum Mischungsfeld sich in einen Abstrom wandelt 


(vgl. Abschn. IV). 


8. Zusammenfassung. Auf der Grundlage eines vom Verfasser in einer vorangegangenen Arbeit 
aufgestellten Gleichungssystem fiir turbulente Strémungen von Gasen stark verdnderlicher Dichte 
wird die turbulente Vermischung ebener Heifluftstrahlen, die zunachst durch eine Scheidewand 
getrennt sind und sich dann in einem Winkelraum mischen, berechnet. Die Lésung des nicht- 
linearen Gleichungssystems erfolgt mittels eines Iterationsformalismus. Berechnet werden die 
Verteilungen der Geschwindigkeit und der Temperatur fiir eine Reihe von Parametern der Ge- 
schwindigkeits- und Temperaturdifferenz. Es resultiert, daB der Winkelraum der Mischung sich 
mit wachsender Temperaturdifferenz zum heiBeren Strahl hin dreht. Weiterhin ergibt sich, dab 
der Zustrom zum Mischungsraum sich beim Uberschreiten der Schwelle gleicher Impulsstrémungen 
beider Strahlen in einen Abstrom wandelt. Der Vergleich mit vorliegenden Messungen an Heifbluft- 
strahlen in bewegter Luft ergibt fiir die Temperaturverteilungen gute Ubereinstimmung sowohl 
hinsichtlich der Mischbreite wie auch des Verhaltnisses von AuBen- zu Innenwinkel. Hinsichtlich 
der letztgenannten Charakteristik ist fiir die Geschwindigkeitsverteilungen befriedigende Uberein- 
stimmung erst bei Einfiihrung eines effektiven Diisenradius, der dem Impulsverlust der Strémung 
an der Diisenwandung Rechnung tragt, zu erzielen. 


(Eingegangen am 24. April 1956) 


Anschrift des Verfassers: Dr. W. Szablewski, Berlin W 8, Jagerstr. 22—23, Deutsche Akademie der Wissen- 
schaften, Forschungsinstitut fur Mathematik. 
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Ein Beitrag zur Plastizitatstheorie * ** 
Von H. Jung 


Herrn Prof. Dr. Siebel zum 65. Geburtstage gewidmet 


1. Einleitung. Wird ein homogener isotroper Kérper durch Oberflachenkrafte belastet, so 
deformiert er sich. Entlasten wir den Kéorper, dann ist ein Teil der Deformation reversibel, wahrend 
der Rest irreversibel ist. Der Kérper wurde durch das an ihm angreifende Kraftesystem elastisch- 
plastisch deformiert. Es soll gezeigt werden, wie derartig elastisch-plastische Zustande einfach 
dargestellt werden kénnen, wenn man zur Darstellung den durch die Hauptspannungsrichtungen 
gegebenen dreidimensionalen Raum heranzieht. 

Der Spannungstensor wird in Ziff. 2 und der Dehnungstensor in Ziff. 3 behandelt. Ziff. 4 gibt 
das Stoffgesetz und seine Darstellung im Raum der Hauptspannungsrichtungen. Die Minimal- 
bedingungen der Plastizitatstheorie werden in Ziff. 5 gegeben. 

Eine Literaturzusammenstellung tiber den derzeitigen Stand der Untersuchungen elastisch- 
plastischer Zustande findet man in dem Buch von R. Hill. 


2. Der Spannungstensor. Wir geben in einem Punkt A des Kérpers ein kartesisches Koordinaten- 
system vor und schneiden zum Zeitpunkt t aus dem Kérper ein Volumenelement dx, dy, dz heraus. 
Die an diesem Volumenelement angreifenden Spannungen sind durch Abb. 1 gegeben. Sie sind die 
Komponenten eines symmetrischen Tensors 


yx 
tee Ox y Ovy Ony ) (61; = Oj: fiir i+ j)- (1) 
Ox x Oy Oz x 
Der Tensor (1) hat die Invarianten 
P,= | | 
0 3 (Ge I Oyy ! Ozx)s 
One} 0, 0, (Ore 
ee yx ye Cay an One 
in P, <4 ale ae ? 
Ox y Oyy yz Oz Ox2 Oxx (2) 
On x Oy x Oxx 
Po = |Oxy Oyy Oxy 
Abb. 1. Volumelement mit Spannungen. Oxx Oy x Oxy 


Die Komponenten des Tensors miissen die Gleichgewichtsbedingungen erfiillen, wenn Massen und 
Volumenkrafte vernachlassigbar klein sind, 


Gx” , Wxy 00x z | 
ae dy eae 
OGy x doy y Oy 2 00; j 
4 OS ay rae Ve ( oie ae 0) ‘ (3) 
00x x 00zy 00x z 
= = ==()). 
Ox dy Oz J 


Ks laBt sich im Punkt A ein neues gedrehtes Koordinatensystem so angeben, daB auf den Flichen 
eines mit diesem Koordinatensystem gebildeten Volumelements nur Normalspannungen 06), 62, 03 
angreifen. Diese Normalspannungen nennt man Hauptspannungen und ihre Richtungen Haupt- 
spannungsrichtungen. Die Richtungscosinusse des neuen Systems gegen das (x, y, z)-Koordinaten- 
system bestimmen sich als Wurzeln der Sakulargleichung 

%— 3 Py) —yv P,— P, = 0. (4) 


* Vortrag auf der Arbeitstagung des VDI-FachausschuB Verfahrenstechnik(Rheologie) in Braunschweig 1956. 


** Vorliegende Untersuchung wurde mit Unterstiitzung der Badischen Anilin- u. Sodafabrik durchoefihrt 
1 R. Hill, The Mathematical Theory of Plasticity, Oxford 1950. “arte Se eae 
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Verwenden wir die Hauptspannungen zur Beschreibung des Spannungszustandes im Punkt A, so ist 


6; .0 20 
ie HOR cope Ot (la) 
OO a, 


Verwenden wir das durch die Hauptspannungsrichtungen gegebene Koordinatensystem zur Dar- 
stellung des Spannungstensors, so stellt sich in diesem System der Spannungstensor als Vektor dar. 


Den Tensor (1) spalten wir in einen Kugeltensor 8 und einen Deviator $$ auf, wobei 


1 
Cy: (Ceo Crue Ora) 


bedeutet: 
dg 0-0 Ing 0, Coe One ! 
P= P+ R=(0 o O}]4 Ons Oe Oy or. 
OF 30:56; One oe O,—= Oy { (5) 
% 9 0 01 — 0% 0 0 | 
=|0 o 0 }+ 0 02 — 0 \ | 
OF 0 G5, 0 0 03 — 09 ) | 


Im Raum der Hauptspannungsrichtungen bildet sich der Kugeltensor auf einer Geraden ay ab, die 
beziglich der Hauptspannungsrichtungen dieselben Richtungscosinusse besitzt. Der Deviator bildet 
sich als Vektor in einer durch den Ursprung des Koordinaten- 
systems gehenden zur Geraden a, senkrechten Ebene ab. 
Abb. 2 zeigt die vektorielle Darstellung von Tensor, Kugeltensor 
und Deviator im Koordinatensystem der Hauptspannungen. 

Durch die hier gezeigte Darstellung wird jedem Punkt des 
physikalischen Raumes (x, y, z) eindeutig ein Punkt des Haupt- 
spannungsraumes (0,, 0, 03) zugeordnet. 

3. Die Dehnungen und Dehngeschwindigkeiten. Die an- 
greifenden Oberflachenkrafte deformieren den Kérper. Zur 
Zeit t habe der Punkt Q die Koordinaten x, y, z. Die Koordi- 
naten des Punktes Q zur Zeit t — dt seien x,y, z. Bezeichnen 
wir mit u; die Komponenten des Verschiebungsvektors 1, so ist 

t=*£—u,, YSHy—u,- 2 =s—U,. (6) 
Ein vom Punkt Q ausgehendes Linienelement ds habe zur Zeit 
t—dt die Langeds. Die Dehnungen in der Zeit dt werden 
definiert durch 


<_ Abb. 2. Der Spannungstensor 
ds 1 ‘d?s — ds im Hauptspannungsraum. 
= In = = ] | = ome +|- 1 . (7) 


Entwickeln wir in (7) den natiirlichen Logarithmus in eine Potenzreihe nach (ds — ds?)/ds®, so 
erhalten wir mit ds = dx; dy; dz die Dehnungen und Verzerrungen 


_ 0Ux 1 [/dux\2 A duy\2 i uz 2 
xx ~ Oy 2| ax ax ax} |? 


duy OUx | OUx duy ; duy : uz : os 


(8) 


Aus (8) erhalten wir die Dehnungsgeschwindigkeiten und Verzerrungsgeschwindigkeiten 
d =< 06; 7 OE; j 0€ij 08; j 9 
ee ate iy ay a (9) 


In Gleichung (9) sind die v; (i, J =, Y, z) die Komponenten des Geschwindigkeitsvektors. ) im 
Punkt Q. Vernachlassigen wir in (8) die quadratischen Glieder, so werden die Dehnungen 


ik ykohns Ou j 10 
erat 3 5) Gy 
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Bezichung (10) in (9) eingefithrt, ergibt in ausfithrlicher Schreibweise z. B. 


— d (Oux\ ux dx | Cux dy ! ux, dz ux (11) 
cox = Fil ae = ox? di! dxdy dt | dx dz dt ' Ox dt’ 
Es ist 
dx . dy _ Wy 
dt SNS) Ch Sat aed i LO 
Put 


In Gleichung (11) sind die Glieder = aj in den Ubergangsgebieten vom plastischen zum elastischen 


Verhalten des Materials nicht mehr vernachlassigbar klein. 

Aus der Theorie der ideal-plastischen ebenen Deformation ist bekannt, daB entlang der Trenn- 
flachen, die das elastische vom plastischen Gebiet trennen, die Normalkomponente der Defor- 
mationsgeschwindigkeit stetig ist, waihrend die Tangentialkomponente unstetig tibergeht. Da der- 
artige Unstetigkeitsflichen in Wirklichkeit eine endliche Ausdehnung haben, so ist in diesen Uber- 


: Oui s 5 Oi 
gangsgebieten Bay tion derselben GréB®enordnung wie aa ° 
Sind die im plastischen Gebiet auftretenden Deformationen von derselben GréBenordnung, wie 
: : ‘ s oe are Pui Pui 
die elastischen Deformationen, so diirfen in vielen Fallen die Glieder et v; gegen ae vernach- 
lassigt werden. 
Die Dehnungsgeschwindigkeiten werden damit 
1 /dv; Ov; 
Meson ee (12) 
Die durch (8) definierten Dehnungen sind die Komponenten des Tensors 
Ex x Eyx Ez x 
O=(ey ly, Exsy | (G7 = ej; fir 17). (13) 
Ex x Ey 2 €xz 
Durch (9) wird der Dehngeschwindigkeitstensor eingefiihrt: 
ex, x Cy x e, x 
GF =(e,, eyy ey | (43 = &; fit U7). (14) 
Ex x Cyz Czy 
Aus (13) erhalt man die Volumenanderung je Volumeinheit zu 
O a Ex x = Ey y Sia Ezz S (15) 
Mit (15) wird 
dO 
We = Cx =f Cvy ae Cruz (16) 
Setzen wir (9) in (16) ein, dann wird in ausfiihrlicher Schreibweise 
dO a) 0 a) 
dt ax (Exe 4 yg, b Ex) Pa 4 oy (xx 4 Eyy Ezs) Vy Bina (Ox ae Eyy ey ay 
0 
== OL (ex ag Evy Tr =) © J (17) 


4. Das Stoffgesetz. Die am Kérper auftretenden Deformationen lassen sich in einen elastischen 
Anteil und einen plastischen Anteil aufspalten. Damit ist 


u; = uj! + uj!. (17) 


Aus (17) ergibt sich, dafS auch der Dehnungstensor und Dehngeschwindigkeitstensor in einen elasti- 
schen und plastischen Anteil aufgespalten werden kann: 


6 =CH 4 SF . 


S* == G*#! a G*Pl ’ (18) 
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Es wird zunachst angenommen, da auch der Schubmodul G eine Funktion des Ortes ist. In diesem 
Fall muB das Elastizitatsgesetz in der Form angesetzt werden, wenn m die Querdehnungszahl ist, 


* 1 df, m—2 
Sh — 5 (B+ ew). (i?) 


Gleichung (19) verlangt, daB der Dehngeschwindigkeitstensor koaxial zum Spannungszuwachs- 
tensor ist. Fiir konstantes G laBt sich Gleichung (19) sofort integrieren und geht in die bekannte 
Form des Hookeschen Gesetzes iiber 


the ee m — 2 
oye (x | m + | 0) : Oo) 
Abb. 3 zeigt die geometrische Darstellung des Elastizitatsgesetzes im Hauptspannungsraum. 


Bleibende Deformationen im Kérper treten auf, wenn die Spannungen die Grenzbedingung 
(Plastizitatsbedingung) erfiillen 


(0, — 02)? + (6, — a5)? + (03 — 0)” = 6 K?, 


GE! 


(20) 


(Ge Gane (855 Oz ? (or Ge) fF 6 (Oy a Ose aE a: :) = 6 K?. 


Im Hauptspannungsraum stellt die Bedingung (20) 
einen Kreis mit dem Halbmesser K dar, dessen Mittel- 
punkt auf der Geraden a, liegt. Die Gerade a, ist Lot 
auf die Kreisebene. 

Ein Korper wird ideal elastisch-plastisch genannt, 
wenn die Spannungen durch die Grenzbedingung (20) 
beschrankt sind. Im Hauptspannungsraum wandert 
der Endpunkt des ,,Spannungsvektors“ auf der Ober- 
flache eines Kreiszylinders mit der Achse ay und dem 
Halbmesser K, wenn elastisch-plastische Deformatio- 
nen auftreten. Rein elastische Deformationen treten 
auf, wenn der Endpunkt des ,,Spannungsvektors“ in 
das Innere des Zylinders wandert, wobei der absolute 
Betrag des Vektors zu- oder abnehmen kann. 

Im Innern eines ideal-plastischen Kérpers kénnen 
sich bei Deformationsvorgingen elastische Gebiete 
aufbauen, die den Deformationsmechanismus wesent- 
lich beeinflussen!. Ein derartiges elastisches Gebiet 
konnte bei der theoretischen Untersuchung des Walz- Abb, 3. Das Elastizitiatsgesetz im Hauptspannungsraum. 
vorgangs nachgewiesen werden. 

Verfestigt sich das Material wahrend der Deformation, dann gilt, wenn die Spannungen die 
Grenzbedingung iiberschreiten, 


= [(, — 92)? + (2 —95)* J 6. ay" = F( loess u). (21) 


V2 


0 


Gleichung (21) sagt aus, daB der Kreisradius K eine Funktion der Belastungsvorgeschichte wird. 
Im Hauptspannungsraum stellt die Gleichung eine Drehflache mit der Achse a) dar, wenn die Funk- 
tion W> definiert ist. Nach R. Schmidt setzt man 
dWp * Pl 4 
ee pe Ce ie (22) 
Bestimmen wir von dem plastischen Dehngeschwindigkeitstensor die Hauptrichtungen und lassen 
diese Hauptrichtungen mit den Hauptspannungsrichtungen zusammenfallen, so stellt auch G*?! 
im Hauptspannungsraum ein Vektor dar. Gleichung (22) ist dann das skalare Produkt von ,,Span- 
nungsvektor‘ und plastischem ,,Dehngeschwindigkeitsvektor“ (skalares Tensorprodukt). 
Aus (21) ergibt sich, daB der Halbmesser K eine Funktion der nicht riickgewinnbaren Defor- 
mationsarbeit ist. 


1H. Jung, Uber die Konstruktion von Gleitliniennetzen; Vortrag auf der Tagung der Deutschen Rheolo- 
gischen Gesellschaft 1956, Oynhausen. 
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Fiir den plastischen Dehngeschwindigkeitstensor wird der Ansatz gemacht 


PL me 2b hare 7 
pis =o furs 9) (23) 
@*Pl — jh 3 : 


Der Ansatz (23) erfiillt die Bedingung, daB bei einer plastischen Deformation keine Volumanderungen 
auftreten. Das elastisch-plastische Stoffverhalten laBt sich mit (19) und (23) beschreiben durch 


0; 1 d{y , m—2 a 

; = 7, {B+ 2a aie (24) 
Der Parameter J ist so zu bestimmen, daB Glei- 
chung (21) erfiillt wird. Abb. 4 zeigt die geo- 


metrische Darstellung des Stoffverhaltens im 
Hauptspannungsraum. 


5. Die Minimalbedingungen der Plastizitats- 
theorie. Ebenso wie in der Elastizitatstheorie 
lassen sich auch in der Plastizitatstheorie Minimal- 
bedingungen angeben.! Diese Minimalbedingungen 
sollen im folgenden erweitert und in vereinfachter 
Form hergeleitet werden. Mit der Erweiterung 
des Gaussschen Integralsatzes von Ostrogradsky ist 


[Bx G*dV=f[T,v, 40. (25) 
Vv 0 


In (25) sind die T; die Oberflachenkrafte und 


die v; die Deformationsgeschwindigkeiten an der 


Abb. 4. Die Plastizitatsflache im Hauptspannungsraum. 


mit (24) ein und erhalten 
[pxorav=[{i, 4|Bx B +22 Bx B] +28 Bay. (26) 
V V 
Wir setzen nun 
G=L+16P, A=A+A (27) 
und verlangen, daB die Variationen 0J§ an den Punkten der Oberflache, wo die Spannungen vorge- 


geben sind, zu Null werden. Gehen wir mit (27) in (26) ein, so wird, wenn die Komponenten yon $f die 
Gleichgewichtsbedingungen befriedigen 


eed 


[log 3 (Bx oP + BSE B xop) +22 Bx dP +200Bx Blav—o. a 


m™m 
V 


Aus der Plastizitatsbedingung und der Verfestigung folgt, da 2 und 6A stets gréBer als Null sein 
mussen. 


Damit eine plastische Deformation auftreten kann, mu mit fortschreitender Deformation / | 


eine stetig zunehmende Funktion sein. Wird 0/ negativ, so findet eine Entlastung statt, die laut 
Voraussetzung rein elastisch ist. Bei einer rein elastischen Deformation ist aber 27 =0. Es kann 
einfach gezeigt werden, da das durch (28) gegebene Extremum ein Minimum ist. Aus (28) folgt 


[8 x 6* dV = Min. (29) 


Da bei der fortschreitenden plastischen Deformation lauter Minima durchlaufen werden, so mu8 in 
Erweiterung der Pragerschen Minimalbedingungen erfiillt sein 


tz 
[ [Bx G* dV dt = Min. (30) 
iV 


' W. Prager u. P. G. Hodge, Theory of perfectly plastic solids, New York 1951. 


Oberflache. In die linke Seite von (25) gehen wir | 
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Die durch (30) gegebene Minimalbedingung eignet sich besonders zur Durchrechnung elastisch- 
plastischer Probleme, wobei aber vorausgesetzt wird, daB die Plastizitatsgrenze (Trennflachen 
zwischen den rein elastischen Gebieten und den elastisch-plastischen Gebieten) bekannt ist. 


6. Zusammenfassung. Ausgehend von der Darstellung des Spannungstensors im Hauptspan- 
nungsraum wird eine anschauliche Darstellung der Plastizitatstheorie gegeben. Insbesonders wird 
gezeigt, daB der Ubergang vom elastischen auf den plastischen Kérper nur méglich ist, wenn bei 
den Dehngeschwindigkeiten nicht lineare Glieder mit beriicksichtigt werden. Eine einfache 


Herleitung der Minimalbedingungen wird gegeben, wobei sich die Minimalbedingungen noch erwei- 
tern lassen. 


(Eingegangen am 28. April 1950.) 
Anschrift des Verfassers: Dozent Dr. H. Jung, Stuttgart N, Robert-Mayer-Str. 66. 
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Direkte Berechnung von Tragfliigelprofilen aus der Druckverteilung 
Von R. Eppler 


1. Einleitung. In einer kiirzlich erschienenen Arbeit des Verfassers! wurde ein auf konformer 
Abbildung beruhendes Verfahren angegeben, mit welchem Tragfliigelprofile aus der Druck- oder 
Geschwindigkeitsverteilung berechnet werden kénnen. Dabei erlaubten eingefiihrte Parameter, 
die Profildicke genau festzulegen und durch Uberlagerung ganze Profilserien zu berechnen. Das 
Verfahren hatte aber noch den Nachteil, daB die berechneten Profile nicht exakt diejenige Geschwin- 
digkeitsverteilung bekamen, mit welcher man in die Rechnung eingegangen war. Diese wurde 
vielmehr nachtraglich mit einem Faktor cos @ multipliziert, wobei @ der von vornherein nicht 
bekannte Steigungswinkel des Profils ist. An den meisten Stellen des Profils hatte dieser Faktor 
zwar keinen groBen EinfluB, so daB man sich entweder mit der entstehenden Geschwindigkeits- 
verteilung begniigen oder mittels eines einzigen Korrekturschrittes eine ausreichende Verbesserung 
erzielen konnte. Nur in der Nahe der Nase, wo cos @ gegen Null geht, war zwischen der vorge- 
gebenen und der wirklich entstehenden Geschwindigkeit keine Ubereinstimmung zu erzielen. Dort 
mute man die Geschwindigkeit im wesentlichen so nehmen, wie es der cos-Faktor angab. 

Dieser Nachteil haftet fast allen zuriickliegenden Arbeiten an?. Lediglich eine Arbeit yon 
M. J. Lighthill® beseitigt ihn, indem sie auf den Logarithmus der Ableitung des komplexen Poten- 
tials iibergeht. Die Allgemeinheit dieses Verfahrens ist aber beschrankt, da nur besonders einfache, 
geschlossen auszuwertende Geschwindigkeitsvorgaben behandelt werden. Im folgenden wird nun 
ebenfalls der genannte Ubergang auf den Logarithmus verwendet, es wird jedoch im Ansatz fiir 
die entstehende Funktion noch eine Singularitat abgespaltet, wodurch eine einfache numerische 
Behandlung méglich wird. Die Vorgabe der Geschwindigkeit unterliegt dann keinen Beschran- 
kungen mehr aufer den selbstverstandlich auch jetzt zu erfiillenden SchlieBungsbedingungen. Diese 
bereiten jedoch keine Schwierigkeiten, sie lassen sich an der Geschwindigkeitsvorgabe leicht nach- 
prifen und erfillen. 

Die Vorgabe der Geschwindigkeit bei stiickweise verschiedenem Anstellwinkel, wie sie in der 
vorangehenden Arbeit angewandt wurde, und wie sie auch bei Lighthill vorhanden war, ist auch — 
jetzt wieder méglich. Sie zeitigt bei dem vorliegenden allgemeinen Verfahren besondere Vorteile. — 
Es werden schon bei der Behandlung der bisher bekannten Aufgaben des Profilentwurfs groBe Vor-— 
teile erzielt und dariiber hinaus noch einige neue, praktisch wichtige Aufgaben ebenso vorteil- — 
haft gelést. | 

Der numerische Aufwand wird beim einzelnen Profil etwas kleiner als beim ersten Verfahren. 
Insbesondere ist bei allen praktischen Beispielen keine Korrektur oder Iteration notwendig. Die — 
nachtragliche Festlegung und Veranderung der Profildicke ist jedoch nicht mehr méglich. Das alte 
Verfahren ist also dann angebracht, wenn man ganze Profilserien mit genau festgelegter Dicke 
berechnen will, bei welchen die Geschwindigkeit in der Nahe der Fliigelnase ziemlich unkontrolliert 
gegen Null gehen darf. Das neue ist am Platz, wenn man Einzelprofile berechnen will, an deren | 
Geschwindigkeitsverteilungen lings des ganzen Profils einschlieBlich Nase hohe Genauigkeits- 
anforderungen gestellt werden. 

Man ist vielleicht geneigt, der kleinen Umgebung der Nase nicht viel Einflu8 zuzuschreiben, | 
solange dort keine unbeabsichtigten Saugspitzen entstehen, und man zieht deshalb vielleicht das | 
alte Verfahren allgemein vor. Es hat sich bei ihm jedoch gezeigt, da bei dickeren und starker | 
gewélbten Profilen die Umgebung der Nase nicht geniigend klein bleibt, in der der cos-Faktor groBen 
EinfluB hat. Wir werden hier auferdem zeigen, daB auch die sorgfaltige Behandlung kleiner Be- | 
reiche der Nasenumgebung gréere Vorteile bringt als bisher allgemein angenommen. | 


1 R, Eppler, Ing.-Archiv. 23 (1955) S. 436. 

2 Vel. auch die Literaturhinweise der Arbeit von FuBnote 1. 

° M. J. Lighthill, A New Method of Twodimensional Aerodynamic Design, Aeron. Res. Council, Rep. and | 
Mem. No 2112, London 1945. | 

* Die Arbeit von Lighthill, FuBnote 3 von S. 32, war mir leider bei der Veroffentlichung der vorangehenden | 
Arbeit nicht bekannt, soda ich dort falschlich die angewandte Art der Geschwindigkeitsvorgabe fiir neu hielt. | 


Inzwischen hat mich Herr Professor Dr.-Ing. A. Weise, Stuttgart, dankenswerterweise auf die Arbeit von. 
Lighthill hingewiesen. } 


| 
| 
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2. Der allgemeine Ansatz. Auch das vorliegende Verfahren beruht auf der Methode der kon- 
formen Abbildung. Wir gehen aus von einer in der ¢-Ebene mit 

C=E+in=reP (1) 

gegebenen unendlichen Parallelstrémung um den Einheitskreis (Abb. 1), deren Zirkulation J’ so 


bestimmt ist, daB der hintere Staupunkt nach ¢ = 1 fallt. Diese Strémung wird durch das komplexe 
Potential 


FQ) =@+iP= Glee + oa) 
He 


Mara 


Ing (2) 


mit der Zirkulation 
Ta a Cain o (3) 


definiert. Wir bilden diese Strémung 
durch eine Funktion z(¢) auf die z-Ebene 
mit 


2=x-+1y (4) 
ab. Soll dort als Bild ebenfalls eine un- 
endliche Parallelstrémung vom Anstell- 
winkel ~ um ein beliebiges Profil ent- 
stehen, dann muB «(f) das AuBere des 
¢-Kreises konform auf das AuBere des 
Profils abbilden, wobei € = 1 in die 
Hinterkante des Profils ibergehen muB, 
wenn dort glatter AbfluB entstehen soll, eats 


und es mu gleichzeitig 

=reell (5) 

gelten. Man kann also 2() in der Form = ee 
AQ=Al+ SR" ©) 


(67 0, teell) 


mit auf |¢| => 1 konvergenter Potenzreihe darstellen. Ist die Abbildungsfunktion 2(¢) bekannt, 
dann ist die Geschwindigkeit v in der z-Ebene 


oo) =e, ( 


Abb. 1. Skizze der eingefiihrten mathematischen Ebenen. 


dF 

dF| [dt 
Sha Soar ) 

| dc | 


als Funktion von £ zu berechnen. Wir wollen die Abbildungsfunktion 2(¢) so bestimmen, daf langs 
des dazugehérenden Profils eine vorgeschriebene Geschwindigkeitsverteilung entsteht. Man 
denkt gewéhnlich bei der Vorgabe der Geschwindigkeit an eine tiber x vorgegebene Verteilung, 
obwohl man sie streng genommen iiber der Bogenlange des Profils vorgeben miiBte. Wir haben 
schon in der vorangehenden Arbeit gesehen, da man an Stelle der Vorgabe v(x) ohne Nachteil 
auch eine Funktion v(y) langs des ¢-Einheitskreises vorgeben kann, da man mit 

xR - cos PD (8) 
(1 = Lange des Profils) immer eine Naherung fiir den Zusammenhang x(~) beim Ubergang von v(x) 
auf u(y) zur Verfiigung hat. Beim Profil der ebenen Platte, das aus (6) mit 2 — Bi pas Or 
0, 2, 3,...) entsteht, gilt (8) mit dem Gleichheitszeichen; bei benachbarten Profilen, wie es Fliigel- 
profile immer in gewissem Maf sind, ist (8) eine Naherung. Meist kennt man aus Beispielen sogar 
noch viel bessere Naherungen fiir x(y) als (8). Erweist sich trotzdem einmal bei irgend einem Profil 
der Unterschied zwischen der Ausgangsnaherung und der am Profil entstehenden Funktion x(q) 
als zu groB, dann kann man gerade das berechnete x(g) als verbesserte Naherung fiir den Ubergang 

3 
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von v(x) auf v(p) verwenden und damit von vorn anfangen. Bei der Einfachheit der spater be- 
schriebenen numerischen Rechnung ware dies durchaus méglich. Bei allen bisherigen praktischen 
Rechnungen hat sich diese Verbesserung allerdings als vollig iberfliissig erwiesen. 

Die Aufgabe, zu einer vorgegebenen Verteilung u(y) das Profil zu finden, wurde in der voran- 
gehenden Arbeit des Verfassers ebenso wie in den meisten zuriickliegenden Arbeiten gelést, indem 
man fiir die Funktion z(C) einen Ansatz durchfihrte. Dadurch kam bei der in (7) benétigten Gribe 
|dz/d¢| eine unbequeme Wurzel in die Rechnung. Um dies zu vermeiden, setzen wir nun hier nicht 
z(€) an, sondern wir bilden mit Lighthill* 


m = =In—ln = (9) 
Fiihren wir den Geschwindigkeitsvektor der z-Ebene mit 
De gee (10) 
ein, dann ist bekanntlich 
oa = ve—i0, (11) 
womit nach (9) 
In = =—Inv +10 (12) 


wird. Wir bekommen also direkt —In v als Realteil. Nun ist die Funktion In dz/dF’, wie man leicht 
nachprift, im AuBeren des €-Einheitskreises I] > Lregular. Sie kann berechnet werden, wenn man 
ihren Realteil auf dem Rand kennt?, wenn also die Geschwindigkeit v(y) vorgegeben ist. Ist In dz/dF 
bekannt, dann ist nach (9) und (2) auch dz/dé gegeben und man kann 2(¢) durch Integration finden. 

Wir miissen also fiir Indz/dF einen Ansatz finden, der auf eine im Sinn der Gleichung (6) 
zulassige Funktion 2(¢) fihrt, und der die Funktion In dz/dF aus ihrem Realteil auf dem Rand 
— In v (e'”) einfach zu berechnen erlaubt. Bei diesem Ansatz nehmen wir abweichend von Light- 
hill? auf die am Rande |¢| = 1 auftretenden Singularitaten Riicksicht. Sie riihren einerseits von 
den Nullstellen von dF/d¢ her, denen die Staupunkte am Kreis entsprechen, und die nach der aus 
(2) und (3) folgenden Gleichung 


dF : 1 i . 
are oS G(s) (= teats (13) 


bei € = 1 und bei ¢ = —e?'* liegen. Andererseits hat aber auch dz/dé bei € = 1 eine Nullstelle 
(vgl. Abb. 1), wenn im Bild dieses Punktes die Hinterkante des Profils entstehen soll. Soll die 
Strémung das Profil mit endlicher Geschwindigkeit verlassen, dann verschwindet dF /dz bei € = 1 
nicht, die beiden Singularitaten von dF'/d£ und dz/dé heben sich dort in (9) also auf. Der Hinter- 
kantenwinkel des Profils mu verschwinden. Endlicher Hinterkantenwinkel hatte in (9) bei € = 1 
eine leicht zu tibersehende Singularitat zur Folge. Wir beriicksichtigen sie im folgenden nicht, wir 
geben vielmehr in Ziffer 5d) an, wie man auch im Rahmen des mathematisch verschwindenden 
Hinterkantenwinkels praktisch geniigend vollige Profilhinterseiten erreicht. Wir machen danach 


fiir (9) den Ansatz 


dz a etifall ; — : os 
In =—InC—In le G + | tS (dp, tt bp) Bm" (14) 
m=0 
Er ergibt mit einer auf dem ganzen Rand |f| = 1 konvergenten Potenzreihe, wie wir sie fiir die prak- 
tische Durchfiihrung spater benétigen, bei € = 1 endliches also glatten AbfluB8 und verschwin- 
a 


denden Hinterkantenwinkel. Aus(14) folgt durch Vergleich mit (9) und (13) 
dz 1 ~ : 
In dé == In ( =) > (4, oie U bn) erie 0 (15) 
Diese Funktion hat bei € = 1 eine Singularitat, da dort beim Ubergang von € nach z Winkelver- 


doppelung stattfindet, am vorderen Staupunkt dagegen keine mehr. Dort wird die Abbildung 
konform, das Profil bekommt als Bild des Kreises dort stetige Tangente. 


1 FuBnote 3 von S. 32. 


2 A. Hurwitz u. R. Courant, Funktionentheorie, 2 li 5) 
itn ee 1onentheorie, 2. Auflage, Berlin 1925. 
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Die aus (15) folgende Funktion 


(16) 


mu mit (6) vertraglich sein. Sie ist ebenso wie (6) in |f| > 1 analytisch, da dies fiir die Reihe im 
Exponenten der Fallist. Man kann (16) also nach fallenden Potenzen von ordnen und das Ergebnis 
mit der Ableitung von (6) vergleichen. Dieser Vergleich verlangt je 


ds 1 1 mao (am = ibm) > ii 
dé ee 


ine a 


eae (17) 
a, = 1, (18) 
by (19) 


Die erste Bedingung sorgt dafiir, daB im Unendlichen keine Drehung stattfindet, (18) und (19) 
bringen das Glied €-! zum Verschwinden. Damit ist gleichzeitig gesichert, daB bei der Integration 
von (16) ein geschlossenes Profil entsteht. Wir nennen (18) und (19) deshalb die SchlieBungs- 
bedingungen. 

Wir sorgen noch dafiir, das die Strémungsgeschwindigkeit im Unendlichen gleich der Einheit 


wird. Nach (14) verlangt das 


V = veo) = Ce = 1, (20) 

was wir mittels 
Cle (21) 
a, = 0 . (22) 


erreichen. Danach kénnen wir die Geschwindigkeit v immer mit dem allein wesentlichen Geschwin- 
digkeitsverhaltnis v/V verwechseln, das in der Bernoullischen Gleichung 


Pie? v \2 
foe) == 1\ 
ee, S 


bei der Berechnung des Drucks p auftritt. 
Die gestellte Aufgabe verlangt nun, die noch freien a,,, 6, so zu bestimmen, daB fiir ¢ = e'” die 
nach (14) und (12) berechnete Geschwindigkeit v(p) einen vorgegebenen Verlauf hat. Wir gehen 


dazu mit € = e'? in (14) ein und setzen 


P(¢) = wish Am COS mM YP + ln sin m P 5 (24) 
m=(0 : 

Q(v) = 5) —a, sinm@p + bn cosmo. (25) 
m=0 


Dann ergibt die Trennung von Real- und Imaginarteil 


; | ' : 
—In o(y) + i OG) =—In2 jos & - x] + Py) +1 E + O(p) — {{ah}], 
0 0 a eS : 2 «) | | | 28) 
; ee =P =e +260 

as tay} (#+2a05p 52m). 
Die Richtung @ der Geschwindigkeit springt im Staupunkt, wo die Geschwindigkeit verschwindet, 
um z. Der Realteil von (26) gibt bei bekannten a,,, 5,, und einem Anstellwinkel « die Ge- 
schwindigkeit v(p), die bei diesem « auftritt. Ist andererseits eine Funktion v(y) vorgegeben und 
dazu ein x, bei welchem v(y) auftreten soll, dann kann man wegen 


P(y) = In 2 |cos & — “| — In v(¢) (27) 


sofort durch harmonische Analyse die zugehérigen Koeffizienten dm, bm gewinnen. Die fiir dy, bo, 
a, und b, abgeleiteten Bedingungen tibertragen sich mittels (24) und (27) auf die Vorgabe v(9). 
Wahrend b, in P(g) nicht auftritt, (17) also mit der Vorgabe u(y) nichts zu tun hat, folgt aus (22), 
(18) und (19) mit den Formeln fiir die Fourier-Koeffizienten 


27 : 
[ Py) dp =0, (28) 
0 

3* 
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[Pq cosydyp =2, (29) 


0 

2a 

i P(p) sing dp = 0. (30) 
0 


Hat man zu einem Anstellwinkel « eine Verteilung v(y) gefunden, die (28)—(30) erfiillt, dann kann 
man auch die iibrigen am, 6, bestimmen, mit ihnen in (16) eingehen und durch Integration die 


Abbildung 
ZS Gmtibnlt-™ 
1 i 
2(C) = [1—4}< ‘ dé (31) 


und damit die ganze Strémung der z-Ebene gewinnen. Das Profil dieser Strémung ist das Bild des 
Einheitskreises. Geht man mit ¢ = e’” in (31) ein, dann folgt mit (24), (25) und (27) nach der Zer- 
legung in Real- und Imaginarteil entsprechend (4) 


x) = | 4.sin £ cos (4 x] ie cos (4 | o)] do , (32) 
yip) = | —Asin © loos [# —a) 2 sin (% + OW) ap (33) 


Man muB also aus der vorgegebenen Verteilung v(p) nur die Funktion Q(y) berechnen und be- 
kommt damit durch (32) und (33) das zugehérige Profil in Parameterdarstellung. Dabei ist Q(¢) 
als konjugiert harmonische Funktion von P(p) durch das Integral 


QP) = 5s | P(o) etg 2" do (34) 


definiert?. : 
Der Auftrieb des Profils ist mit der Zirkulation J’ nach (3) wegen des Kutta-Joukowskischen 
Satzes (0 = Dichte des Strémungsmediums) 


A=0 VI = 6-477 sin &, (35) 


und der Auftriebsbeiwert ist dann 


Cai = SP Shige (36) 
0 ys l 
rast 
wenn / die Lange des Profils ist. Bei der ebenen Platte ist bekanntlich 
CH= 20, SLI, (37) 
sie muf also in unseren Formeln mit der Lange 
L=4 (38) 


enthalten sein. Bei Profilen mit nicht verschwindender Dicke ist c, etwas eréBer als nach (37). 
Sie fiihren also auf wenig kleinere J- Werte als (38), und es ist 
l 
q™ IL. (39) 
Fiir « = 0 verschwindet c,, das Profil bekommt also eine zur x-Achse parallele Nullauftriebs- 
richtung. 
Das Moment der Luftkrafte beziiglich des Ursprungs ist nach H. Glauert®, in mathematisch 


negativem Sinn positiv gezahlt, 
1 dF \2 
M =~0 Re p (a Z a (40) 
wo der Weg des Integrals das Profil in mathematisch positivem Sinn umschlieBen mu, Fiir « = 0, 
also fiir verschwindenden Auftrieb folgt daraus das Nullmoment, wenn man mit unserem Ansatz 


* I. E. Garrick, Conformal mapping in aerodynamics, with emphasis on the method of successive conjugates ; 
Proceedings of a Symposium: Construction and Application of Conformal Maps. Washington 1952. 

° W. Weinberger, Luftfahrtforschung 17 (1940) S. 3. 

° H. Glauert, Die Grundlagen der Tragfliigel- und Luftschraubentheorie, Berlin 1929. 
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eingeht, 
My = 0 47by, (41) 
und der Beiwert desselben wird 
eee ey (42) 
sever 
wo 6, als Koeffizient der Fourier-Reihe P(q) 
b, = 4 | Pio) inland (43) 
0 


ist. Fiihren wir den Ubergang von P(g) auf Q(y) mittels (34) ohne die a,,, b», durch, dann muB 6, 
fiir das Moment gesondert berechnet werden. 

Unser Ansatz ist insgesamt gleichbedeutend mit einem Ansatz der Form (6) und ermiglicht es, 
zu einer Geschwindigkeitsverteilung v(q) ein Profil samt Auftrieb und Moment zu berechnen. Es 
entsteht sicher ein geschlossenes Profil, wenn (29) und (30) erfiillt sind. Durch (22) oder (28) wird 
garantiert, dal} v(y) zu der Anstrémgeschwindigkeit V = 1 gehért, was durch Multiplikation von 
v(p) mit einem konstanten Faktor stets erreicht werden kann. Derselbe beeinfluBt (29) und (30) 
nicht. 


3. Die Einfiihrung der Geschwindigkeitsverteilung. Mit den abgeleiteten Formeln ist das eigent- 
liche Problem der Profilberechnung aus der Geschwindigkeitsverteilung im Prinzip gelést. Man 
kann ein x oder wegen (36) ein c, vorgeben, bei welchem eine ebenfalls frei gewahlte Geschwindig- 
keitsverteilung v(p) auftreten soll und daraus das Profil berechnen. Mittels (28) bis (30) kann man 
von vornherein ausprobieren, ob ein geschlossenes Profil entsteht und ob v(p) zu einer Anstrém- 
geschwindigkeit V = 1 gehért. Damit ist aber noch nicht allgemein garantiert, daB die gewahlten 
Vorgaben auch auf ein physikalisch sinnyolles Profil fiihren. Dies ist erst mit Hilfe des numeri- 
schen Ergebnisses nachzupriifen. 

Wir wollen hier nun v(q) in einer giinstigen Form mit geeigneten Parametern einfiihren. Diese 
Form wird sehr allgemeine Profilentwurfsaufgaben in einfacher Weise erfassen. Die Bedingungen 
(28) bis (30) werden mit den Parametern ohne Probieren erfiillt. Die wichtigste Méglichkeit un- 
sinniger Profile wird von vornherein ausgeschlossen. 

Wir haben schon erwahnt, daB wir bei der praktischen numerischen Auswertung der Formeln auf 
die Konvergenz der in (14) enthaltenen Reihe auf dem ganzen Rand angewiesen sind. Dies mui 
durch eine entsprechende Vorgabe der Funktion v(y) erreicht werden. Wir verlangen von ihr, dah 
sie P(y) nach (27) zu einer stiickweise analytischen Funktion macht, die als héchste Singularitaten 
Knicke, also endliche Sprungstellen in der ersten Ableitung hat. Dann werden die am, bm von der 
Ordnung 1/m?1, was numerisch bewaltigt werden kann, und was gleichzeitig geniigend glatte Pro- 
file garantiert. : 

Diese Voraussetzung ist noch sehr allgemein. Die einzige durch sie bedingte wesentliche Ein- 
schrankung ist die Stetigkeit von v(y) auch an der Hinterkante y = 0. Diese Stetigkeit ist jetzt 
nicht mehr zu umgehen, so wie dies bei der Wirbelverteilungstheorie und in gewissem Sinn auch 
bei der vorangehenden Arbeit? méglich war. Die durch eine Unstetigkeit von P(g) entstehende 
Singularitat ware jetzt zu stark. Im iibrigen verlangt unsere Voraussetzung nach (27) beig=2-+ 20 
ein Verschwinden von v(g), also dort den vorderen Staupunkt, und sonst iiberall nicht verschwin- 
dendes u(y), insbesondere auch an der Hinterkante y = 0); es entsteht der erwihnte Hinterkanten- 
winkel 6 = 0. 

Auch im Rahmen dieser Einschrankung kénnen noch unsinnige Profile definiert werden, bei- 
spielsweise solche mit Selbstiiberschneidungen. Wir schliefien sie nicht mehr durch allgemein for- 
mulierte Bedingungen aus. Sie sind praktisch leicht zu vermeiden (vgl. auch Ziff. 5d). 

Bei der endgiiltigen Einfiihrung der Geschwindigkeitsverteilung wollen wir noch einen Vorteil 
unserer Formeln ausniitzen. Bis jetzt sind sie so abgeleitet und so aufzufassen, daB ein fester An- 
stellwinkel « gewahlt wird und dazu eine Verteilung v(¢), die (28) bis (30) erfiillt. Daraus wird mittels 
(32) und (33) ein Profil berechnet, dessen Geschwindigkeitsverteilung v(q, «) wegen (26) bei dem 
gewahlten Anstellwinkel « mit der Vorgabe v(y) iibereinstimmt. Nun stellt die praktische Aero- 


1 H. v. Mangoldt — K. Knopp, Finfiihrung in die héhere Mathematik, Bd. 3, S. 508, Leipzig 1944. 
2 R. Eppler, FuBnote | von S. 32. 
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dynamik aber noch eine etwas allgemeinere Aufgabe. Es ist haufig nicht entscheidend, v(p) bei 
einem festen x oder c, am ganzen Profil vorgeben zu kénnen. Man will vielmehr oft ganze Anstell- 
winkelbereiche mit giinstigen Eigenschaften haben. 

Das wichtigste diesbeziigliche Problem tritt bei der Berechnung von sogenannten Laminatr- 
profilen auf. Dort wird beispielsweise in einem ganzen c,-Bereich zwischen zwei Werten Cay < Cay 
auf Ober- und Unterseite vom vorderen Staupunkt bis zu einer gewissen Stelle hin nicht abnehmende 
Geschwindigkeit verlangt, weil anderenfalls die laminare Grenzschicht instabil wird. Im Idealfall 
ist dann die Geschwindigkeit in dem betreffenden Stiick der Oberseite bei c,,, in demjenigen der _ 
Unterseite bei c,) konstant!. Man muB also nach der bisherigen Auffassung unserer Formeln v(9) 
etwa bei dem c,, entsprechenden Anstellwinkel «, vorgeben und hat dann auf der Oberseite ein 
Stiick weit konstantes v(p) zu nehmen, auf der Unterseite ein veranderliches v(y), aus welchem erst 
mit der durch die Anstellwinkelanderung bis c,) bedingten Geschwindigkeitsanderung zusammen 
ein Stiick konstanter Geschwindigkeit entsteht. 

Dieses gleichzeitige Mitverfolgen eines zweiten c, oder x, das frither viel zusatzliche Mihe 
machte, gelingt hier sehr einfach. Nach (15) bestimmen die a,,, b, die Abbildung 2(), also die 
Form des Profils. Nehmen wir ein festes Profil, also feste a, b» und andern nur in (3) und (2) den 
Anstellwinkel der Ausgangsstrémung in der ¢-Ebene, dann wird dadurch die Abbildung 2(¢) nicht 
beeinfluBt. Es andert sich nur wegen (12) und (14) die Geschwindigkeit in der z-Ebene, also auch 
am Profil. Diese Anderung ist am einfachsten mit (27) zu tibersehen. Dort darf sich P(y) wegen 
(24) bei festen a, b» nicht andern, also gilt fiir zwei verschiedene Anstellwinkel x, und « 


608 & a >| (44) 
Man kann damit an einem durch v(q, x) und « definierten Profil fiir jedes ¢— leicht v(g, «,) berechnen, 
ohne das Profil zu kennen. Umgekehrt kann man auch in einem beliebigen g-Bereich direkt v(g, x) 
wiahlen und daraus mittels (44) ausrechnen, wie man v(g, x) vorgeben mu, damit beim Anstell- | 
winkel x, die dort gewiinschte Geschwindigkeit v(q, x,) auftritt. Geht man mit (44) in (27) ein, dann | 
entsteht die alte Gleichung, nur mit «, anstatt «. : 

Allgemein kann man ebenso wie in der vorangehenden Arbeit? ~« = «*(g) als (praktisch stiick- | 
weise konstante) Funktion von g vorgeben und dazu mit einer Verteilung v*(¢) | 


— In v(¢, «,) + In 2 lees fe 2) = —In v(g, x) -- In2 


cos ir) —In v*(9) (45) | 


| 
bilden. Auch dadurch sind wegen (24) gewisse am, bm und ein konkretes Profil definiert. Dasselbe | 
berechnet sich, wie man leicht verfolgt, aus (32) und (33), wenn man auch dort «*(g) statt « und | 
v*(p) statt v(p) einsetzt, also formal ohne weitere Anderung. Die Geschwindigkeitsverteilung des | 
entstehenden Profils ist bei einem festen « wegen der aus (27) folgenden Bezichung 


P(@) = in 


v(p,.0) = 2 [cos (ga e—P(o) (46) 


cos i. = z| 


cos Le —«a-(p) 


mit (45) 


v(y, «) = v*(g). (47) | 


Sie stimmt also dort mit der Vorgabe v*(q) tiberein, wo das gewahlte « mit der Vorgabe o*(p) tiber- 
einstimmt. Man wablt also mit «*(p) in gewissen p-Bereichen und damit auf gewissen Profilstiicken 
den jeweiligen Anstellwinkel ~ oder Auftriebsbeiwert c,, bei dem die Vorgabe v*(p) dort als Ge- | 
schwindigkeit auftritt. 

Die Funktion P(y) nach (45) darf entsprechend unserer Voraussetzung keine Sprungstelle | 
haben. Springt beispielsweise an einer Stelle gi der Anstellwinkel «*(p) von a; nach O41, dann | 
muB dert wegen 


g 
cos (5 mila Oj | cos a ae xi 


' Vgl. die spateren Abbildungen 4 und 7, oder auch I. H. Abbott, A. E. v. Doenhoff, L. S. Sti 
of Airfoil Data. NACA-Report 824, Washington 1945. nT eee ae 
2 FuBnote 1 von S. 32. Vel. auch FuBnote 4 von S. 32. 


lim [— In v*(p + Aq)] + n2 (48) | 


Ag—>0 


= lim [— In v*(g — Ag)] + In2 


Ag—0 
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oder 


lim peneAG)) 2 


FO v *( — Ag) | end (e a x, 
5 i 


(49) 


auch v*(p) eine Sprungstelle haben. Mit (49) ist gleichzeitig die Stetigkeit von v(~, «) nach (47) bei 
jedem « gesichert. 

Auf diese Weise kann man die erwahnte Aufgabenstellung der Laminarprofile direkt lésen, in- 
dem man im entsprechenden Bereich der Oberseite «*(p) = ox, wahlt, in demjenigen der Unterseite 
«*(~) = ay und jeweils konstantes v*(p) vorschreibt. Zur Behandlung solcher und auch noch allge- 
meinerer Aufgaben machen wir fiir «*(p) und v*(p) einen speziellen Ansatz. Wir fiihren ihn zu- 
nachst formal durch und erfiillen alle bisher abgeleiteten Bedingungen. In Ziffer 4 beschreiben wir 
den Gang der numerischen Rechnung, in Ziffer 5 zeigen wir, was mit dem Ansatz zu leisten ist. Dies 
wird ziemlich viel sein, er ist aber trotzdem speziell und kann wenn noétig leicht durch einen anderen 
ersetzt werden. 

Wir teilen das Intervall (0, 2 z) in sechs Teile ein, fiir deren Grenzen 


P=9< 9, <P. < G9 < G1 <5 <M = 20 6°) 


gelten soll. In den einzelnen Abschnitten soll jeweils 


a*(y) =o, =comst Mik y<y) (51) 
sein. Dazu setzen wir v*(y) in der Form 
v*(y) =vwiilp) GaSe S*#i) (52) 


an, wo die v; und K; freie Parameter sind. In den vier mittleren Abschnitten werden wir mit kon- 
stantem v*(p~) auskommen. Wir kénnen also 

Sho Ko Ke = 0 (53) 
setzen. Die beiden anderen Abschnitte ergeben die Umgebung der Hinterkante. Dort ist bei Pro- 
filen mit brauchbarer Dicke immer ein Geschwindigkeitsabfall vorhanden, iiber dessen Form mit 
den Funktionen w,(g) und w,(~) verfiigt werden kann. 

Bei der Wahl von w,(y) und w,(p) beriicksichtigen wir innerhalb ihrer Definitionsintervalle 
jeweils die von P(y) geforderten Stetigkeitsbedingungen. AuBerdem diirfen diese beiden Funk- 
tionen nirgends verschwinden. Wir miissen dann mit den in (52) freien Gréfen nur noch die 
Stetigkeit von 


cos ($ a; ine Kiinwlo)-ein?, @i1=9Se); @=1,2,-.-6) 64) 
an den Stellen q; erreichen und die Bedingungen (28) bis (30) erfiillen. Wir setzen fiir w,(p) und w,(¢) 
wi(P1) = Wel(Ps) = 1 (55) 


voraus. Dann verlangt die Stetigkeit von P(p) nach (54) an der Stelle yy = 0 bezichungsweise 
Yg =2X 


P(p) =In 


—Inv, — K, In w, (22), (56) 


In |cos x,| — In v, — K, In w,(0) = In [cos ag 


an den anderen Stellen 


cos & x) In v; = In 
Wir iiberdecken mit dem stiickweise.konstanten v*(y) auch dic Fliigelnase y ~ x. In dieser 
Gegend kann P(g) unendlich werden, wenn in einem der Intervalle 


In 


cos (Fats) i cei le ee) ee oT) 
| 


p=at+ 2a; (58) 


wird. Dies mu nach den eingangs dieser Ziffer fiir P(p) getroffenen Voraussetzungen vermieden 
werden. Wir verteilen zweckmaBig unsere p-Abschnitte so, da® drei davon auf die Oberseite, drei 
auf die Unterseite des Profils filhren. Wegen (8) ist dann py, ~ 7. Damit nun (58) in (54) nicht vor- 
kommen kann, verlangen wir 

Xe > Oy (59) 
und 


tt + 204 <3 <4 +24. (60) 
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Mit (56) und (57) haben wir sechs Gleichungen fiir die sechs v;. Legen wir in (50) und (51) alle 
yi; und o; von vornherein fest, dann sind in P(y) nur noch K, und Kg zur Erfiillung von (28) bis (30) 
frei, was nicht ausreicht. Wir miissen noch eine der GréBen 9;, x; freihalten. Dafiir ist ~, am besten 
eeeignet. Durch die Wahl von v*(p) und «*(~) ist wegen (47) fiir jedes « eine ganze Verteilung 
v(@. x) und damit ein Auftrieb gegeben. Derselbe muB mit dem wegen (35) durch « gegebenen tiber- 
einstimmen, was durch (28) bis (30) garantiert wird, denn fiir geschlossene Profile gilt ja gerade (35). 
Wir kénnen also (28) bis (30) dann gut erfiillen, wenn die zusatzlich freigehaltene GréBe den Auf- 
trieb gut beeinflussen laBt. Genau das gelingt mit ~; besonders giinstig, denn nach (57) hangt da- 
yon v,/v, und damit das Geschwindigkeitsverhaltnis zwischen Ober- und Unterseite sehr empfind- 
lich ab. (Im Rahmen von (60) kénnen mit ¢; alle positiven Werte v/v, erreicht werden.) Alle 
anderen qj; und die «; waren weniger geeignet, auberdem werden sie spater zur Einfiihrung von Ent- 
wurfsmerkmalen benétigt. Wir lassen also 3 frei. 

Wir gehen danach mit (54) in (29) und (30) ein. Diese Integrale sind, zunichst bis auf die Anteile 


Pj 
= i cos y In wg) dy = Wi, (61) 
Pi-] 


Pj 
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— i sin p In wy) dp = Wi, (62) 


iI 


iiber die wir ohne die wp) noch nichts aussagen kénnen, geschlossen auszuwerten. Es ist namlich | 


dp = (sing + sin2«;) In 


cos ( >, 
cos b ss » 


[ees g In os (F —ai) 


und 


| sin y In ea G — x 


(63) 


+ 5 sin 2 x; + cos gy) + Const. 
(64) 


dy = — (cosy + cos 24;) In 


Damit wird insgesamt in (29) 


2% Pi 
4 6 COS | 6 ; 
| eos gy P(p) dp = K, W., + K, We, + Da sin 2 «; In |— ( aL = (~: — Vi_-1) cos 2 0; 
0 au cos (eS re 
Laem , : : ; 

oe = (sin p; — sin y;_1) + sin g; |In |cos G » In v, sin Pi—1 Im cos ts — x == hn »,| 

| J j 
(65) 


Hierbei heben sich bei der Summation der vorletzte Term mit i und der letzte mit i + 1 wegen (57) 
jeweils weg. Es bleiben von diesen beiden Termen nur der letzte mit 1 = 1 und der vorletzte mit 


i = 6 iibrig, die wegen der sin-Faktoren verschwinden. Auch die Glieder des dritten Terms fallen 


bei der Summation weg, so dab 


ad 6 eos |! ia: 
[ PO? 08 9 dp = Ky Wer + Ky Weg + >) join 2x; i +> (Gi — Gin) cos Zaz} (66) 
i Jal cos (A «i 


wird. Wir fihren fiir die Summe die Abkiirzung J, ein und fir die Logarithmen 


cos ie er, ~) 


Dann wird, wenn wir in der Summe noch die Glieder mit yi zusammenfassen, die in (29) bendtigte 


GroéBe 


In 


saline iis (67) 


6 
aa! . vane 5 pa) 
IpH = > prin 20. In (i, t) — sin 2 041 In (i,t + 1) + = ; (cos 2 0; — 008 2.0441) i (68) 


i=0 


wobei formal die bisher noch nicht definierten GréRen 


& =a, = 0 (69) 


| ; (p cos 2 x; — sin p) + Const. 
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gesetzt werden miissen. Ganz analog folgt mit (64) fiir das in (30) bendtigte Integral 


See ee + 5} (Pi — Qi—-1) sin 2 Xi 


eee 
ed cos (PEt, | 


20 4 

ae 28 cos ee eee 
eae) Sa Wa Ky Wig +: > — cos 2 4; In ) 2 ee 
; 


. 1 | Pi | ie 
ae = (cos y; — cos y;_1) — cos gi \In eee (3 — “| —In 7 + cos yy in cos - — x —Inv; 


(70) 
Auch hier haben wir die gleichen Vereinfachungen wie oben, nur verschwindet jetzt der letzte Term 


mit 1 = | sowie der vorletzte mit i = 6 nicht. Bezeichnen wir die Summe von (70) mit J, und fiihren 


die Abkiirzung (67) ein, dann ist 


6 
ae cos 2 0; [In (i, ?) —In  —1,)] +4 (i — Gi) sin2 


i=l 
+ In (0, 1) — In v, — In (6, 6) + In vg. (71) 
Daraus eliminieren wir v, und v, durch fortlaufende Anwendung von (57), was 
6 
% : : 1 
4 2 = (1 + cos 2 «;)[In (i, i) — In (i — 1, i)] 4 5 (pi — Yi_1) sin 2 x (72) 
i 


ergibt. Die zu (68) analoge Umschreibung fiihrt auf 


J,= »| (1 + cos 2 «;) In (i, i) + (1 + cos 2 «;41) In(i,i + 1) eo, (sin 2 x; — sin 20441)! Go) 


i=0 | J 
was wieder unter Beriicksichtigung von (69) in Ordnung ist. Damit lauten die SchlieBungsbedin- 
gungen (29) und (30) 
K, Wa =~ K, Weg = la le, (74) 
ie W =e KK, Ws a uh == (75) 
Aus (56) eliminieren wir ebenfalls durch laufende Anwendung von (57) v, und vg. Dies ergibt, wenn 
wir wieder unter Beriicksichtigung von (69) 
6 


J,= > {—In (i, 1) + In Gi, i + 1)} (76) 
i=0 
einfiihren, 
— K,nw,(0) + K, nw, (22) + J, =0. (77) 


Die Gleichungen (74), (75) und (77) enthalten als einzige Unbekannte K,, K, und @3. Sie lassen 


sich daraus bestimmen. Wir eliminieren zuerst K, und K, und erhalten 


(J-—m) D,—J, Dz + J, D3 =9 (78) 
mit 
D, = W,,\nw,(2 2) + Wee In w,(0) , (79) 
D, = W., In w,(2 2) + We, In w,(0) , (80) 
Ds = Wey Ws 7%, Wie Wy S (81) 
Gleichung (78) ist fiir y, eine transzendente Gleichung der Form 
a + b log cos (a5) + clog cos (Fe —ay)| + dos = 0, (82) 


deren Konstanten a, b, c, d aus (68), (73), (76) und (79) bis (81) berechnet werden kénnen. Wir ver- 
folgen dies in der nachsten Ziffer im Einzelnen. Aus (82) wird gs, mittels bekannter Methoden be- 
stimmt. Damit lassen sich die J,, J, und J, nach (68), (73) und (76) berechnen, mit denen 


K, = Ss (— In We(2 qt) Ss + W ap) ? (83) 
K, == x (— ln w,(0) ab aa W's1 J) (84) 


folgt. 
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Nun benétigen wir fiir die Bestimmung der v; zu den fiinf Gleichungen (57) hin noch eine weitere. 
Sie folgt aus der bisher nicht beniitzten Gleichung (28) zu 


Q; 


[pov = 3 | m9 


Pi-1 


K; In w(p)| dp | + 22 (In2—Inv,) =0. (85) 


Hierin kénnen die Geschwindigkeitsverhaltnisse aus (57) entnommen werden, gy, und die K; sind 
schon bekannt. Es ist also ein Wert v, definiert, aus dem mit (57) die iibrigen v; berechnet werden 
kénnen. 

Unser Ansatz gibt also insgesamt zur Wahl der w,(y), x; und ¢; auBer pz eindeutig eine zulassige 
Funktion P(g). Dazu ist mit (32) und (33) ein Profil definiert. Wir beschreiben anschlieBend noch 
den numerischen Rechnungsgang zur Gewinnung dieses Profils und zeigen in Ziffer 5, was der An- 
satz fiir die praktischen Aufgabenstellungen leistet. 


4. Die praktische numerische Rechnung. Die praktische Rechnung unseres Verfahrens besteht 
aus drei wesentlich verschiedenen Teilen. Es mu erst die Funktion P(g) nach (54) mit allen Kon- 
stanten berechnet werden, sodann aus ihr Q(p) mittels (34), und damit ist die Integration (32) 
und (33) durchzufihren. 

Wir teilen das Intervall (0, 2 7) in 2 N Teile und arbeiten an den Stellen 

=o = (y= 0 leer ee Nae (86) 
Die grade Zahl N darf nicht zu klein sein, wenn die Genauigkeit des Verfahrens gut sein soll. Es ist 
dann keine Einschrankung, wenn wir die Stellen gy; (auBer p,) von (50) an gewissen Stellen t, wahlen. 
AuBer den q; seien uns die x; von (51) und die w;(~) von (52) gegeben. 

Im Verlauf unserer Rechnung benétigen wir an den Stellen t, die folgenden Funktionen, die wir 
gleich zu Beginn aufschlagen: 


Ge 


b) cos G _ x fiir alle x; und alle t, fiir die Gleichungen (32), (33) und (47), 


a) log jeweils nur zwischen ~;_; und q; fiir die Bildung von P((), 


c) sin ™ fiir (32) und (33), 


d) log w,(y) fiir P(p) und die Integrale (61) und (62). 


Die dekadischen Logarithmen reichen aus. Wir miissen nur ganz vereinzelt durch Multiplikation 
mit In 10 auf die natiirlichen iibergehen. Die Stelle gy ist noch nicht bekannt. Man begrenzt daher 
die Bereiche mit x; und «, vorerst durch einen benachbarten Wert T,, den man wegen (60) festlegen 
kann. An den anderen Intervallgrenzen entstehen die fiir die In (i, 7) nach (67) nétigen Logarithmen, 
die wir entsprechend (67) mit 


log lene (3 a) SORE. F) (87) 


bezeichnen. 
Die Integrale W,; und W,; nach (61) und (62) kénnen je nach Wahl der Funktionen wip) ge- 


schlossen ausgewertet werden oder mittels der Interpolationsmethoden aus den Werten 
cos T, log w,(t,), sin T, log 1;(T,) 


gewonnen werden', Die Werte cost, und sin, sind in denjenigen von sin T,/2 enthalten. Beim Er- 
gebnis mu®B noch der Faktor — In 10 beriicksichtigt werden. 

Mit den W,;, W.; und den anderen, bis auf den Faktor In 10 schon aufgeschlagenen GréBen 
der Gleichungen (79) bis (81) bilden wir die dortigen D;. Die transzendente Gleichung (82) wird 


zweckmabig in der folgenden Reihenfolge aufgestellt. Man berechnet erst die GréGen der Gleichung | 
J.—% = a, + b, log (3, 3) + ¢, log (3, 4) + d. gg (88) | 


wo a, die Summe (68) ohne das Glied mit i = 3 ist, b,c, und d, die Anteile des Gliedes i — 3 sind. | 


Analog sind nach (73) fiir J, die GréBen a,, b,, c, und d, zu berechnen, nach (76) die a,, b, und ¢,; d, 


verschwindet. Alle diese GréBen kénnen mit den schon vorhandenen log (i,j) und den wenigen 


1 C, Runge und H. Kénig, Vorlesungen iiber numerisches Rechnen, Berlin 1924. 
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noch aufzuschlagenden Werten sin 2 ,, cos 2 x; und @; (im Bogenmab 


ae os leicht berechnet werden. 
Mit ihnen werden die in (82) benétigten Grofen a, b, esd 


a— DD, 4.— D,a, + D 

und analog die b, ¢, d. ey oe 

:. Die transzendente Gleichung (82) fiir ~; wird am einfachsten mittels der Regula falsi aufgelést. 

Anhaltspunkt fiir die Ausgangsnaherung fiir y, gibt (60). Der Arbeitsaufwand ist dabei 
ring. 

Mit der Lésung von (82), die automatisch auch die log (3,3), log (3,4) und @, gibt, berechnet man 
entsprechend (88) die J,—a, J; und J, und mit ihnen die K, und K, nach (83) und (84). Zur Kon- 
trolle geht man mit den gewonnenen Zahlen noch in (74) ein. 

Mit der Lésung des Gleichungssystems (74), (75), (77) berechnet man miihelos die Werte 


z, ) 
cos | > Oi 


fiir welche die wesentlichen Zahlen schon aufgeschlagen sind. Inshesondere kénnen wegen (57) die 


Pr P(r, )— log 


log - K; log w;(t,) , (90) 


log “+! = log (i, i + 1) —log (i, i) (91) 


vi 
aus den log (i, 7) berechnet werden. Der Zusammenhang der Zahlen aus (90) mit P(g) ist 
P, = P(z,) = In 10 (P(z,) + log 2 — log »,). (92) 


Er mu8 aber gar nicht explizit berechnet werden, da beim Ubergang (34) auf Q(~) eine additive 
Konstante keine Rolle spielt. Die multiplikative Konstante nehmen wir bei der nun folgenden Be- 
rechnung von Q(¢~) mit in (34). 

Gleichung (34) werten wir in bekannter Weise naherungsweise aus!, indem wir die P(r,) durch 
ein trigonometrisches Polynom mit 2 N Koeffizienten interpolieren, bei diesem auf die konjugierte 
Form entsprechend (24) und (25) tibergehen und dann Q(y) wieder an den Stellen p = 7, summieren. 
Dieser ProzeB 1aBt sich in der einfachen Form 


N-—-1 
1 LTE 
OF OG) == ND (Pia Pe) CES oa (u ungerade) (93) 


darstellen. Wir benétigen in (32) und (33) Q(y) im GradmaB, weil hierfiir bessere Tafeln vorhanden 
sind?. AuSerdem beniitzen wir an Stelle der P, die P, von (90), also folgt 


UNG es Ee 
Q, = 2 (Prin os Pa lene : (Paton a Pr) (94) 
a 


worin die Faktoren 
180° In 10 wn 
ctg = 
— IU N 2) 2 N ‘ 
0 (u gerade) 


(uw ungerade) 


ky (95) 
fiir ein einmal gewahltes N ein fiir alle Male berechnet werden kénnen. In Tabelle 1 sind sie fiir 
N = 24 und N = 30 angegeben. 

Die Berechnung der Q, nach (94) ist der Hauptaufwand des Verfahrens. Die Anzahl N darf 
nicht zu klein gewahlt werden, damit am Profil geniigend Punkte berechnet werden. Nimmt man 
etwa N = 24, was schon eine Mindestforderung ist, dann erfordert (94) 48 Produktsummen zu je 
12 Summen und 12 Produkten. Wenn méglich kénnen hicr gréfere Rechenautomaten vorteilhaft 
eingesetzt werden. Das IBM-Lochkartenverfahren hat sich beispielsweise gut bewahrt®. 

Als Probe fiir die Q, dient die Umkehrung von (94) 


= 180°Inl0\?/— 1 15 
(Or7 oem Oe) ky, ir Ao (», ~~ d ies ° (96) 
y= o= 


1]. E. Garrick, FuBnote 1 von S. 36. 

2 Am besten bewahrt hat sich diejenige von A. Danjon, Tables des fonctions trigonometriques, valeurs 
naturelles & 6 décimales de centiéme en centiéme du degré nonagésimal, Paris 1948. 

3 Diese Arbeiten konnten mit freundlicher Genehmigung des Vorstandes der Allianz-Lebensversicherungs- 
A.-G. in der dortigen Holerith-Abteilung mit einem elektronischen Rechenlocher IBM 604 ausgefiihrt werden. 
Der genannten Gesellschaft, insbesondere Herrn Dr. Kanngiefer, der die notwendige Schaltung hergestellt hat, 


bin ich auBerordentlich dankbar. 
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Die letzte Summe ist darin fiir alle geraden und ungeraden y jeweils gleich, denn die P, sind mit 
2 N periodisch. 

Uber die Genauigkeit von (94) macht Garrick! Angaben. Sie ist dann schlecht, wenn die Fourier- 
Reihe von P(g) nicht oder im praktischen Sinn schlecht konvergiert. Lighthill? geht direkt von der 
Funktion In v(p), die im Staupunkt eine Singularitat hat, auf die konjugiert harmonische tiber, was 
nur von Fall zu Fall mittels des uneigentlichen Integrals (34) und verschiedener speziell tabu- 
lierter Integralfunktionen méglich ist. Unser 
geandertes Verfahren erméglicht die Befriedi- 
gung der vorangehenden Stetigkeitsforderun- 
gen und damit die Anwendung von (94) mit 


Tabelle 1. Die k, fir N = 24 und N = 30. 
N = 24 (Nie) 30 


big kul] i kul] geniigender Genauigkeit. 
1 83,86828 1 83.9147 Fiir die Integranden von (32) und (33) 
3 27,63537 3 27,76544 fehlt nun nur noch die Funktion v*(y) nach 
5 16,19369 5 16,41212 (52). Wir schreiben sie in der Form 
7 11,14685 7 11,45617 
9 8,22687 9 8,03080 Pye GL 
ul 6.26814 in 6.77173 SCI Beis ee 
13 4,82075 13 5,43060 : ee : ; : 
15 3.67299 1s 4.39761 Die Geschwindigkeitsverhaltnisse sind wegen 
17 2,71083 17 3,56112 (91) in den bekannten log (i, 7) auszudriicken. 
19 BOE ms 2,85585 Die wki kénnen logarithmisch aus den schon 
‘ Ronee 53 eae vorhandenen Funktionen log w; berechnet 
: 95 1.17834 werden. Die GréBe v, geht in (32) und (33) 
27 0,69651 nur multiplikativ ein, reguliert also den Mab- 
29 0,23047 stab der z-Ebene. Derselbe ist nur wegen der 


GréBe 1/4 von (36) bedeutsam, die nicht sehr 
genau bendtigt wird. Wir werten (85) naherungsweise aus. Wegen (92) kénnen wir auch hier an- 


statt der P, die P, verwenden. Wir erhalten dann statt (85) 


a_i 
log v, = log 2 + Set P,, (98) 


v= 1 


wenn wir das Integral mit der einfachen Rechteckformel auswerten. Dies ist gleichbedeutend mit 
der Interpolationsintegration beliebig hoher Ordnung, da bei der Integration einer periodischen 
Funktion iiber ihre Periode die Differenzen héherer Ordnung keinen EinfluB haben®. 


Es sei bemerkt, daB (97) und (98) auch schon vor dem Ubergang (94) auf die Q, berechnet wer- 
den kénnen. Wenn man sich also tber die aus den gewahlten AnsatzgréBen entstehenden Vertei- 
lungen u(y, x) einmal nicht sicher ist, kann man dieselben vor dem hauptsachlichen numerischen 
Aufwand schon aus (97) mit (47) berechnen. Sieht man ihnen irgendwelche Nachteile an (vgl. auch 
Ziff. 5), dann kann man notfalls die Rechnung mit neuen Eingangsgréen neu beginnen. Andert 
man dabei nicht alle GréBen, dann kann man manchmal einen Teil der vorhandenen Rechnungen 
wieder verwenden. Andert man beispielsweise nur die w,(y), dann bleiben die a,, b,, ¢., de von (88) 
und die analogen GréBen von J, und J, erhalten. Halt man dagegen die w;(p) fest, dann bleiben 
auch die D; von (79) bis (81) gleich. 


Mit (97) und (98) kénnen wir die Integranden von (32) und (33) vollends punktweise berechnen. — 


Wir fiihren diese Rechnung aber noch nicht durch, sondern nehmen bei ihr noch auf einige prak- 
tische Belange Ricksicht. 


Als erstes lassen wir in (32) und (33) den Faktor 4 weg. Das Profil erhalt damit eine Lange 
I, hase (99) 


mit der aus (36) 


1 


folgt. 


1 FuBnote! von S. 36. 
2 FuBnote > von S. 32. 
* Ygl. auch spater bei der Integration von (32) und (33). 


wee | 
Ca = 7 2 Hsin (100) | 
| 
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Das sich somit ergebende Profil hat fiir « = 0 verschwindenden Auftrieb, wie wir schon frither 
festgestellt haben. Stark gewélbte Profile, die fiir positive « oder c, berechnet sind, bekommen 
deshalb eine stark nach unten geneigte Vorderseite. Dies kénnen wir bei der Integration schon ver- 
meiden, indem wir in (32) und (33) den Steigungswinkel iiberall um einen Winkel 6 verkleinern. 
Wie gro8 6 ungefahr sein muB, damit Nase und Hinterkante gleich hoch zu liegen kommen, kann 
aus den berechneten Steigungswinkeln entnommen werden. Das Profil hat in der neuen Lage bei 
einem Anstellwinkel ~« = —f verschwindenden Auftrieb. 

Wir verwenden demnach anstatt (32) und (33) fiir die Integration an den Stellen 7, die GréBen 


fone «(z,) —-=— sin 2 es (F », 5) cos i + Q,— p) > (101) 
F . eatonl U 1 : 
f= y(t.) = —~ sin, ke ( 5 a) oe) sin (3 Oye p) ? (102) 


die wir mit unseren seitherigen Ergebnissen leicht berechnen kénnen. 

Aber auch die Integration fithren wir noch nicht durch, sondern verfolgen see zunichst ohne 
Rechnung und sehen dabei, welcher Rechnungsgang fiir das endgiiltige Ergebnis am einfachsten 
wird. Wir bilden zu der bekannten Integration durch Interpolation! nach Tabelle 2 mit den x, von 
(101) die ersten und zweiten Differenzen. Damit wir die letztere auch fiir y == 0 bekommen, ver- 
langern wir das Schema riickwarts mit » = 2 N— 1. Nun ist, wenn wir x(0) = 0 setzen, 


x, = x(q) wh Ogee) Se (Ate 4%) 
trainee : () 
Xp = X(T) wh E (x, + %) 94 (A?x, + Aris) + Xs 
und so weiter, wo h die Differenz des Arguments 
pe (104) 


N 


ist. Diese Integration ist leicht fortlaufend zu bilden, wobei man erst mit den Funktionswerten ganz 
durch rechnet, dann mit den zweiten Differenzen. Die jeweils konstanten Faktoren kénnen dann 
in der Maschine bleiben. 


Die entstehende Funktion x(y~) muB periodisch sein. Tabelle 2. Differenzenschema fiir die 
Es muB, wenn sich das Profil schlieBen soll, x) = xy Integration (103). 

sein. Dies wird nicht genau erfillt sein, da wir bei der - x As Wx 

Berechnung von Q(¢p) und bei der Integration Naherungs- 
verfahren angewandt haben. Die entstehende Differenz 2 ee 1} %.n-4 A ke os XyN—4 

Xo ° Xo 

&, = Xon — Xo (105) 1 xy vee Ae 

ae! 2 ky %2 ke 

ist einerseits ein Ma fiir die Genauigkeit der Rech- 3 E phe we . 

nung, andererseits eine véllig unabhangige Rechnungs- 4 ity Ag A? x, 

kontrolle. Ist sie klein genug, so daB wir nicht auf einen 5 ve 5 Arne 


Rechenfehler hingewiesen werden, dann gleichen wir sie 
noch auf einfache Art und Weise aus. Wie wir an den Differenzen leicht feststellen, tragen die 


héheren nichts zu ¢, bei, es muB vielmehr 
2N 
S, = + %, = 0 (106) 


y=1 


sein. ; 
Damit legen wir unseren Rechnungsgang folgendermaBen an. Wir bilden zuerst mit den GréBen 


aus (101) nach (106) S,. Der daraus folgende Fehler ist nach (103) 
é, = W Se. (107) 


Wir beurteilen an ihm, ob er der Rechenungenauigkeit zugeschricben werden kann. In Anbetracht 
der Grenzschichtdicke an der Hinterkante kénnen wir als sicher unbedeutenden Fehler 


&, = 0,003 1 (108) 


1 C, Runge und H. Kénig, FuBnote 1 vonS, 42. 
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i 
i} 
zulassen, also bei einem Profil von 1 m Lange eine waagrechte Verschiebung von 3 mm. Ist &€ leit 
genug, dann gleichen wir die %, aus, indem wir von allen x, nach (101) auBer dem verschwindenden 
Pedic GréBe 
at ere 109 
Osta eor ee 
abziehen. Dann ist (106) sicher erfiillt und die entstehende Funktion x(p) wird periodisch. 
Nun erst bilden wir die Differenzen und berechnen die x, nach (103), wobei wir x) = 0 setzen. 
Dann geht x(p) von Null bis zu einem Minimum 


min <= aa a (110) 
und wieder zuriick nach Null. Die GréBe x,in wird im allgemeinen nicht an einer Stelle x, ange- 
nommen. Wir interpolieren deshalb «(p) in der Nahe des Minimums durch eine Parabel und nehmen 
deren Minimum. Sei x, das kleinste der entstehenden x,, dann folgt 


(X41 ae Kia)" (111) 


ie ‘—1 =x, } 8 (2 Xn — Xyi4 — Xy1) 


Damit fiihren wir leicht vollends die Transformation 
x 100 (7 n 1 (112) 
a 


durch. Die Funktion x(y) geht dann wie iblich von 100 bis 0. 

Die zugehérige Funktion Y(p) entsteht ganz analog. Wir berechnen erst S, analog zu (106) | 
und daraus ¢, und 6,. Beziiglich ¢, miissen wir, was praktisch auch meist zutrifft, etwas héhere 
Genauigkeit verlangen als beziiglich ¢,, etwa 


ey < 0,002 1. (113) | 


Bei der Integration nach (103) nehmen wir an Stelle von h den Faktor 100 2 , so daB sofort die Y, | 
1 } 

in der den X, entsprechenden GréBe entstehen. 
Nach (47) kénnen wir vollends die Verteilungen v(g, x) und v(¢, x) fiir alle « berechnen. Fiir die 
in x«*(~) beniitzten a; sind die UmrechnungsgréBen schon aufgeschlagen. Das Nullmoment ist nach | 
(42) und (43) zu berechnen. Anstatt (43) nehmen wir die Naherung 


eevee init S32 Dede Pah (114) | 


vA 


welche in der Genauigkeit der Gleichung (93) entspricht. 

Das gesamte Verfahren bringt also auch bei der numerischen Auswertung keinerlei Schwierig- | 
keiten. Ein eingearbeiteter Rechner kann es mit allen Einzelheiten ohne Verwendung grofer 
Rechenautomaten in etwa 30 Arbeitsstunden erledigen, wenn er mit N = 24 arbeitet. 

Die Wahl der GréBe N hangt von der ESI: Genauigkeit ab. Man wird sie aber nicht ohne 
weiteres gréBer als N = 30 wahlen, da sonst der Rechenaufwand zu sehr steigt. Man kann viel- | 
mehr die Genauigkeit des Ergebnisses noch auf eine andere Art und Weise verbessern. Die gréfRte_ 
Ungenauigkeit unseres Verfahrens entsteht durch einen starken Knick der Funktion P(g) bei. 
~ =z der durch ein trigonometrisches Polynom nicht gut interpoliert wird (vgl. auch Ziff. 5, | 
Abb. 5). Da diese Interpolation in (93) enthalten ist, entstehen in den Q, kleine Schwankungen. | 
Sie werden bei der Intergation vollstandig geglattet und stéren am Ergebnis nicht mehr, sie be- | 
dingen aber hauptsachlich die Fehler ¢, und ¢,. Will man sie verkleinern, dann spaltet man von | 
P(p) vor dem Ubergang auf Q(¢) einen geschlossenen Anteil ab, der den starken Knick enthalt und 
fiir den (34) geschlossen auszuwerten ist. Dies kann in der Form 


PO) (115) 


mit 


1 'p  [ | | 
2 tal0 te(2 %s] = X4), (116) 


geschehen, was auf | 


| 
Oy) = — a In 10 sin es In ak QM(y) (117) | 


| 
| 


<XV. Band 1957 Eppler: Direkte Berechnung von Tragfliigelprofilen aus der Druckverteilung 47 


fiihrt. Dabei hat P®(p) bei p = ep, stetige Tangente und QM(~) kann daraus mit den abgeleiteten 
Formeln (94) und (95) berechnet werden. Diese Abspaltung ergibt im allgemeinen mit weniger 
Aufwand giinstigere Verbesserungen als eine Erhéhung von N. Bisher hat sie sich aber bei keiner 
Beispielrechnung als nétig erwiesen. 


Unser Ansatz fihrt also insgesamt mit ertraglichem Aufwand zum Ziel. Wir zeigen anschlieBend, 
was er bei der Profilberechnung leistet. 


3. Die erfafien Méglichkeiten des Profilentwurfs. Je nachdem die GréBen a, ~; und wg) 
der vorangehenden Ziffern gewahlt werden, entstehen Profile mit verschiedenen Kigenschaften. 
Wir beschreiben in dieser Ziffer die wichtigsten in unserem Ansatz enthaltenen Méglichkeiten des 
Profilentwurfs qualitativ. In Ziff. 6 geben wir ein numerisches Beispiel. 

Wir haben durch (50) den Kreis der ¢-Ebene in sechs Teile eingeteilt. Wegen (8), was auf die 


Normalkoordinate X transformiert 


X wx 50 (1 + cos ¢) (118) 


ergibt, verfiigen wir mit den Grenzen ~; von (50) ungefahr iiber die Grenzen X; = X(g;) der ent- 
sprechenden Bereiche am Profil. Nur gy; haben wir im Verlauf unseres Verfahrens festgelegt. Die 
Stelle X; liegt ungefahr bei der Profilnase X,; ~ 0. Wir haben also auf Ober- und Unterseite je 
drei Bereiche, deren Grenzen wir willkiirlich festlegen kénnen. 

In den einzelnen Bereichen wahlen wir mit (51) wegen (47) diejenigen Anstellwinkel, bei denen 
das v*(~) des Ansatzes (52) dort als Geschwindigkeit auftreten soll. Abgesehen vom ersten und 
letzten Bereich in (50) wird v*(p) = v; wegen (53) konstant. In den beiden Nachbarbereichen der 
Nase und den daran anschlieBenden wahlen wir also lediglich gewisse Anstellwinkel «;, bei denen 
dort konstante Geschwindigkeit herrschen soll. Die Gré8e v; derselben wird durch das Verfahren 
segeben. 

Benétigen wir auf Ober- und Unterseite nur jeweils einen solchen Bereich mit konstanter Ge- 
schwindigkeit, dann wahlen wir 


Xy = O3, X= Os > (119) 
was mit (57) 
v= ee, CP (120) 


ergibt. Die beiden jeweiligen Bereiche haben dann beim gleichen Anstellwinkel die gleiche konstante 
Geschwindigkeit, flieBen also in einen einzigen zusammen. Die Stellen ~, r0,(@) 
und ~, haben dann keinerlei Einflu8 mehr und miissen nicht mehr fest- I 
gelegt werden. 7 

Uber die genaue Form der w,(y) wird noch zu sprechen sein. Da sie ! 
bei normal dicken Profilen an der Hinterseite immer einen gewissen | 
Geschwindigkeitsabfall erzeugen miissen, denken wir uns w,(~) zunachst | 


. r 5 : Pp 
in der Form der Abb. 2, wo auBer (55) giinstig noch 5 yes 
7 
dw, = (12 1) Abb. 2. Ungefahrer Verlauf einer 
‘do a = brauchbaren Funktion w,(g). 
\ /p=9 


ist, damit beim Ubergang auf X keine senkrechte Tangente in v(X) auftritt. Ganz analog kénnen 
wir auch w,(y) einfiihren, 
Fiir die Anstellwinkel x, und a, kénnen wir fast immer 


a (122) 


voraussetzen, was nach (57) 
Un Ons Use (123) 


gur Folge hat. Wir legen dadurch nur die durch w,(p) und w,(P) erfaBte Form des Geschwindigkeits- 
abfalls an der Hinterseite bei denjenigen Anstellwinkeln fest, bei denen in den Bereichen davor 
konstante Geschwindigkeit herrscht. Der Geschwindigkeitsabfall wird dabei nach (52) proportional 
gu den v;, was bei der Wahl von w((p) zu beriicksichtigen ist. ook ¥ 
Symmetrische Profile erhalten wir durch eine weitere Spezialisierung. Bei ihnen ist bei « = 0 
die Geschwindigkeit auf Ober- und Unterseite gleich, auBerdem mu die Abbildung 2(C) symmetrisch 
sein und X(y) = X(— ¢) gelten. Damit wird v(p, 0) und nach (27) auch P(g) eine gerade Funktion. 


In (24) verschwinden alle 6,,. In unserem Ansatz erfordert dies insgesamt wegen der Gleichungen 
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(50), (51), (52), (61) und (62) 


oy =—%, Gr =2A—-Y; w,(~) = we(27 — @) » | 
on mmm Xo Po =2N—G » W., =— Wg, (124) _ 
a see A 8 ee Wey = Wee - 


Im einzelnen wollen wir damit die folgenden Méglichkeiten des Profilentwurfs unterscheiden. 
a) Kinfache,symmetrische Laminarprofile. Wir wahlen unter Beriicksichtigung von (124) | 
iy =n = Ae (125) 
eine GréBe ,, die nach (118) den Beginn X, des Geschwindigkeitsabfalls festgelegt, und ein w,(~) 
etwa entsprechend Abb. 2 mit (121). Der Rechnungsgang der Ziff. 4 zur Bestimmung von P(¢) 
und v*(~) lauft damit automatisch. Es wird wegen (57) 
Vv; =U, (126) 


fir alle i. Die Gleichungen (75) und (77) sind mit p,; = a und K, = K, erfiillt. Aus (74) erhalten | 
wir ein K, und aus (98) ein v,. Nach (52) entsteht die in Abb. 3 skizzierte Verteilung v*(p) (dicke | 


Oberse/te 


Unterse/te 


+i, TWre0C, 


Abb. 3. Die Funktionen v*(g), v(~, «,) und v(~, «,) bei einem einfachen, Abb. 4. Die Funktion v(X, «,) bei einem einfachen, 
symmetrischen Laminarprofil (schematisch). symmetrischen Laminarprofil (schematisch). 


—= 


Linien), die symmetrisch zu g = ist, und die mittels (47) die gleichfalls eingezeichneten Vertei-+ 
lungen v (gy, X,) (durchgezogene Linien) und v (¢, xg) (gestrichelte Linien) liefert. Gehen wir vollend& 
auf X iiber, dann entsteht der in Abb. 4 skizzierte Verlauf v(X, «,), der genau der Forderung de 
bekannten Laminarprofile entspricht!. Beim Anstellwinkel a, ist die Geschwindigkeit auf deri 
Oberseite von der Nase bis X, konstant. Fiir x, = —«, entsteht das gleiche Bild, nur sind die Linieni: 
fiir Ober- und Unterseite vertauscht. Bei allen Anstellwinkeln « zwischen a, und —o, steigt diey 
Geschwindigkeit v(X, «), wie man mit (47) leicht nachpriift, auf Ober- und Unterseite an bis zur) 
Stelle X,. 
Fiir die endgiiltige Form des Geschwindigkeitsabfalls an der Hinterseite (X, < X < 100) is 
auBer der gewahlten Funktion w,(y) auch die sich ergebende GréBe K, mafgeblich, die ja als: 
Exponent auftritt. Vom grenzschichttheoretischen Standpunkt aus? darf die durch w,(p) ge+ 


gebene Form nicht zu sehr gestért werden. Man ist also daran interessiert, daB K, ~ 1 wird. 


Dies kann man erreichen, wenn man w,(~) von vorn herein etwa so groB wahlt, wie es der erfah-. 
rungsgemaB notwendige Geschwindigkeitsabfall mit K, = 1 verlangt, was nach einigen Beispieler 
leicht gelingt. Entsteht je einmal ein K,, das w{(p) zu sehr verformt, dann bemerkt man dies jat 
schon nach dem ersten Rechnungsgang, der Auflésung von Gleichung (82), also vor der hauptsich-: 
lichen Rechenarbeit. Man kann dann leicht mit einem geanderten w,(y) neu anfangen. 
Die insgesamt mit der Wahl (125) berechenbaren symmetrischen Profile entsprechen ungefahui 
denjenigen des NACA-Report 8243. Hier kénnen wir nur mit w,(y) die Form des Geschwindigkeits+ 
abfalls noch variieren, und auBerdem geht die konstante Geschwindigkeit bei « = x, genau bis 
zur Nase m =, wahrend dies bei den NACA-Profilen nicht ganz erreicht werden konnte. Diesen) 
letzte Unterschied erscheint vielleicht unwesentlich. Gehen wir jedoch auf die Funktion Pp 
nach (27) iiber, dann sieht der Einflu8 wesentlich gréBer aus, weil die nachste Umgebung der Nase 


1 NACA-Report 824, FuGnote 1 von S. 38. | 
2 F. X. Wortmann, Z. Flugwiss. 3 (1955) S. 333. | 
3 FuBnote 1 von S. 38. 
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schon einen beachtlichen y-Bereich erfiillt. Bei der in unserem Verfahren beniitzten Funktion P(g) 
nach (54) entsteht bei y = z infolge der bis dorthin konstanten Geschwindigkeit v(y,x,) ein spitzen- 


artiger Knick (Abb. 5), denn das Argument von cos ra, geht nahe an 7/2 heran. Diese Spitze 


hiefert einen erheblichen Beitrag zum Integral K;. Sie geht verloren, wenn v (@, o,) gegen die Nase 
abfallt (gestrichelt in Abb. 5). Dann miissen entsprechend die K und damit der Geschwindigkeits- 
abfall an der Hinterseite gréBer sein, damit (74) erfiillt ist. Man kann also durch sorgfaltige ean 
tung der Nase einigen Druckanstieg auf der Hinterseite sparen. 

Die Spitze von P(y) wird zwar bei unseren numerischen Rechnungen durch die Interpolation 
auch nicht exakt ausgeschépft, aber wir verwenden in ihrem Bereich immer mehrere Punkte, so daB 


wir jedenfalls dem Idealfall sehr 

nahe kommen. Gegebenenfalls AFP) 
k6énnen wir immer mit der Ab- 4 
spaltung (115) noch genauer 
arbeiten. 

Der geringere Druckanstieg 
an der Hinterseite entspricht 
physikalisch einem diinneren 
Profil. Wir gewinnen also hier 
bei vorgegebener Aufgabenstel- 
lung das diinnstmégliche Profil. 

Beiden NACA-Beispielenist - 
die Form des Profils, insbeson- 
dere seine Dicke durch die Auf- 
gabenstellung und das_ Be- 
rechnungsverfahren festgelegt. 
Grobe «a,-Werte, also grobe 
Ca-Bereiche mit Druckabfall- ' ; , z 
Meeker kleine Cpanel fraheron Vetahion Gaucher), 
groBe X, und lange Druck- 
abfallstrecken ergeben dicke Profile und umgekehrt. Auch bei unserem Verfahren ist das Profil 
mit «,, ~, und w,(~) eindeutig festgelegt und der EinfluB von a, und ¢, bleibt der gleiche. Die 
Berechnung des diinnstméglichen Profils zu einer Aufgabenstellung bedeutet also gleichzeitig, dal 
unter Profilen mit gleicher Dicke und gleicher Lange der Druckabfallstrecken das hier berechnete 
den gréBten Wert «,, also den gréBten c,-Bereich mit Druckabfallstrecken hat. 

Die Berechnung von einfachen symmetrischen Laminarprofilen laBt sich also mit dem vorlie- 
genden Verfahren vollstandig ausschépfen. 


b) Einfache, unsymmetrische Laminarprofile. Unsere Formeln lassen ohne weiteres 
die Ausdehnung der unter a) getroffenen Voraussetzungen auf den unsymmetrischen Fall zu, der 
damit genau so erschépfend zu behandeln ist wie der symmetrische. Wir bendtigen wieder auf 
Ober- und Unterseite jeweils nur einen Bereich mit konstanter Geschwindigkeit, wahlen also wegen 
(122), (119) und (59) 

OX = Xy = Wy > Ky = Os = Oe - (127) 


Anstatt (124) sei nun die dort verschwindende GréBe 
IL 
= (x +05) >0. (128) 


Dazu wablen wir p, und g;, wodurch wir wegen (118) iiber den Beginn des Druckanstiegs auf Ober- 
und Unterseite verfiigen. Die Funktionen w,(p) und w,(p) kénnen ahnlich wie bei a) (Abb. 2) voraus- 
gesetzt werden. 


Nach (57) wird mit (127) 
Up Ue gs 


Die GréBen gs, K,, Kg, v2 und v, ergeben sich aus unserem Rechenverfahren. We calc ea 

ahnliche Geschwindigkeitsverteilungen erhalten wie bei a), nur ist jetzt die Symmetrie von v (~) 

zu y = 7 gestort. Wegen (128) ist der Betrag des Auftriebs bei « = a, gréBer als bei « = a5, es 
A 
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mu& deshalb 
Ue Uz (130) 


| 


sein. Entsprechend wird der Geschwindigkeitsabfall auf der Oberseite gréBer sein als auf der Unter-— 


seite, denn die Geschwindigkeit v(g, x) ist an der Hinterkante wegen (56) stetig. Ks entsteht im 
Prinzip der in Abb. 6 skizzierte Verlauf von v*(p) (dicke Linien), und mittels (47) lassen sich auch die 
Verteilungen v(p, x9) (durchgezogen) und v((, x5) (gestrichelt) leicht tibersehen. Gehen wit auf xX 
iiber, dann entstehen ungefahr die Verteilungen v(X,«.) und v(X,«;) von Abb. 7, die wieder 
genau der in a) fiir symmetrische Laminarprofile beschriebenen Forderung entsprechen. Bei allen 
Anstellwinkeln zwischen x, und x, steigt die Geschwindigkeit von der Nase ab auf der Oberseite 
bis X,, auf der Unterseite bis X, an. 


Up) = lp) 


is | 
Val Tos Ps x i TG 
Abb. 6. Die Funktionen v*(p) und v(p, ;) bei einem einfachen, unsymmetrischen Abb. 7. Die Funktionen v(X, 7) bei einem einfachen, 
Laminarprofil (schematisch). unsymmetrischen Laminarprofil (schematisch), 


Auch iiber die Form der entstehenden Profile kann man von vornherein einiges aussagen. Kin 


groBer Wert (128) fiihrt auf hohe Auftriebsbeiwerte, also stark gewélbte Profile. Der GréBe ag, 
die in a) wichtig fiir die Profildicke war, entspricht nun die halbe Differenz der Anstellwinkel 
1/2 (x,—«a5), die mit (124) ina, iibergeht. Ist sie groB, und ist gp, und 27—q; klein, dann entstehen 
wieder dicke Profile und umgekehrt. Alle berechneten Profile sind im gleichen Sinn optimal wie 
beia). 


Solche unsymmetrischen Profile wurden bisher nur naherungsweise berechnet. Nach der voran- 


gegangenen Arbeit! waren sie zwar mit konformer Abbildung zu berechnen, die Umgebung dex | 
Profilnase bringt dabei aber, vor allem bei dickeren Profilen, ein unbeabsichtigtes Abfallen von v(@,0a), | 


was wieder der nicht ausgeniitzten Spitze von Abb. 5 entspricht. Mit dem vorliegenden Verfahren 
ist auch diese Aufgabe erschépfend zu lésen. 


c) Zweifache Laminarprofile. Die nach a) und b) zu berechnenden Profile haben einen 


gewissen Nachteil. Ihre Dicke hangt eindeutig von der Aufgabenstellung ab. Profile mit groBer | 


Dicke haben immer einen groBen Abstand «1, — a; zur Folge, also einen groBen c,-Bereich mit lami- | 


narer Grenzschicht. Nun ist aber ein solch groBer Bereich haufig gar nicht erforderlich, beispiels- | 


weise bei Motorflugzeugen, die im Reiseflug fast nur bei einem einzigen c, fliegen. Verwendet man 
hierfiir die einfachen Laminarprofile von b) dann kommt als eine Méglichkeit ein relativ diinnes 


Profil mit kleinem «1, — a, in Frage. Abgesehen von den baulichen Nachteilen solcher Profile treten | 
bei ihnen fiir « >a, an der Nase rasch groBe Geschwindigkeiten auf, die keine hohen Auftriebsbei- | 
werte ohne Stroémungsablésung erlauben. Die andere Méglichkeit ware ein dickes Profil mit groBem | 


X%2—;. Solche Profile bedingen aber naturgemaB an den Flanken schon hohe Geschwindigkeiten 


und haben dadurch Nachteile. AuSerdem sind auch hier die «.-Werte noch nicht so groB, daB | 


giinstige Maximalauftriebe erreicht werden. 
Um diese Nachteile zu vermeiden, haben wir in unserem Ansatz auf Ober- und Unterseite je | 
zwei Bereiche konstanter Geschwindigkeit eingefiihrt. Wir erlautern die hierdurch entstehenden 


Méglichkeiten des Profilentwurfs zunachst an einem symmetrischen Beispiel. Unter Beriicksich- | 


tigung von (124) wahlen wir 


ie = 0 = = 06 (131) | 
1 R. Eppler, FuBnote 1 von S. 32. | 


q 
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as a 1, das wie bei a) den Beginn des Druckanstiegs festlegt, und ein q,, dessen nach (118) 
ae netes X» nahe bei der Nase liegt, etwa zwischen X, =5 und X, = 10. Nach (118) mue 

ann QP, ungefahr zwischen 140° und 155° liegen, also von x noch einigen Abstand haben. Die 
Funktion w,(p) wahlen wir wieder wie bei a) er 


Wegen (57) wird nach (131) 


. 


Vg > Vg = 0,, (132) 
und unser Verfahren gibt sonst a i ie bei ' : 
(py (dick), ca a gil en i dem gleichen Weg wie bei a) etwa das in Abb. 8 skizzierte Bild 
A? \ ey me mittels (47) die Verteilungen v(~,x;) leicht ermitteln lassen. Der Uber- 
gang auf X ist in - 9 dargestellt. Die Linien fiir x, und 5 sowie fiir x, und «, tiberdecken sich 
wieder mit vertauschten Ober- und Unterseiten. ‘ 

Ein solches Profil hat n reierlei Streck stei mdiokee i 
Pepence co eee etrecken ansteigender Geschwindigkeit. Bei allen Anstell- 
ae 2 te “ = — 4, Ssteigt die Geschwindigkeit auf Ober- und Unterseite bis 

i allen . mit « ; 5 } 3 5 
e a m7 winkeln & mut Pa > Os = — M2 Steigt sie zusidtzlich auf beiden Seiten 
1; Vie zweiten Strecken sind die langeren. Sie werden die hauptsachlichen Laminareffekte 


<i 
ve fl SS ocean 
Oberse/te aN 


\ 
CASS ee an \ 
a 


x 
Y1 Pe Pu Ps ie 700 
Abb. 8. Die Funktionen v*(g) und v(9, a;) bei cinem zweifachen, Abb. 9. Die Funktionen 0(X, o ;) bei einem 
symmetrischen Laminarprofil zweifachen, symmetrischen Laminarprofil 
(schematisch). (schematisch). 


bringen. Man kann aber jetzt x, so klein machen, wie es die jeweilige Aufgabenstellung erfordert, 
ohne zu diinne Profile und deren Nachteile befiirchten zu miissen. Um sie zu vermeiden, steht 
jetzt noch x, und @, oder X, zur Verfiigung. Mit diesen GréBen kaun man die Profildicke genau so 
regeln wie bei a) mit x, und X,. GroBe a, und X, ergeben divke Profile. Im allgemeinen wird man 
dabei x, wesentlich gréBer nehmen als x, bei den einfachen Laminarprofilen, so daB man schon mit 
kleinem X, auf brauchbare Profildicke kommt. Der Hauptvorteil der vorderen Druckabfallstrecken 
liegt namlich nicht in dem Gewinn der kleinen Strecke laminarer Grenzschicht. Vielmehr hat die 
maximale Geschwindigkeit, die am Profil in dem grofen x-Bereich zwischen x, und —«, auftritt, 
den Wert vg, der relativ klein ist, da er auf einem ganzen Stiick des Profils angenommen wird (vel. 
das Beispiel der nachsten Ziffer). Bei allen bisherigen Profilen sind in vergleichbar groBen «-Berei- 
chen an der Profilnase langst Geschwindigkeitsspitzen mit viel héheren Werten vorhanden. Wir 
gewinnen also durch den vorderen Bereich ansteigender Geschwindigkeit hohe Auftriebe ohne hohe 
Geschwindigkeitsspitzen. Diese Méglichkeit besteht immer dann, wenn die eigentliche Aufgabe so 
kleine x,-Werte, also so kleine c,-Bereiche mit langen Druckabfallstrecken verlangt, daB dadurch 
bei den einfachen Laminarprofilen von a) zu diinne Profile entstiinden. Diese Falle kommen haufig 
vor. Gerade bei Motorflugzeugen, wo man wie erwahnt vornehmlich bei einem einzigen c, auf 
geringen Widerstand angewiesen ist, hat man bei Landung und Start durch hohe Maximalauftriebe 
groBe Vorteile. 

Dabei erfolgt schon dann ein geniigender Abbau der Nasen-Saugspitze, wenn die Profildicke 
durch den zweiten Laminarbereich nur um wenige Prozente (2—4%) gréfer gemacht wird, als sie 
beim entsprechenden einfachen Laminarprofil (mit gleichem a.) ware. Mit dieser geringen Dicken- 
zunahme, die auch nur eine kleine Zunahme des Druckanstiegs auf der Hinterseite und damit des 
Widerstandes mit sich bringt, gewinnt man schon so hohe Auftriebsbeiwerte, wie sie bei ein- 

4* 
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fachen Laminarprofilen héchstens mit extrem hoher Dicke (21%) zu erreichen sind (vgl. das Bei- 
spiel von Ziff. 6). | 

Vergleicht man andererseits einfache und zweifache Laminarprofile mit gleicher Dicke und | 
gleichem X;,, also gleich langen Druckfall- und Druckanstiegsstrecken, dann hat das zweifache ein 
kleineres x, also einen kleineren Bereich mit geringem Widerstand, der Wert des geringsten Wider- | 
standes ist jedoch ungefahr der gleiche. Allein durch den Verzicht auf einen Teil des «-Bereiches 
mit kleinem Widerstand kann aber ein hohes c, gewonnen werden, es kénnen die Vorteile der 
diinnen und der dicken Profile miteinander kombiniert werden. 

Derartige ,,zweifache Laminarprofile“ kénnen mit unserem Verfahren ohne Schwierigkeiten 
auch unsymmetrisch berechnet werden. An Stelle des hohen qs tritt ein groBer Wert 1/2 (x3; — a4), 
an Stelle des kleinen , ein kleiner Wert 1/2 (x,.—«a;). Dadurch wird lediglich eine zusatzliche Wélbung 
der Profile erreicht, und die Geschwindigkeiten der Unterseite werden kleiner als diejenigen der 
Oberseite, genau so wie dies beim Ubergang von a) nach b) bei den einfachen Laminarprofilen 
geschah. Mit x; wird auch der Maximalauftrieb gegeniiber den unsymmetrischen Profilen noch 


10] 


0 Py L H Ps Ps ent A iy Xs 700 ' 
Abb. 10. Die Funktionen v*(g) BES v(p, X,) LO ono speziellen zweifachen Abb. 11. Die Funktionen v(X, o;) bei einem 
unsymmetrischen Laminarprofil speziellen zweifachen, unsymmetrischen ; 

(schematisch). Laminarprofil (schematisch). 


vergréBert. Mit x, und «, kann man denjenigen Anstellwinkelbereich begrenzen, in dem man mit 
langen Druckabfallstrecken und damit langen laminaren Grenzschichten geringen Widerstand 
erzeugen will, mit x; und «, erhalt man die Grenzen der tiberhaupt noch ohne gréBere Strémungs- 
ablésungen méglichen Anstellwinkel. 

Ist man dabei beispielsweise nur auf hohe Anstellwinkel angewiesen und an den niedrigen (Sturz- 
flug) nicht interessiert, dann kann man «, =a; wahlen und nur auf der Oberseite den zweiten | 


Laminarbereich einfiihren. Da auch durch die Wahl «, < «a; die GréBe = (x3; —o,) und damit die | | 


Profildicke vergréBert wirde, kann man bei Verzicht auf diese Forderung den Effekt auf der Ober- 
seite vermehren, also etwa x3 oder X, gréBer wahlen, um wieder zur gleichen Profildicke zu gelangen. 
Unser Verfahren bindet ja an keinerlei Symmetrie. 

Ein spezielles derartiges Profil entsteht, wenn man zudem noch a, = «;, also insgesamt 


wahlt. Mit «1, —a, = 0 hat man dann nur einen einzigen Anstellwinkel, bei dem auf Ober- und _ | 
Unterseite nicht abfallende Geschwindigkeit vorhanden ist, und man erzeugt nur mit groBem « | 
hohen Maximalauftrieb. Genau wie bei den vorhergehenden Beispielen evn man fiir tess Wahl 
die in Abb. 10 und 11 skizzierten Verteilungen v(g, «;) und v (X, «;) ibersehen. 

Man ist, wie in Ziff. 5a), auch bei den zweifachen Laminarprofilen daran interessiert, zu einer 
Wahl der Formen w,(y) und w,(p) des Geschwindigkeitsabfalls an der Hinterseite ede K,x1 
zu bekommen. Im symmetrischen Fall kann dies nun leicht mittels x; oder X rogubiere aenien | 
denn diese GréBen regulieren die Profildicke, die im wesentlichen fiir die GréBe iiés Drdckaneces aul | 
der Hinterseite verantwortlich ist. Im unsymmetrischen Fall kann es dagegen vorkommen, daB bei | 


i 
der Fes ron — (x Sage : d Shee 
estlegung von 5 (x; —«,), X_, und X, die eine der beiden GréKen K,, Ky wesentlich gréBer, | 
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die andere wesentlich kleiner als 1 bleibt. Wie dies ausgeglichen werden kann, wird in Ziff. 5 d) be- 
schrieben. 

Allgemein lassen sich mit dem zweifachen Laminarprofil eine Menge von Aufgabenstellungen 
lésen, und zwar durchweg optimal im vorher erwahnten Sinn. Die unsymmetrischen Profile, die 
frither groBe Schwierigkeiten bereiteten, sind in den Formeln automatisch enthalten. 

Man kann auch ganze Serien von Profilen berechnen, die ahnliche Kigenschaften, aber verschie- 
dene Dicke haben, indem man beispielsweise bei festgehaltenen aj, Y,, 5: w,(y) und w,(~) nur ~, und 
(oder) g, andert. Dabei spart man 
bei der wiederholten Berechnung von 4721/9) 7U5(P) 

P(y) einigen Rechenaufwand. Der 1 - i= 
eigentliche Ubergang auf das Profil 
mu jedoch jedesmal durchgefiihrt 
werden. hy 

d) Die Umgebung der Hin- cil 
terkante. Die Geschwindigkeits- 
verteilung in der Nahe der Hinter- te why 
kante wird durch die Funktionen w, (7) s a. oe, : % ana 
Eid we() bestimmt, ee cintealior Abb. 12. Die Wahl von geese) cole hoher Geschwindigkeits- 
ungefahr entsprechend Abb. 2 voraus- 
gesetzt haben. Nach den vorangehenden Abb. 5, 7, 9 und 11 nimmt dadurch die Differenz der 
 Geschwindigkeit v(x) zwischen Ober- und Unterseite gegen die Hinterkante langsam ab bis zum 
Wert Null. 

Bei den mit Naherungsverfahren berechneten Profilen hat es sich als giinstig erwiesen, wenn man 
bis zur Hinterkante mit nicht verschwindendem Geschwindigkeitsunterschied rechnete. Dies tritt 
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zwar in Wirklichkeit nicht auf, aber die so berechneten Profile zeigten doch bis nahe an die Hinter- — 


kante groBe Geschwindigkeitsunterschiede und lieferten deshalb ebenso weit groBen Auftrieb. 


Oberseite 


a f PEs Ps au 0 | a 7 700 
Abb. 13. Die Folgen von h, und h, im Verlauf v(¢) eines einfachen, unsymme- Abb. 14, Die Folgen von hy und h, im Verlauf »(X, 2) 
trischen Laminarprofils. eines einfachen, unsymmetrischen Laminarprofils, 


Mit unserem Verfahren kann wegen (56) an der Hinterkante kein Geschwindigkeitssprung vor- 
gegeben werden. Wir kénnen aber, ebenso wie es in Wirklichkeit auftritt, den Geschwindigkeits- 
unterschied zwischen Ober- und Unterseite bis fast an die Hinterkante auf einem grofien Wert halten 
und dann rasch gegen Null gehen lassen. Dies gelingt, wenn wir w,(p) und w,(y) etwa gemaB Abb. 12 
wahlen. Wir fiihren bei w,(y) in der Nahe von y = 0 eine rasche Abnahme der GréBe h, ein. Wir 
vermeiden dabei einen Knick und wahlen stetige Tangente. Bei ~ = 0 soll sie nach (121) waage- 
recht sein. Analog erhalt w,(y) bei p = 2 a eine rasche Zunahine der GroBe hg. Damit entstehen bei 
dem Beispiel des einfachen, unsymmetrischen Laminarprofils von b) (Abb. 6 und 7) die in Abb. 13 
und 14 dargestellten Funktionen v(y) und v(X, «;), wenn man gegen vorher nur die Funktionen 
wp) andert. Abb. 14 zeigt genau die beabsichtigte hohe Differenz zwischen Ober- und Unterseite 
bis fast zur Hinterkante. 

Die GréBe h, bedingt auf der Oberseite einen steilen Geschwindigkeitsabfall in der Nahe der 
Hinterkante, also einen Druckanstieg, der leicht zu Strémungsablésung fihrt. Dagegen erzeugt hg 
einen villig ungefahrlichen Druckabfall. Man wird also versuchen, hy klein nad he grok zu machen 
(punktiert in Abb. 14). Eine solche Anderung hat ungefahr die in Abb. 15 skizzierte Anderung 
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AP(p) cur Folge. Nach (34) folgt daraus die in Abb. 15 gleichfalls skizzierte Anderung AQ(¢), 
wie man qualitativ leicht priift. Mit (26) wird also auf der Oberseite (p > 0) die Steigung O griBer, 
auf der Unterseite (p < 0) kleiner. Das entstehende Profil wird an der Hinterkante diinner und 
spitzer. Die Grenze dieser zusatzlichen Zuspitzung ist die Selbstiiberschneidung. 

Wir kénnen also in der Nahe der Hinterkante einerseits durch kleines h, und groBes hg den Druck- 
anstieg der Oberseite vermindern, wobei die Profile hinten dinner und spitzer werden, anderseits 
durch die umgekehrte Tendenz Profile mit zu diinner (oder sich iiberschneidender) Hinterseite 
korrigieren. 

Auch die Folgen einer gleichzeitigen VergréBerung oder Verkleinerung von h, und h, lassen sich 
iibersehen. Im ersten Fall erhalten wir (abgesehen von dem raschen Ausgleich der Geschwindigkeit 

unmittelbar bei der Hinterkante, 


P(g) Ag?) der erfahrungsgem4B keine groben 
gs¢ g 

p g  Nachteile bringt) den hauptsach- 

a lichen Geschwindigkeitsabfall auf 

der Unterseite (vgl. etwa Abb. 14), 

im zweiten auf der Oberseite (etwa 

Abb. 15. Die Anderungen 4 P(p) und AQ(¢) bei VergréBerung der Geschwindigkeit Abb. 7). Wenn alle anderen GréBen 


an der Hinterkante. . 5 , 
gleich bleiben, insbesondere der 


iibrige Verlauf von w,(p) und w,(y), wird sich im ersten Fall K, verkleinern und K, vergréBern, 
im zweiten umgekehrt. Also auch hieritiber kénnen wir ohne Schwierigkeiten verfigen. 

Unser Verfahren ist insgesamt fiir viele wichtige Aufgabenstellungen geeignet. Es folgt an- 
schlieBend noch ein Zahlenbeispiel, an dem die Einzelheiten noch einmal zu verfolgen sind. 


6. Beispiel. Wir belegen unser Verfahren, ohne auf vereinfachte Spezialfalle einzugehen, mit dem 
Beispiel eines unsymmetrischen zweifachen Laminarprofils, dem eine Aufgabenstellung des Segel- 
fluges zugrundeliegt. Dort wird im c,-Bereich zwischen c, = 0,4 und 0,6 extrem niedriger Wider- 
stand verlangt, und gleichzeitig sollen hohe Auftriebsbeiwerte ohne wesentliche Strémungsablésung 
méglich sein. Wir arbeiten mit 2 N = 60 und wahlen 


O10, = 4954 Cy) =i, = 425 
One La,20-, | Yn = Te, = 162°, 
Orit 00 = OLE Oy Og = oa 


Dem Bereich zwischen «, und «, entspricht wegen (100) gerade ungefahr der c,-Bereich zwischen 
0,4 und 0,6. In diesem Bereich steigt dann die Geschwindigkeit auf der Oberseite bis X,, fiir das 
nach (118) X, ~ 90 gilt, auf der Unterseite bis X,; ~ 60. Wir werden also in diesem Bereich hohe La- 
minareffekte erhalten. Die c,-Werte unter 0,4 interessieren nicht, dagegen diejenigen iiber 0,6 noch 
sehr weit. Wir verzichten deshalb mit «, =o, auf den Abbau der Saugspitze unterhalb von «; und 
nehmen nur mit dem hohen «,-Wert auf die hohen ¢, Riicksicht. Bei ihnen fallt die Geschwindigkeit 
auf der Oberseite schon von X, ~ 5 an ab. Die Grenzschicht wird also schon frith turbulent und 
hat in der Nahe der Hinterkante schon groBe Dicke. Damit dadurch keine wesentlichen Ablésungen 
entstehen, haben wir X, gréBer als X; gewahlt und den hauptsachlichen Geschwindigkeitsabfall, 
der fiir die SchlieBung des Profils natiirlich vorhanden sein mu, auf die Unterseite verlegt. Auf der 
Oberseite fiihren wir nur mit 


1 — cos 9, 


w,(p) =1—hy een oye (134) 


h, = 0,36, einen kurzen Druckanstieg ein, wie er fiir den Druckausgleich an der Hinterkante nach 
Ziff. 5 d) erforderlich ist. Fiir die Form des Druckanstiegs auf der Unterseite wahlen wir 


we(?) = (1 L og cos ~ — cos a {1 he ee == GOS a. | 


1 — cos 9; 1 — cos g¥ 


(135) 


mit 

(ji==1 fir p= 9F 
und % = 11,26, h, = 0,2025, p¥ =1;, = 330°. Der erste Term gibt in groBen Ziigen den Geschwin- 
digkeitsabfall wieder, wie ihn F'. X. Wortmann} empfiehlt. Der zweite Term in (135) gibt den Anteil 


der Unterseite zum Druckausgleich an der Hinterkante. Er wurde mit einem relativ kleinen h,-Wert 
relativ klein gewahlt. Entsprechend wird der Anteil der Oberseite zum Druckausgleich gréBer 


1 FuBnote 2 vonS. 48. 
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werden, was nach Ziff. 5d) eine véllige Hinterkante erwarten laBt. Die Druckausgleichsglieder 
wurden mit Riicksicht auf die Rechnung multiplikativ eingefiihrt, wodurch sich die Bedeutung der 
GréBen h, und h, gegeniiber Abb. 12 unwesentlich verschiebt. 


Der Rechnungsgang verlauft nun genau entsprechend Ziff. 4. Die W.;, W,; von (61) und (62) 
lassen sich geschlossen auswerten. Es gilt hierfiir 


y 2% 
= y | c “t 
cos ta 1 eae cos Pp cos y Zs : cos ie Bac cos wy i” = 
| yp ie if Y cos g In|1 nt See i lp = —Ap—siny 
0 2a—p 
ares Tene j@—1 
+ 2 yA?— 1 are a tg eh (136) 


1 1 
A = & 1) cosy. 


Fir negatives m wird A < 1, so daB man beim letzten Glied in bekannter Weise auf ArLq ibergehen 
mu8B. AuBerdem ist 


yp 2a 
te | cos Y — cosy Be seit cos p — cos p 
| sin @ In y [Sree ‘i | dy = sing In |1 —m eon 
22—p 
= — (1 — cos y) f te | In (1 m)| : (i37) 
Damit wird, wenn man jeweils fiir y die @; und ~;, fiir m die V/h; und x, eingefiihrt, nach (61) und (62 
J ee? Ys 6 C1Dg 
W., = 0.197262 , W,, = =.0,034860 , In w,(0) = — 0,446287, 
Wo 0427 715,, W,, = — 0,305324, In w,(2 a) = — 0,350185 
und nach (79) bis (81) 
D, = 0,124048, D, =—0,245954, D3 = — 0,062935. 


In der transzendenten Gleichung (82) wird 
a = 0,242120, 6b =—0,801684, c= 0,951164, d = 0,034488 , 
woraus 
= 193> 58! 48 == 3.385590 


folgt. Dieser Wert wird zweckmaBig sehr genau bestimmt, weil er sehr empfindliche Einfliisse 
austiibt. Mit ihm wird entsprechend (88) 


J.— 1% = — 0,544720 , I-20 256317. ; J; = — 0,071977 
und nach (83) und (84) 
K, = 0,90074 , Ke = 0594239", 
Der Druckanstieg in (135) und die dadurch im wesentlichen festgelegte Profildicke erweist sich als 
genigend groB, so daB bei unserem kleinen Wert S (d, —d,) = 0,825° noch geniigend Spielraum 


fiir die vordere Laminarstrecke blieb. Mit den angegebenen Werten fiir y, und a; wurde gerade 
K, ~ lerreicht. Dies gelang nicht sofort, da aber der Einflu8 von g, auf Kg leicht zu iibersehen ist 
(kleineres g, ergibt gréBere Profildicke und gréBeres Ky), konnte Pz, leicht festgelegt werden. 

Nun berechnen wir die P, nach (90), v, nach (98), v*(y) nach (52) und evtl. einige v(p, d;) nach (47). 
An ihnen priifen wir nach, ob durch die v; und K; keine unbeabsichtigten Wirkungen auftreten. 
Wir kénnen auf die tabellarische Wiedergabe dieser Rechnungen verzichten, da wir am Ergebnis 
alles wesentliche iibersehen. Auch der Ubergang (94) auf die Q, sowie die Berechnung der x,, y, 
bereitet keinerlei Schwierigkeiten. Wir arbeiten in (101) und (102) mit B = 8°. Dieser Wert kann 
auf Grund des in den x, und y, schon enthaltenen Profilsteigungswinkels Of). leicht so gewahlt 
werden, daB das Profil im X, Y-System ungefahr in der iiblichen Lage erscheint. Die Fehler ¢, 
und ¢, nach (107) werden 


&, = — 9,00082 , i= 0,00149 , 
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sind also bei einer Lange von 1, ~ 0,9 geniigend klein. Nach dem Ausgleich derselben und dem 
Ubergang auf X und Y ergeben sich die in Tabelle 3 dargestellten Profilkoordinaten, sowie die gleich- 
falls angegebenen Werte von 
ate kal (138) 
L 2% dx 
nach (100) und von ¢,,) nach (42) und (43). Der Wert von c¢,,) weist auf ein relativ hohes kopf- 
lastiges Moment hin. Dieser Nachteil spielt aber bei den schlanken Fliigeln von Segelflugzeugen 
gegeniiber den erzielten Vorteilen keine Rolle. Die Langsstabilitat ist auch so noch leicht zu 
erreichen, und die Drehsteifigkeit von Laminarfliigeln ist ebenfalls immer ausreichend, da sie 
infolge der erforderlichen hohen Oberflachengiite ohnehin meist als ganze Schalen ausgebildet 
werden. 
v 


06,=13,2° 
6,=4,95° au Oberseitex 
a eS a ee ee 


0,-3.3) 


7 Unterseite 


0 50 100 
Abb. 16, Das Profil der Beispielrechnung mit seinen Verteilungen v(X, a). 


In Abb. 16 sind fiir die drei «;-Werte die Verteilungen v(X,«;) angegeben, die fiir die Be- 
rechnung der Grenzschicht nétig sind. Die Linien v*(~) sind verstarkt gezeichnet. AuBerdem ist 
noch das berechnete Profil wiedergegeben. 

Das Beispiel gibt einen Begriff von den in den zweifachen Laminarprofilen enthaltenen Entwurfs- 
méglichkeiten. Bis hinauf zu « = a; = 13,2°, was wegen (100) c, ~ 1,6 entspricht, ist die héchste 
auftretende Geschwindigkeit vn, = 1,98. Dabei bleibt die Profildicke so klein, daB zwischen «, 
und «;, wegen (100) zwischen cg = 0,4 und c, = 0,6, auf der Oberseite bis X, = 90, auf der Unter- 
seite bis X; = 50 Druckabfall vorhanden ist. Der entstandene Wert von X, hat sich gegeniiber (118) 
etwas verschoben. Diese Verschiebung ist bei allen gewélbten Profilen auf der Unterseite in gewis- 
sem Ma} vorhanden und war auf Grund ahnlicher Beispiele vorauszusehen. An der Hinterkante 
ist der Druckanstieg auf der Oberseite wesentlich gréBer als der Abfall auf der Unterseite. Dadurch 
ergibt sich fast genau ein endlicher Hinterkantenwinkel (Abb. 16), obwohl er mathematisch ver- 
schwindet. Man kénnte zwar h, noch etwas kleiner machen und h, gréBer, also Druckanstieg auf der 
Oberseite sparen. Viel ist hierbei jedoch nicht zu erreichen, da die dabei nach Ziff. 5 d) schlanker 
werdende Hinterkante rasch an die Grenze der Selbstiiberschneidung kommt. Schon das darge- 
stellte Beispiel hat keine sehr villige Hinterkante, eine ahnliche Rechnung zeigte schon eine kleine 
Selbstiiberschneidung, ohne da® der Druckanstieg auf der Oberseite entscheidend yermindert 
werden konnte. Die Auftriebserhéhung in der Nahe der Hinterkante nach Ziff. 5d) bringt also 
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potentialtheoretisch immer einen Druckanstieg auf der Oberseite mit sich, der durch Grenzschicht- 
und Ablésungseinfliisse abgebaut wird. Dies hat sich bei den niherungsweise berechneten Profilen! 
nicht nachteilig ausgewirkt. Man kann also vielleicht bei weiteren Beispielrechnungen h, gréBer 
und hg Kleiner wahlen und erhalt dadurch je nach Anforderung noch weniger schlanke Hinterkanten. 
Ebenso wie an der Hinterkante lat sich die Geschwindigkeitsverteilung auch noch beziiglich 
anderer Hinzelheiten sehr feinfiihlig variieren. Dies ist jedoch Sache weiterer Bec eenacacen: 
die tiber den Rahmen des grundlegenden Verfahrens hinausgehen. i 


Tabelle 3. Die Koordinaten des als Beispiel berechneten Profils. 


x Y x Y oe Y 
100,00 0,00 29,51 10,02 V7.9 1 28:27 
99,67 0,13 25,06 9,77 21,44 | —3,66 
98.73. | 0,53 20,87 9,43 25,64 | —3,97 
97,29 1,20 16,98 9,00 30,04 | —4,15 
95,45 2,08 13,43 8,48 34,58 | —4,21 
93,28 3,05 10,26 7,87 39:90 9) ==4.10 
90,78 4,02 7,52 7,19 43,84 | —3,80 
87,92 | 4,92 5,24 6,42 TYRE al an Oy 
84567 | 2-°5,75 3,42 5,50 5a 0) eet TO 
81,05 6,53 1,97 4,48 58,65 | —0,34 
77,09 7,25 0,90 3,45 64,40 0,82 
72,84 7,91 0,23 2,47 70,18 1,71 
68,35 8,49 0,01 1,65 75,82 2,27 
63,66 9,01 0,41 0,97 81,16 2,53 
58,82 9,44 1,43 0,31 86,06 2,48 
53,88 9,78 2,99 —0,35 90,39 2,17 
48,89 10,02 5,05 == ,0% 94,00 1,63 
ZOE ei 1017 7,51 —1,64 96,74 0,99 
38,98 | 10,22 10,50 == 9 94, 98,60 0,45 
34,17 10,17 13,82 —2,79 99.66 0,11 

100,00 0,00 

2 = 1,1097 Cp = —0,2097 

ut 


7. Zusammenfassung. Es werden einfache Formeln abgeleitet fir die direkte Bestimmung eines 
Fligelprofils aus seiner Geschwindigkeitsverteilung, die gleichzeitig leicht nachzupriifen erlauben, 
ob eine vorgegebene Geschwindigkeitsverteilung iitberhaupt zu einem sinnvollen, geschlossenen Profil 
gehort. Die Einfachheit der Formeln erméglicht einen speziellen Ansatz fiir die Geschwindigkeits- 
vorgabe, der auf eine praktisch wichtige Verallgemeinerung der bisher bekannten Laminarprofile 
fiihrt. Das Berechnungsverfahren wird bis in alle numerischen Einzelheiten ausgebaut und an 
einem Beispiel belegt. 


(Eingegangen am 5. Mai 1956.) 
Anschrift des Verfassers: Dr. Richard Eppler, Stuttgart W, Leibnizstr. 84. 


1 Vgl. beispielsweise NACA-Report 824, FuBnote 1 von S. 38. 
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Thermodynamik und rheologische Probleme 


Von H. Ziegler 


1. Einleitung. Die vorliegende Arbeit stellt einen Versuch dar, die Thermodynamik zur 
Abklarung gewisser rheologischer Probleme heranzuziehent. Dabei soll insbesondere folgendes 
gezeigt werden: 

1. Die sogenannte viskoelastische Analogie, die von T. Alfrey? entdeckt und seither in verschie- 
denen Richtungen verallgemeinert, aber stets nur fiir bestimmte (insbesondere isotrope) Kérper 
bewiesen worden ist, laBt sich bei Beschrankung auf lineare rheologische Gesetze thermodynamisch 
als Sonderfall der Onsagerschen Theorie linearer irreversibler Prozesse® in voller Allgemeinheit 
begriinden. 


2. Die von N. J. Hoff* aufgezeigte Méglichkeit, die viskoelastische Analogie auch auf bestimmte 
nichtlineare Kérper auszudehnen, legt eine Erweiterung der Onsagerschen Theorie auf nichtlineare 
irreversible Prozesse nahe. Ob eine solche Verallgemeinerung physikalisch haltbar ist, wird wohl 
nur mit den Mitteln der statistischen Thermodynamik abgeklart werden kénnen. Hier soll nur 
gezeigt werden, dafs aus dieser Erweiterung die unbeschrankte Giiltigkeit der viskoelastischen Ana- 
logie folgen wiirde. 


3. Die Misessche Theorie des plastischen Potentials> weist in ihrer Begriindung noch einige 
Liicken auf. Diese lassen sich schlieBen, wenn man annimmt, daB der idealplastische Kérper als 
Grenzfall eines viskosen Kérpers aufgefaBt werden darf, auf den sich die viskoelastische Analogie 
anwenden 1aBt. 


2. Rheologische Grundlagen. Wird die Deformation eines Kérpers durch die Formainderungen 
x, beschrieben, und ist die von den auBeren Kraften bei der Anderung dx; der Deformation geleistete 


Arbeit ® 
dA XG dee (2.1) 


dann kénnen die rechterhand auftretenden Beiwerte als Spannungen bezeichnet werden. Die 
verschiedenen Kérper, die in der Rheologie untersucht werden, unterscheiden sich in der Art des 
Zusammenhangs zwischen den Deformationen und den Spannungen, und die Art dieses Zusammen- 
hangs wird das rheologische Gesetz des betreffenden Kérpers genannt. Im Falle eines Raum- 
elementes werden als Formanderungen zweckmabig die Komponenten ¢,, des Verzerrungstensors 
gewahlt; die Spannungen sind dann die Komponenten o;; des Spannungstensors, und das rheolo- 


gische Gesetz vermittelt den Zusammenhang zwischen den Verzerrungs- und den Spannungskom- 
ponenten. 


Das Typische eines rheologischen Gesetzes tritt vielfach schon im eindimensionalen Fall, 
d.h. dann zutage, wenn die Deformation durch eine einzige Formanderung x beschrieben werden 
kann. Es 1a8t sich besonders anschaulich mit Hilfe eines rheologischen Modells darstellen. 


So hat das rheologische Gesetz des elastischen oder Hooke-Kérpers im eindimensionalen Fall 


(linearer Spannungszustand, Zustand reinen Schubes usw.) die einfache Form des Hookeschen 
Gesetzes 


Xt eee (2.2) 


wobei mit x* die (elastische) Formanderung, mit X die Spannung und mit a® die einzige elastische 
Konstante (Elastizitats-, Schubmodul usw.) bezeichnet ist. Als eindimensionales Modell des Hooke- 


Kérpers hat sich die Schraubenfeder (Abb. 1) mit der Verlangerung x°, der Belastung X und der 
Federkonstanten a° eingebiirgert. 


1 Der Verfasser ist den Professoren Dr. W. Prager, Dr. H. Wéiffler und Dr. K. Bleuler, sowie seinem Assi- 
stenten, Herrn Dipl.-Math. A. Schopf, fiir wertvolle Anregungen zu Dank verpflichtet. 
> T. Alfrey, Quart. appl. Math. 2 (1944) S. 113. 
° L. Onsager, Phys. Rev. 37 II (1931) S. 405 und 38 II (1931) iS: 2265. 
* N. J. Hoff, Quart. appl. Math. 12 (1954) S. 49. 
: oe v. Mises, Z. angew. Math. Mech. 8 (1928) S. 161. 
iu 


r untenstehende Zeiger soll hier und im folgenden stets die Summationsregel angewandt werden. 


oO 
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Ahnlich lautet das theologische Gesetz des Newton-Kérpers im eindimensionalen Fall 
x: — q" an. (2.3) 


wobei jetzt x" die (unelastische) Formanderungsgeschwindigkeit und a" die Zahigkeits- 
konstante ist. Das zugehérige Modell ist in Abb. 2 wiedergegeben und besteht aus einem mit vis- 
koser Fliissigkeit gefiillten Zylinder, in dem ein durchlécherter Kolben gleitet. 

Das eindimensionale Modell des Maxwell-Kérpers (Abb. 3) wird dadurch erhalten, daB man 
sich eine elastische Einheit gema8 Abb. 1 und eine viskose Einheit gemaB Abb. 2 hintereinander- 
geschaltet denkt. Die Formianderung setzt sich hier gemaB 


je st aed all gl 


aus einem elastischen Anteil x° und einem unelastischen x“ zusammen, und da diese durch (2.2) 
und (2.3) gegeben sind, lautet sein rheologisches Gesetz im eindimensionalen Fall 


x 
Abb. 1. Modell des Abb. 2. Modell des Abb. 3. Modell des Abb. 4. Modell des 
Hooke-Kérpers. Newton-KG6rpers. Maxwell-Kérpers. Kelvin-Kérpers. 


Ahnlich erhalt man durch Parallelschalten einer elastischen und einer viskosen Einheit (Abb. 4) 
das eindimensionale Modell des Kelvin-Kérpers. Sind x seine Formanderung und X°, X" die Span- 
nungen in den beiden Zweigen, dann gilt 

Xa Xe XK, (2.6) 
und da zufolge (2.2) sowie (2.3) 
Xt = 4 und X= a & (257) 
ist, lautet das rheologische Gesetz im eindimensionalen Fall 
B= Ox ax. (2.8) 

Mit dem Hooke-Kérper ist der Kelvin-Kérper als fest zu bezeichnen, mit dem Newton-Kéorper 
der Maxwell-Kérper als fliissig. Die auffalligsten Eigenschaften des Maxwell-Korpers sind das 
Kriechen (namlich die fortschreitende Deformation unter einer konstanten Spannung) und die 
Relaxation (der Abbau der Spannung bei festgehaltener Deformation), wahrend der Kelvin- 
Kérper als elastischer Kérper mit Retardation aufgefaBt werden kann (nach einem raschen 
Lastwechsel stellt sich die nach dem Hookeschen Gesetz der neuen Spannung entsprechende Defor- 
mation nur asymptotisch ein). All diese Eigenschaften sind yon M. Reiner’ ausfihrlich diskutiert 
worden, wie auch die Méglichkeit, durch weitere Kombination der behandelten Elemente noch 
wirklichkeitsgetreuere Modelle fiir bekannte Kérper zu erhalten. _ 

Der elastisch-idealplastische oder Prandtl-Reuf-Kérper besitzt im eindimensionalen Fall das 
in Abb. 5 wiedergegebene Modell, das aus einer Schraubenfeder und einem starren Kérper auf 
rauher, trockener Unterlage besteht. Auch hier gilt die Zerlegung (2.4), wobei die elastische Form- 
anderung durch (2.2) gegeben ist, wahrend die unelastische (plastische) Formanderungsgeschwin- 


1M. Reiner, Twelve Lectures on Theoretical Rheology, Amsterdam 1949. 
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digkeit aus 


ppberiite pi 


folgt, wenn k die sogenannte FlieBgrenze und A eine beliebige nichtnegative GréBe ist. Der 
Prandtl-ReuB-Kérper hat also fiir Spannungen, die unterhalb der FlieBgrenze liegen, die Kigen- 
schaften eines festen (elastischen) Kérpers und fiir Spannungen an der FlieBgrenze (hohere Span- 
nungen kommen nicht vor) diejenigen einer (nichtviskosen) Fliissigkeit. 

Die vorstehenden Uberlegungen geben zu folgenden Bemerkungen AnlaB: 

1. Beim Maxwell- und beim Prandtl-ReuB-Kérper 1aBt sich die Formanderung gemaB (2.4) in 
einen elastischen und einen unelastischen Anteil zerlegen. Dabei wollen wir vorderhand unter dem 
Vorbehalt, diese Definition in Ziff. 4 durch eine thermodynamische zu ersetzen, unter dem elasti- 
schen denjenigen Anteil verstehen, der sich bei einer Entlastung zuriickbildet, wahrend der unela- 
stische der auch nach der Entlastung zuriickbleibende Anteil ist. Der Hooke- und der Newton- 
Korper stellen Grenzfalle dar, in denen der eine der 
beiden Anteile verschwindet. Beim Kelvin-Kérper 
dagegen (wie auch beim Poynting-Thomson- sowie 
beim Burgers-Korper) ist die Unterteilung in ela- 
stische und unelastische Formanderung in dieser ein- 

VAL LLTMM IIL fachen Weise nicht méglich. Kérper der zuletzt ge- 
Abb. 5. Modell des Prandtl-Reup-Kérpers. nannten Art sollen daher im folgenden ausgeschlos- 
sen bleiben. 

2. Von den oben beschriebenen Kérpern gehorchen alle bis auf den Prandtl- Reuf-Kérper rheo- 
logischen Gesetzen, die sich zunachst im eindimensionalen, dann aber auch im mehrdimensionalen 
Fall durch lineare Beziehungen zwischen den Formanderungen, den Spannungen und den zeit- 
lichen Ableitungen dieser GréBen beschreiben lassen. Beim Prandtl-Reufi-Kérper dagegen ist dies 
nicht méglich; sein rheologisches Gesetz ist demnach als nichtlinear zu bezeichnen. Selbstverstand- 
lich kénnen auch die tibrigen Kérper so verallgemeinert werden, daB sie nicht mehr linear sind. 
Man hat zu diesem Zweck nur die elastischen Konstanten ihrer elastischen Elemente als Funk- 
tionen der elastischen Formanderungen oder/und die Viskositétskonstanten ihrer viskosen Ele- 
mente als Funktionen der unelastischen Formanderungsgeschwindigkeit anzusetzen. Die prak- 
tisch wichtigste Abweichung vom linearen Gesetz liegt indessen wohl beim plastischen Kérper vor. 
Dieser soll im folgenden in den Vordergrund gestellt werden, nicht zuletzt deshalb, weil er ur- 
spriinglich den Ausgangspunkt fiir diese Untersuchung bildete. 

3. Die sogenannte viskoelastische Analogie beruht in ihrer einfachsten Form auf der Ahn- 
lichkeit im Aufbau der rheologischen Gesetze (2.2) beim Hooke- und (2.3) beim Newton-Kérper. 
Man zieht hieraus zunachst den trivialen SchluB, daB sich im eindimensionalen Fall beim Newton- 
Kérper die (unelastische) Formanderungsgeschwindigkeit in der gleichen Weise aus der Spannung 
berechnen lat wie beim Hooke-Kérper die (elastische) Formanderung. Interessant wird die Ana- 
logie freilich erst bei mehrdimensionalen Kérpern (d. h. dann, wenn die Deformation nur durch 
mehrere Formanderungen x; beschrieben werden kann). Daf sie auch dann noch gilt, ist aber 
keineswegs selbstverstandlich und la8t sich nur dadurch begriinden, daB man sich 

a) nach dem Vorgang von T. Alfrey!, H. S. Tsien?, H. H. Hilton? und N. J. Hoff* auf Kérper 
mit mehr oder weniger speziellem rheologischem Gesetz beschrankt, 

b) der Mittel der Thermodynamik oder 

c) bestimmter Modellvorstellungen 
bedient. 

4. Um beim Hooke-Kérper auch im mehrdimensiona'en Fall iiber ein einfaches Modell zu ver- 
fiigen, liegt es nahe, an ein System elastischer Einheiten gema® Abb. 1 zu denken. Diese Vor- 
stellung, die eine genaue Kenntnis des molekularen Aufbaus ersetzt, ist naturgem4B einigermaBen 
willkiirlich. Kine elementare Untersuchung, die entweder rein mechanisch (Ziff. 3) oder thermo- 
dynamisch (Ziff. 4) durchgefiihrt werden kann, zeigt aber, daB diese Modellvorstellung sowohl im 
linearen wie im allgemeineren Fall zweckmaBig ist. Baut man in analoger Weise das Modell des 


1s, FuBnote 2, S. 58. 
* H. S. Tsien, Quart. appl. Math. 8 (1950) S. 104. 
_° H. H. Hilton, Analytical Studies of Thermal Stresses in Media Possessing Temperature-Dependent 
Viscoelastic Properties, Wright Air Development Center, Tech. Rep. 53—322 (1953). 
* s, FuBnote 4, S.58; vgl. auch N. J. Hoff, Stress Distribution in the Presence of Steady Creep, Proc. 
Conference on High-Speed Aeronautics, Polyt. Inst. of Brooklyn, S. 271, New York 1955. s 
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mehrdimensionalen Newton-Kérpers aus viskosen Einheiten gemaB Abb. 2 auf, dann ist natiirlich 
die viskoelastische Analogie fiir beliebig viele Dimensionen, ja sogar fiir nichtlineare Korper ganz 
allgemein gesichert. Diese Modellvorstellung ist indessen noch ein Stiick willkiirlicher als die- 
jenige des Hooke-Kérpers; ihre allgemeine Brauchbarkeit l4Bt sich namlich mechanisch iiberhaupt 
nicht und thermodynamisch nur im linearen Fall (Ziff. 5) einsehen. 


3. Mehrdimensionale Kérper. Die mit den letzten Bemerkungen angedeuteten Probleme treten 
sofort klar zutage, wenn man die rheologischen Gesetze von Ziff. 2 auf mehrere Dimensionen zu 
ibertragen versucht. 

Im Falle eines Gesetzes, das im eindimensionalen Fall linear ist, la8t man sich von der Forde- 
rung leiten, daB auch das mehrdimensionale Gesetz linear sein soll. Man verallgemeinert daher das 


rheologische Gesetz (2.2) des Hooke-Korpers mit 


Me Sg X. (3.1) 

Aus physikalischen Griinden darf man annehmen, da sich (3.1) umkehren 1a8t, d. h. daB 
a® = |a®| + 0 ist. Die Umkehrung lautet 

a= bi, X; mit bf, =—, (3.2) 


wobei A; der Kofaktor von aj, in der Determinante a’ ist. 
Analog erhalt man fiir das rheologische Gesetz (2.3) des Newton-Kérpers die Verallgemeinerung 


sae (3.3) 
mit der Umkehrung 


Cee x . uo 
x — bX; mit by, 


(3.4) 
Beim Maxwell-Kérper kommt mit (3.2) und (3.4) 

Spee a eter 6) XG, Sere OX (3.5) 

und fiir einen verallgemeinerten Newton-Kérper, dessen rheologisches Gesetz im eindimensionalen 


Fall durch 
X=cx (3.6) 


gegeben ist, wird man im zweidimensionalen Fall 


hohe ie 1 ne A ee eh 
Xy = Cy Xf + Cy XE Xq + Cy Xy XQ oe | 


2 23 
2 C82 


(3.7) 


Xy = 5X} + Cg Xf Ky + C7 X 
ansetzen. 

Damit ist aber das Ubertragungsproblem noch nicht erledigt. In all diesen Fallen ist namlich 
noch zu untersuchen, ob die Beiwerte, die im mehrdimensionalen Gesetz auftreten, voneinander 
unabhangig oder durch irgendwelche Beziehungen miteinander verkniipft seien. Dabei sollen 
Symmetrien, die durch den speziellen Aufbau oder die Isotropie des betrachteten Kérpers bedingt 
und fiir die in Ziff. 2 unter Bemerkung 3 a) zitierten Beweise der viskoelastischen Analogie wesent- 
liche Voraussetzungen sind, unberiicksichtigt bleiben. 

Im Falle des Hooke-Kérpers kann man bekanntlich die Spannungen X; gemaB 


oD 
Lae a 3.8 
aus einem Potential, namlich der (positiv definiten) Formanderungsenergie 
1 
cD) = oe ai, xe oy (3.9) 


herleiten, deren Beiwerte ohne Beschrankung der Allgemeinheit noch der Symmetriebedingung 
Ai, = Aig (3.10) 


unterworfen werden kénnen. Aus (3.8) bis (3.10) folgt dann (3.1); gleichzeitig ist aber festgestellt, 
daB die Koeffizientenmatrix (a%,) von (3.1) symmetrisch ist. Auch beim verallgemeinerten nicht- 
linearen Hooke-Kérper gilt die Darstellung (3.8), die sich auch hier auf die moglichen Ansatze ein- 
schrankend auswirkt. Beim Newton-Kérper dagegen liegt mechanisch kein Grund fiir eine ahnliche 
Symmetrie vor, und damit erscheint die viskoelastische Analogie in einem auferst fragwirdigen 


Lichte. 
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Im Grunde dasselbe Problem, wenn auch in veranderter Gestalt, stellt sich bei der Erweiterung 
des rheologischen Gesetzes (2.9) des Prandtl-Reuf-Kérpers auf mehrere Dimensionen. Bekanntlich 
wird diese Erweiterung durch die Misessche Theorie des plastischen Potentials! geleistet, die sich 
in der verallgemeinerten Pragerschen Gestalt? wie folgt formulieren laBt: 


Ist die FlieBgrenze im Spannungsraum X; eine geschlossene, konvexe, den Ursprung enthaltende 
regulare F'lache mit der Gleichung 


F(X) =c, (3.11) 


in der c eine Konstante und die Funktion F so gewahlt ist, daB sie im Inneren der FlieSflache 
kleiner als ¢ ist, dann sind die plastischen Formanderungsgeschwindigkeiten xj, welche zu einem 
gegebenen Spannungszustand X, an der FlieBgrenze gehéren, durch 

ewe oF 

x, =A aX, (3.12) 
gegeben®, wobei A ein beliebiger nichtnegativer Proportionalitatsfaktor ist. 

Diese Theorie, in der bemerkenswerterweise wieder ein Potential — das sogenannte plastische 
Potential F(X,) — eine Rolle spielt, wird insbesondere auch auf starr-plastische Kérper und auf 
solche mit Spannungsverfestigung angewandt; im letzten Falle muB man lediglich zulassen, daB 
die Funktion F' auch von der plastischen Vorgeschichte abhangen kann. Ferner ist die Theorie 
durch W. T. Koitert und W. Prager® auf irregulare, d.h. aus mehreren regularen Teilstiicken 
F(X;) = 0 zusammengesetzte FlieBflachen verallgemeinert worden. 

Die Misessche Theorie weist indessen zwei Liicken auf. Erstens postuliert sie eine konvexe 
FlieBflache. Unzahlige Anwendungen haben in der Tat stets konvexe FlieBflachen ergeben, aber 
ein physikalischer Grund hiefiir hat sich bis heute nicht aufzeigen lassen. Zweitens entbehrt auch 
die Beziehung (3.12) der physikalischen Begriindung und stellt damit nur eine Hypothese dar. 
R. v. Mises hat sich denn auch darauf beschrankt, die formale Zulassigkeit des Ansatzes (3.12) so- 
wie seine Brauchbarkeit in einigen Sonderfallen nachzuweisen. W. Prager’ hat gezeigt, da} der 
Ansatz (3.12) eine hinreichende Bedingung dafiir ist, daB die Spannungsinkremente in einem end- 
lichen Kérper mit Spannungsverfestigung eindeutig durch die Inkremente der Oberflachenkrafte 
gegeben sind. Auch R. Hill’ beruft sich zur Begriindung von (3.12) auf die erst damit gegebene 
Méglichkeit, gewisse Eindeutigkeits- und Variationsprobleme fiir endliche Kérper zu formulieren. 
Es ist klar, daB diese Argumente wie auch die Erfahrung die Misessche Theorie weitgehend sichern. 
Sie tun dies aber auf indirektem Wege und dirfen nicht dariiber hinwegtauschen, daB eine direkte 
Begriindung noch fehlt, von der aus umgekehrt (und analog wie dies in der Elastizitatstheorie 
getan wird) die Eindeutigkeits- und Variationssdtze erschlossen werden kénnten. 

Es liegt nahe, fiir die Beantwortung der hier aufgeworfenen Fragen die Thermodynamik her- 
anzuziehen. Der zweite Hauptsatz diirfte dazu freilich nicht ausreichen, da aus ihm nicht die Rich- 
tung erschlossen werden kann, in der ein gegebener Vorgang verlauft. sondern nur sein Richtungs- 
sinn. Die Thermodynamik verfiigt aber in den Onsagerschen Beziehungen® und der hierauf beruhen- 
den Theorie® iiber ein Mittel zur Bestimmung der Ablaufrichtung irreversibler Vorgange. 


4. Thermodynamische Grundlagen. Bei den Kérpern, auf die wir uns gema Bemerkung 1 in 
Ziff. 2 beschranken, kénnen die Formanderungen x, in elastische Anteile xf und unelastische x¥ 
zerlegt werden, wobei wir jetzt aber von der in Ziff. 2 verwendeten Abgrenzung keinen Gebrauch 
machen und uns vorbehalten, diese unten neu zu definieren. Ist } die Temperatur, so kénnen die 
GroBen x;, x, J als unabhangige thermodynamische Zustandsvariablen des Kérpers aufgefabt 
werden. Sollen sie seinen thermodynamischen Zustand vollstandig festlegen, dann mu8 man sich 
im folgenden auf geniigend langsame (quasistationaére) Zustandsanderungen oder auf Raumele- 
mente beschranken, so daB zum Beispiel stets mit einem vdlligen Temperaturausgleich tiber den 
ganzen Kérper gerechnet werden kann. 

ls, FuBnote 5, S. 58. ; 

2 W. Prager, Probleme der Plastizitatstheorie, S. 18 ff., Basel 1954. 

3 Die hier und im folgenden auftretenden Funktionen werden, soweit nichts anderes vermerkt ist, als stetig 
differentiierbar vorausgesetzt. 

* W. T. Koiter, Quart. appl. Math. 11 (1953) S. 350. 

° W. Prager, J. appl. Mech. 20 (1953) S. 317. 

° W. Prager, J. appl. Phys. 20 (1948) S. 235. 

? R. Hill, The Mathematical Theory of Plasticity, S. 50, Oxford 1950. 

8s, FuBnote 3, S. 58. 


° Vel. auch S, R. de Groot, Thermodynamics of Irreversible Processes, Amsterdam 1952. 
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Die Arbeit der auBeren Krafte bei einer Zustandsanderung dx, dx, d9 ist durch einen Aus- 
druck der Form (2.1), d. h. durch 


dA = X, (dxt + dx) (4.1) 


gegeben, und die GréBen X; kénnen wieder als Spannungen bezeichnet werden. Ist ferner dQ die 
zugefiihrte Warmemenge, dann gilt nach dem ersten Hauptsatz der Thermodynamik 


A ee dO d (4.2) 


: aT Sas Stef : : ; noe : 
wobei U(x;,, xf, 0) eine eindeutige Zustandsfunktion, nimlich die innere Energie des Kérpers 


ist. Aus (4.1) und (4.2) folgt sodann 


a aU . , (au as: 
0 =4d0 dA = teas x; due + Ee ms x.) deg + So dd. (4.3) 


Ferner ist nach dem zweiten Hauptsatz der Thermodynamik 
dQ@<¢ds, (4.4) 


wobei S(xj,, «;, ?) wiederum eine eindeutige Zustandsfunktion, namlich die Entropie des Kérpers 
ist und das Gleichheitszeichen fiir reversible, das Ungleichheitszeichen fiir irreversible Zustands- 
anderungen gilt. 

Man kann nun die Abgrenzung zwischen den xf und den x durch die Definition treffen, daB die 
elastischen Formanderungen reversible, die unelastischen irreversible Zustandsanderungen dat- 
stellen sollen. Dann ist auf Grund von (4.3) und der Beziehung (4.4), in der jetzt fir dx¥ = 0 das 
Gleichheitszeichen gilt, 


aU Ui 
ie oe: X,] dx, + 5.49 = 9 dS (4.5) 
und daher 
aS). 1o7aU aS 1 aU 
Ta eee rere eo) 


wahrend iber die partiellen Ableitungen 0S/0x} noch nichts ausgesagt werden kann. 
Fihrt man wie tiblich mit 


F=U—d#8S (4.7) 


eine weitere eindeutige Zustandsfunktion F(«j, xi, 0) ein, die als freie Energie des Korpers 
bezeichnet wird, dann folgt aus (4.6) und (4.7) 

oF oF 

De = on S cae ao . 


La One ? 


(4.8) 


Die abhangigen Zustandsvariablen X; und —S lassen sich also als partielle Ableitungen der freien 
Energie nach den zugehérigen unabhangigen Variablen gewinnen, und damit rechtfertigt sich fir F 
auch die Bezeichnung als thermoelastisches Potential. Zufolge (4.8) ist das thermoelastische 
Verhalten des Korpers (seine elastischen Konstanten, thermischen Formanderungskoeffizienten und 
spezifischen Warmen) durch die Funktion F vollstandig bestimmt und mit ihr im allgemeinen auch 
von den unelastischen Formanderungen x; abhangig. 

Es ist in der Rheologie iiblich, sich auf Kérper zu beschranken, deren thermoelastisches Ver- 
halten von den unelastischen Formanderungen unabhangig ist. In diesem Falle ist die freie Energie 
nach (4.8) bis auf eine additive Funktion f(x}) eine Funktion F(x;, 0) der elastischen Formande- 
rungen und der Temperatur allein. Ferner ist dann nach der zweiten Beziehung (4.8) die Entropie 
von der Form S(xé, 9), d. h. ebenfalls von den xj unabhangig, und nach (4.7) gilt auch fiir die 
innere Energie bis auf die additive Funktion f(xj) die Darstellung U(x;, 7), Geht man die Be- 
ziehungen (4.1) bis (4.8) einzeln durch, dann stellt man fest, daB sich die additive Funktion f(x;) 
nur in (4.3) bemerkbar macht. DaB sie dort auftritt, hat zur Folge, das das unelastisch-thermische 
Verhalten (z. B. die spezifischen Warmen des unelastischen Kérpers) im Aligemeinen von den 
unelastischen Formanderungen x} abhangig ist. Man pflegt in der Rheologie auch dies auszu- 
schlieBen und hat dann f(x%) = 0 zu setzen. Wir machen im folgenden von beiden Vereinfachungen 
Gebrauch, indem wir U = U(xf, 0) und F = F(xj, 0) annehmen. 

Um nun zunachst das elastische Verhalten des Kérpers zu diskutieren, kann man sich (wie dies 
in der Elastizitatstheorie iiblich ist) unter Berufung auf die meist geringe Warmeentwicklung und 
den Warmeaustausch mit der Umgebung auf isotherme Zustandsanderungen beschranken. Die 
Temperatur fallt dann als Zustandsvariable weg, und die freie Energie, die jetzt nur noch von den 
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x¢ abhangt, geht in die gewohnliche potentielle Energie der inneren Krafte, d. h. in die Form- 
ainderungsenergie (xf) tiber, aus der sich die Spannungen gemah (4.8) mit (3.8) berechnen. Setzt 
man die Zustandsanderungen umgekehrt als adiabatisch voraus, so folgt mit dx; = 0 und dQ = 0 
aus (4.1) und (4.2) dU = X;,, dxj. Die innere Energie ist also auch hier von der Temperatur unab- 
hangig, und die Spannungen sind, wenn man U = @P(x~) setzt, auch in diesem Fall durch (3.8) ge- 
geben. 
¢ Es stellt sich also heraus, da® die Beziehungen (3.8) das allgemeinste rheologische Gesetz des 
isotherm- bzw. adiabatisch-elastischen Kérpers reprasentieren. Man kann es dadurch linearisieren, 
daB man die potentielle Energie als quadratische Form (3.9) ansetzt, und kommt damit auf den 
Hooke-Kérper mit seinen Symmetriebedingungen (3.10) zuriick, die ibrigens als Umkehrungen des 
Maxwellschen Reziprozitatssatzes gedeutet werden kénnen. 

Treten auch unelastische Deformationen auf, dann gilt in (4.4) das Ungleichheitszeichen. Man 
hat also nach (4.3) bis (4.5) 


dQ = bdS — X,, dxi < OdS. (4.9) 
Fiir die sogenannte Dissipationsarbeit dA’ = X, dx; ergibt sich hieraus die Ungleichung 
dAt = X,, dxi > 0 (4.10) 
und fiir die Entropiezunahme die Beziehung 
dS == (dQ + dA’). (4.11) 


Aus (4.10) und (4.11) folgt sodann wie erwartet, daB fiir dQ = 0 jede unelastische Formanderung die 
Entropie vergréBert. 
Man kann im AnschluB an die Carnot-Clausiussche Relation 


dS = dS* + dst = . + dsi (4.12) 
die GréBe dS* = dQ/0 als Entropiezufuhr von auBen und den mit (4.11) umgeformten zweiten Term 
ee el 
dS' = ia: — a X;, dxy = 0 (4.13) 
als Entropieproduktion im Inneren des Kérpers auffassen. Mit der Dissipationsleistung 
Li xX, 0 (4.14) 


ergibt sich dann aus (4.13) die spezifische (d.h. auf die Zeiteinheit bezogene) Entropiepro- 
duktion zu 
gil _ly ew A.15 
See pap Ore (4.15) 
Wahrend diese thermodynamischen Uberlegungen, die ausschlieBlich auf den beiden Haupt- 
satzen beruhen, das thermoelastische Verhalten unserer Kérper vollstandig beschreiben und unter 
den erwahnten Vereinfachungen auf den Hooke-Kérper zuriickfiihren, sagen sie tiber das unelastische 
Verhalten und damit insbesondere fiir unsere iibrigen rheologischen Kérper noch recht wenig aus. 
Man kann indessen fiir diese Kérper von (4.15) aus weiter schlieBen, und zwar zunachst fiir den Fall, 
daB die Spannungen X; und die unelastischen Formanderungsgeschwindigkeiten «* durch lineare 
Bezichungen miteinander verkniipft sind, mit Hilfe der Onsagerschen Theorie. 


5. Die Onsagerschen Relationen. Die Onsagersche Theorie! fiihrt wesentlich iiber die beiden 
Hauptsatze der Thermodynamik hinaus. Sie bezieht sich auf beliebige irreversible Prozesse, die 
sich durch ein lineares Gesetz beschreiben lassen. Obschon sie sich nur durch statistische Uber- 
legungen begriinden [bzw. auf eine einfache und plausible Annahme zuriickfithren laBt]2, kann man 
sie rein phanomenologisch folgendermafen formulieren: 

An einem irreversiblen Vorgang seien neben der Temperatur 9 die unabhangigen Zustands- 
variablen x, beteiligt. Ist die spezifische Entropieproduktion durch 


ari 1 . ‘ 
S = 5 Xn Xk» (5.1) 


die Dissipationsleistung mithin durch 


Li = X,, X, (5.2) 


1 5. FuBnote 3, S. 58. 
2s, FuBnote 9, S. 62. 
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gegeben, so kann man die X; (oder auch die X;/9) als die zu den x, konjugierten Zustandsvariablen 
bezeichnen. Wenn nun der Vorgang durch einen Satz linearer Beziehungen 


X; = ay, %, (5.3) 


zwischen den konjugierten Variablen beschrieben werden kann, dann miissen die hierin auftreten- 
den Konstanten den Symmetriebedingungen 


% a, = Gi, (5.4) 
geniigen. 

Die linearen Beziehungen (5.3) werden in diesem Zusammenhang als das phanomenologische 
Gesetz des irreversiblen Vorgangs, die Symmetriebedingungen (5.4) auch als Onsagersche Rezi- 
prozitatsrelationen bezeichnet. 

Ks ist klar, daB die Onsagersche Theorie gerade fiir rheologische Vorgange, wie wir sie hier im 
Auge haben, eine wichtige und auf anderem Wege nicht zwingend zu begriindende Aussage enthalt. 
Die x, sind namlich in unserem Fall die unelastischen Formanderungsgeschwindigkeiten, in denen 
nun nach Erledigung des elastischen Verhaltens der obere Zeiger u fortan weggelassen werden kann, 
und wenn die X; mit den Spannungen identifiziert werden, sind (5.1) und (5.2) mit (4.15) bzw. (4.14) 
gleichbedeutend. Nach der Onsagerschen Theorie gelten nun, wenn das rheologische Gesetz die 
lineare Gestalt (5.3) hat, die Symmetriebedingungen (5.4). 

Damit ist jetzt auf thermodynamischem Wege die zum rheologischen Gesetz (5.3) des Newton- 
Kérpers gehérende Symmetrie 


Ae = ix (5.5). 


physikalisch begriindet und somit auch die viskoelastische Analogie fiir lineare rheologische 
Gesetze, wie sie der Newton- und der Maxwell-Kérper aufweisen. Auch die in Bemerkung 4 
von Ziff. 2, diskutierte Modellvorstellung fiir den (linearen) Newton-Korper, die sich eines Systems 
viskoser Einheiten gemaB Abb. 2 bedient, findet damit nachtraglich ihre Rechtfertigung!. 


Die mit (5.4) ausgedriickte Symmetrie iibertragt sich auch auf die Umkehrung 
i = bX, mit by = (5.6) 
des phanomenologischen Gesetzes (5.3), die aus physikalischen Griinden als existent vorausgesetzt 
werden darf. Mit Riicksicht auf die Symmetrie der Matrix (a;,) ist nimlich A,; = A,, und daher 
auch 


by, = by; (5.7) 


Beziehungen von der Art der Onsagerschen Relationen (5.4) bzw. (5.7) kommen in der Mechanik 
auf den verschiedensten Gebieten vor und sind dort stets der Ausdruck fiir das Vorhandensein eines 
Potentials. Man denke zum Beispiel an die Maxwellschen Reziprozititsrelationen oder an die dazu 
dualen Symmetriebedingungen (3.10) des Hooke-Kérpers. Es liegt daher nahe, auch hinter den 
Onsagerschen Beziehungen ein Potential und damit vielleicht einen allgemeineren physikalischen 
Inhalt zu suchen. 

In der Tat kann man die Beziehungen (5.4) mit dem linearen Gesetz (5.3) in der Form 


OX% __ OXi (5.8) 


ax; Oxy 


schreiben, und hieraus folgt, worauf schon L. Onsager hingewiesen hat, daB sich das phinomeno- 


logische Gesetz (5.3) mit Hilfe des Potentials 


P(X;,) = B Ay, X; %, (5.9) 
in der Gestalt 
eee (5.10) 
i Ox. 


darstellen laBt. Fiir unser rheologisches Problem bedeutet dies, dab der Spannungsvektor X; im 
Geschwindigkeitsraum x, wirbelfrei und damit der Gradient einer Funktion P(x,) ist, die im Folgen- 
den kurz als Geschwindigkeitspotential bezeichnet werden soll. 


1 Wie der Verf. wahrend der Drucklegung dieser Arbeit erfahren hat, ist diese Feststellung oe bite 
Zeit schon von M.A. Biot gemacht und in verschiedenen Richtungen ausgewertet worden. ee % we ; 
J. Appl. Phys. 25 (1954), S$. 1385, Phys. Rev. 97 (1955), S.1463 sowie J. Appl. Phys. 27 ( ), S. 240. 
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Ebenso gilt fiir die Umkehrung (5.6) des phanomenologischen Gesetzes nach (5.7) 


ew 
so daf (5.6) mit dem Potential 
V(X) => bys XX, (5.12) 
auf die Gestalt . 
i) 
Ky, =a af 
EBS, (5-43) 


gebracht werden kann. Dabei ist der Zusammenhang zwischen den Funktionen @(x,) und ¥/(X;) 
nach (5.12), (5.3), (5.6), (5.4) und (5.9) durch 


1 ine 1 Aix ee 1] felts 1 a. ; 
P(X) = > bei IIMB eal mamar hn Cel 8 tad at id = Oi CFM kod faa ey ear ist a P(X;,) (5.14) 


Jj 
gegeben, wenn x, und X; durch das phanomenologische Gesetz verkniipfte Werte sind. In unserem 
Fall heiBt dies, da® der Vektor der Formanderungsgeschwindigkeit x, im Spannungsraum X, wirbel- 
frei und damit der Gradient einer Funktion Y/(X;) ist, die als Spannungspotential bezeichnet 
werden soll. 

Das phanomenologische Gesetz (5.3) vermittelt eine umkehrbar eindeutige Abbildung des 
Raums %, (im rheologischen Fall also des Geschwindigkeitsraums) auf den Raum X; (den Span- 
nungsraum). Die Umkehrung dieser Abbildung wird durch (5.6) gegeben. Nach (5.14) folgt aus 
@® = konst. auch Y = konst. und umgekehrt; die Potentialflachen im einen Raum sind demnach 
die Bilder derjenigen im anderen Raum. 

Beilaufig schlie8t man sowohl aus den Symmetriebedingungen (5.4) wie aus der Darstellung 
(5.9), (5.10) darauf, daB sich das phanomenologische Gesetz und damit auch seine Umkehrung auf 
Hauptachsen transformieren lassen. 

SchlieBlich muB die Dissipationsleistung (5.2), die gemaB (5.3) bzw. (5.6) in jeder der beiden 


Formen 


it@ — Gis, Lp oa = 2} DP, ike — bp X;, X; = 2 WV (5.15) 


geschrieben werden kann, mit der spezifischen Entropieproduktion nach (4.14) oder (4.15), d. h. 
auf Grund des zweiten Hauptsatzes der Thermodynamik, fiir beliebige nichtverschwindende 
Vektoren %, bzw. X; positiv sein. Die Onsagerschen Relationen kénnen also durch die weiteren 
Bedingungen erganzt werden, da die Matrizen 
(a;,) und (b,;) positiv definit sind. Hieraus folgt 
in Verbindung mit (5.15), daB die Potentialflachen 
@® = konst. baw. Y = konst. Ellipsoide sind. 
Abb.6 zeigt die Verhaltnisse bei zwei Dimen- 
x, Sionen im Hauptachsensystem und liefert bereits 
einen Hinweis auf den Zusammenhang mit den 
Fragen, die in Ziff.3 in bezug auf die Misessche 
Theorie des plastischen Potentials aufgeworfen 
worden sind: die Flachen Y = konst. sind konvex, 
und die ihren Punkten zugehérigen Vektoren x, 


der Formanderungsgeschwindigkeit sind normal 
Abb. 6. Durch ein lineares rheologisches Gesetz vermittelte Abbildung, zu ihnen 


‘Y=const 


Auf Grund der nunmehr fiir den linearen Fall gesicherten viskoelastischen Analogie lassen sich 
natiirlich diese Uberlegungen analog auch fiir einen linearen Hooke-Korper anstellen. So verkér- 
pert (5.10), wenn man x, durch x, ersetzt und O(x,) baw. Y(X;) als Formanderungsenergie deutet, 
das Prinzip der virtuellen Leistungen und (5.13) das Prinzip von Castigliano. 


6. Eine Hypothese. Die Onsagerschen Relationen (5.4) haben nur fiir lineare phanomenologische 
Gesetze einen Sinn. Dagegen ist die Onsagersche Theorie in der durch die Beziehungen (5.8) bis 
(5.10) gegebenen Gestalt erweiterungsfahig. Sie fordert eine Verallgemeinerung auf nichtlineare 
Vorgange geradezu heraus, da die Beziehungen (5.10) auch fir Potentiale @(x,) sinnvoll bleiben, 
die nicht quadratische Formen in den %, sind. 

Nimmt man versuchsweise an, daf bei irreversiblen Vorgangen mit beliebigem phanomenologi- 
schem Gesetz 


X; = X,(%) (6.1) 
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die Beziehungen (5.8) giiltig bleiben, so bedeutet dies den Ersatz der Onsagerschen (5.4) durch die 

allgemeineren Relationen ; 
aX, aX; 
ae, Ok, 


(6.2) 


Diese Beziehungen, welche die Wirbelfreiheit des Feldes X,(%,) postulieren, sind mit der Aussage 
gleichwertig, daB im x,-Raum ein Potential 


D(%,) = | X; dx; (6.3) 
existiere, mit dem sich das phanomenologische Gesetz (6.1) in der Form 
_ a@ 
i= on (6.4) 


anschreiben 1a8t. Dabei kann das Potential (6.3) durch Integration langs einer beliebigen Kurve 
im x,-Raum gebildet werden, die einen festen Ausgangspunkt (beispielsweise den Ursprung) mit 
dem gewahlten Aufpunkt verbindet. 

Diese Verallgemeinerung, die im Falle eines linearen phanomenologischen Gesetzes offensicht- 
lich auf die Onsagersche Theorie zuriickfiihrt, ist natiirlich vorerst eine bloBe Hypothese. Wahrend 
namlich in der Mechanik die Wirbelfreiheit eines Kraftfeldes in allen Fallen, in denen die Erhaltung 
der mechanischen Energie feststeht, aus dem Energieprinzip folgt, fehlt in der phanomenologischen 
Thermodynamik ein Prinzip, aus dem sich die mit (6.2) postulierte Wirbelfreiheit des Spannungs- 
feldes ableiten lieBe. Vielleicht laBt sich aber die vorgeschlagene Verallgemeinerung ahnlich wie die 
Onsagersche Theorie auf statistischem Wege mit dem Prinzip der mikroskopischen Reversibilitat 
begriinden. : 


Hier seien nur die folgenden Anhaltspunkte fiir die Richtigkeit der Hypothese notiert: 


Erstens bildet sie, wie schon bemerkt wurde, eine naheliegende Verallgemeinerung der Onsager- 
schen Theorie. 

Zweitens reduziert sie sich in der Rheologie auf die ebenfalls naheliegende Forderung, daB sich 
die in Bemerkung 4. von Ziff. 2, diskutierte Modellvorstellung fiir den linearen Newton-Kérper, die 
sich eines Systems viskoser Einheiten gemaf3 Abb. 2 bedient und durch die Onsagersche Theorie 
nunmehr gerechtfertigt ist, auch im Falle des nichtlinearen Newton-Kérpers verwenden lasse. Beim 
nichtlinearen Hooke-Kérper fiihrt namlich die Vorstellung eines Modells, das sich aus nichtlinearen 
elastischen Einheiten gemaB8 Abb. 1 [mit X = X(x) an Stelle von (2.2)] zusammensetzt, auf das 
rheologische Gesetz , 


oD 
1G ie (6.5) 
mit der Formanderungsenergie 
D(x,) = [ X; dx; , (6.6) 


und wenn man annimmt, daB sich in analoger Weise das Modell des nichtlinearen Newton-Korpers 
aus nichtlinearen viskosen Elementen gemaB Abb. 2 [mit X = X(x) an Stelle von (2.3)| aufbauen 
lasse, dann kommt man von (6.5), (6.6) aus direkt auf das durch (6.3) erginzte rheologische Gesetz 
(6.4). 

Drittens ist innerhalb der Rheologie nach (6.5) und (6.6) unsere Hypothese, d. h. die Annahme, 
daB der nichtlineare Newton-Kérper den Beziehungen (6.3) und (6.4) gehorche, mit dem Postulat 
der unbeschrankten Giiltigkeit der viskoelastischen Analogie gleichbedeutend. Fiir dieses Postulat 
spricht aber einerseits der Erfolg, mit dem N. J. Hoff1 die Analogie fiir immerhin recht allgemeine 
nichtlineare isotrope Kérper etabliert hat, andererseits aber auch die in Ziff. 9 zu besprechende 
Méglichkeit, die Misessche Theorie des plastischen Potentials darauf zu griinden. 


7. Anwendung auf den nichtlinearen Newton- und Maxwell-Kérper. In Ziff. 3 wurde fiir einen 
verallgemeinerten Newton-Kérper, dessen rheologisches Gesetz im eindimensionalen Fall durch 
(3.6) gegeben ist, im zweidimensionalen Fall der Ansatz (3.7) gemacht. Auf Grund unserer Hypo- 
these bzw. der Forderung, daB der betrachtete Kérper der viskoelastischen Analogie geniige, kann 
man aber vom Geschwindigkeitspotential 


D = a, 4 + ay 43 ey 4 Gg XP XZ + ay XX} + as x (7.1) 


15, FuBnote 4, S. 58. 
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ausgehen und erhalt dann statt (3.7) nach (6.4) 
X, = 4a, 3 + 3 a, #2 %, + 2 ag X43 + a, X83, 
X, = a, %? + 2a, 42%, + 3 a, X, 45 hss 
mithin ein rheologisches Gesetz, das im Vergleich zu (3.7) statt acht nur noch fiinf Konstanten 


enthalt. 
Man kann diesen SchluB natiirlich verallgemeinern. Bei einem Kérper mit n Dimensionen, 
dessen rheologisches Gesetz im eindimensionalen Fall die Gestalt 


X=ca? (p = 4, 6,...) (7.3) 


(7.2) 


hat, wird man das Geschwindigkeitspotential mit Hilfe der nichtnegativen Zahlen a,, X, .-- 5 
in der Form 


DE Ny canteen pats iat ne (xo, +o, +... +a, =p) (7.4) 


O12 Kgreery Ky 


ansetzen und das allgemeine rheologische Gesetz hieraus mit (6.4) ableiten. 
Damit ist auch das Vorgehen beim nichtlinearen Maxwell-Kérper gegeben, bei dem selbstver- 
standlich die elastischen und die unelastischen Formanderungen getrennt zu behandeln sind. 


8. Komplementire Potentiale. Mit Hilfe der Legendreschen Transformation! lassen sich die im 
zweiten Teil von Ziff. 5 gewonnenen Ergebnisse fiir beliebige Potentiale und damit insbesondere fiir 
beliebige nichtlineare Newton-Kérper verallgemeinern, fiir welche eine viskoelastische Analogie 
besteht. 


Nach (6.4) 1aBt sich das rheologische Gesetz eines solchen Kérpers in der Form 


oP 
x= ¥, (8.1) 
aus dem Geschwindigkeitspotential @(x,) ableiten. Dabei soll zunachst angenommen werden, dai 
das phanomenologische Gesetz (8.1) eindeutig sei, und zudem kann die Funktion @(*,) so normiert 
werden, daB sie im Ursprung o des %,-Raums verschwindet. Die Beziehungen (8.1) vermitteln 
dann eine eindeutige Abbildung eines Bereiches g des x,-Raums [des Definitionsbereiches von 
@(x,)| auf einen Bereich G im X,;-Raum, wie sie in Abb. 7 fiir zwei Dimensionen angedeutet ist. 
Die Umkehrung des phanomenolo- 
gischen Gesetzes 
i, = &(X,) (8.2) 
soll ebenfalls zunachst als eindeutig 
vorausgesetzt werden (so daB die Funk- 
tionaldeterminante |0X;/0x,| im ganzen 
Bereich g von Null verschieden ist). Sie 
vermittelt ihrerseits eine eindeutige Ab- 
bildung von G auf g. Daraus, daB® nach 
Ziff. 4 irreversible Zustandsanderungen 
ohne Dissipationsleistung ausgeschlossen 
sind, folgt, daB der Ursprung O im X;- 
Raum der Bildpunkt von o sein mu, so 
Abb. 7. Durch ein allgemeineres rheologisches Gesetz vermittelte Abbildung. daB insbesondere g den Ursprung o und 
G den Ursprung O enthalt. 
Ist c eine beliebige Kurve in g, die 0 mit einem Punkt p mit den Koordinaten x, verbindet, und 
ist ferner der Punkt P mit den Koordinaten X,, der Bildpunkt von p in G, dann verbindet das Bild C 
von ¢c in G die Punkte O und P, und es gilt 


[ X, dk, + [ dX; = f d(X;%) =X, x,. 
¢ G 


X3 


(8.3) 
Unter Beniitzung von (8.1) erhalt man hieraus 


{OR aps (8.4) 


1 R. Courant und D. Hilbert, Methoden der Mathematischen Physik, Bd. II, S. 26 ff, Berlin 1937; s. auch 
FuBnote 1, S. 69. é 
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Denkt man sich rechterhand x, vermittelst (8.2) durch X, ausgedriickt, dann sieht man, daB das 
Integral linkerhand eine Funktion 


[ %, dX, = ¥(X,) 
Cc 


(8.5) 
der Koordinaten X, von P ist, und daB 
: ow 
= (8.6) 


k 
gilt. Somit laBt sich die Umkehrung (8.2) des rheologischen Gesetzes von einem Spannungspotential 
P(X;) ableiten, das nach (8.5) und (8.4) bereits so normiert ist, daB es in O verschwindet. 

Man kann ¥(X;) das zu @(x,) komplementare Potential nennen und (8.1), (8.6) in die Fest- 
stellung zusammenfassen, daB® sich auch im nichtlinearen Fall im Geschwindigkeitsraum der Span- 
nungsvektor X; als Gradient des Geschwindigkeitspotentials und im Spannungsraum dex Vektor %, 
der Formanderungsgeschwindigkeit als Gradient des komplementaren Spannungspotentials dar- 
stellen laBt. 

Ersetzt man unter Berufung auf die viskoelastische Analogie «, durch x,, dann stellt. D(x,) 
die Formanderungsenergie und (8.1) das Prinzip der virtuellen Leistungen, Y/(X;) die sogenannte 
komplementire Energie und (8.6) das von F’. Engesser! fiir nichtlineare clastische Kérper ver- 
allgemeinerte Prinzip von Castigliano dar. 

Die Dissipationsleistung L? = X;,, x,,, die sich sowohl im x,- wie im X;-Raum als Skalarprodukt 
aus dem Fahrstrahl und dem Gradienten des Potentials auffassen ]aBt, ist nach (4.14) positiv. 
Hieraus folgt, da die Potentialflachen in beiden Raumen beziiglich der Urspriinge konkav sind, 
und da die Potentiale auf jedem Radius nach auBen monoton anwachsen. Die Potentiale D(x,) 
bzw. ¥/(X,) sind demnach eindeutige Funktionen ihrer Argumente (und zwar offensichtlich auch 
dann, wenn das rheologische Gesetz (8.1) oder/und seine Umkehrung (8.6) nicht eindeutig sind). 

Die Potentialflachen ® = konst. und Y = konst. sind aber im nichtlinearen Fall keine Ellip- 
soide mehr und bilden sich im allgemeinen auch nicht mehr aufeinander ab. Fiir die Transformation 


zwischen @(x,) und Y/(X;) gilt vielmehr nach (8.4) und (8.5) die Beziehung 
D(x,) + Y(Xy) = Xp Xe» (8.7) 


in der die rechte Seite wiederum die Dissipationsleistung L‘ darstellt. Aus ® = konst. folgt nur 
dann Y = konst. und umgekehrt, wenn auf der betrachteten Flache auch L* einen konstanten 
Wert besitzt, wie dies nach (5.15) beispielsweise im linearen Fall zutrifft. 

Von besonderem Interesse ist der Fall, daB alle O-Flachen im x,-Raum ahnlich sind und beziig- 
lich des Ursprungs o ahnlich liegen (wie dies z. B. wiederum beim linearen Kérper zutrifft). Fihrt 
man hier von verschiedenen Punkten x, einer beliebigen ®-Flache aus radiale und dem Abstand 
von 0 proportionale Schritte dx, aus, so das mit festem de 

dt == &, de (8.8) 
gilt, dann ist der zugehorige Zuwachs von @ nach (8.1) und (8.8) durch 
db = 2 di, = X, de (8.9) 
k 


gegeben, und da die neuen Punkte (1 + de) %, voraussetzungsgemaB wieder auf einer D-Flache 
liegen, muB d@ fiir alle Ausgangspunkte %, gleich, mithin X;, x, eine Funktion 


von @ allein sein. Mit (8.10) folgt aber nach (8.7) aus ® = konst. auch YW — konst. und umgekehrt; 
in diesem Sonderfall bilden sich daher die ®- und die ¥/-Flachen aufeinander ab. 
Mit (8.10) geht (8.9) in 
d® = f(®) de (8.11) 
iiber. Die Bildpunkte X,, der x, im X;-Raum, die urspriinglich alle auf der gleichen Y-Flache liegen, 
verschieben sich mit (8.8) nach (8.1) und (8.11) um 
: pS EE 7 ees — ¢(G) oo de = X,f'(P) de, (8.12) 
CMe ae Gq OP 5, J) 0 iC ea he) 


1 F, Engesser, Z. Arch. Ing. Ver. Hannover 35 (1889) S. 733. 
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wobei f'(®) die Ableitung von f(®) nach @ bedeutet. Diese Verschiebungen sind aber radial und 
den Abstanden von O proportional, so da auch alle Y-Flachen im X,;-Raum ahnlich sind und 
beziiglich des Ursprungs O ahnlich legen. 

Mit Riicksicht darauf, daB® das Potential @(%,) auf jedem Radius im *x,-Raum nach auBen 
monoton zunimmt, schneidet ein gegebener Strahl von o aus eine gegebene ®-Flache héchstens in 
einem Punkt p. Es gibt also auf der zugehérigen Y’-Flache héchstens einen Punkt P mit gegebener 
Normalenrichtung. Hieraus folgt aber, daB die Y-Flachen konvex sind, und ganz analog laBt sich 
die Konvexitat der ®-Flachen nachweisen. 


9. Anwendung auf den Pranditl-Reuss-Kérper!. Die FlieBflache (3.11) eines Prandtl-Reuss- 
Korpers, die wir zunichst als regular voraussetzen, ist aus physikalischen Griinden geschlossen und 
enthalt den Ursprung O des Spannungsraums. Fiir Spannungsvektoren X;, deren Endpunkte 
innerhalb der FlieBflache H liegen, treten keine plastischen Formanderungen auf; die durch das 
rheologische Gesetz vermittelte Abbildung ist also derart, da dem Ursprung o im *,-Raum das 
ganze Innere von H im X;-Raum entspricht. Da ferner zu jedem plastischen Flu8 ein Spannungs- 
vektor gehért, dessen Endpunkt auf H liegt, wird der ganze x,-Raum mit Ausnahme des Ursprungs 0 
im X,-Raum auf die FlieBflache H abgebildet. An die Stelle des Bereiches g von Ziff. 8 tritt hier also 
der ganze «,-Raum; der Bildbereich G umfaBt die FlieBflache H und ihr Inneres, und weder das 
rheologische Gesetz noch seine Umkehrung ist eindeutig. 

Es liegt indessen schon aus physikalischen Griinden nahe, den Prandtl-Reuss-Kérper als Grenz- 
fall einer Folge von nichtlinearen Newton-Kérpern K” aufzufassen, deren Viskositat innerhalb H 
mit zunehmendem n iiber alle Grenzen wachst, mit der Annaherung an die Flache H aber auf Null 
abfallt. Setzt man dabei noch voraus, daB die Kérper K” der viskoelastischen Analogie geniigen, 
dann kann man nach Ziff. 8 jeden derselben im X,-Raum innerhalb H durch eine feste Anzah] m 
aequidistanter Y-Flachen dargestellt denken. Diese Flachen schmiegen sich mit zunehmendem n 
von innen her der Flache H immer mehr an, und zwar aus Regularitatsgriinden derart, daB fur 
genugend groBes n die Punkte P; (j =1, 2,...,m) bzw. P, in denen ein beliebiger Radius die 
Flachen Y = konst. sowie H schneidet, beliebig nahe beisammen liegen, und da sich ferner die in 
diesen Schnittpunkten auf den zugehérigen Flachen errichteten Normalen in ihren Richtungen 
beliebig wenig unterscheiden. 

Ks ist dann klar, daB mit wachsendem n auf jedem Radius die Richtung des Gradienten 0Y/dX; 
gegen die Richtung von 0F'/0X; im Punkte P strebt, und daB sein Betrag bei der Annaherung an 
P alle Werte in einem Intervall annimmt, dessen Grenzen gegen 0 und oo streben. Im Grenzfall 


gilt daher nach (8.6) 
(9.1) 


mit A < 0, und damit ist die Misessche Grundbeziehung (3.12) begriindet. 

Setzt man zusatzlich voraus, daB sich der geschilderte Grenziibergang durch Kérper K” reali- 
sieren lassen miisse, deren ®- und Y-Flachen ahnlich sind und liegen, dann folgt in Verbindung mit 
Ziff. 8, daB sie gegen eine konvexe Flache konvergieren. Somit ist auch die Konvexitat der FlieB- 
flachen zwar nicht einwandfrei begriindet, aber doch immerhin plausibel gemacht. 

Da man allfallige irregulare FlieBflachen durch Abrunden von Ecken und Kanten stets regulari- 
sieren und den damit verbundenen Eingriff physikalisch belanglos halten kann, lassen sich die vor- 
stehenden Schliisse auch auf solche Falle erweitern. Auch der Ubertragung auf plastische Kérper 
mit Spannungsverfestigung steht nichts im Wege. 


(Eingegangen am 22. Mai 1956) 
Anschrift des Verfassers: Prof. Dr. Hans Ziegler, Riischlikon (Ziirich), Weiherweg 6. 


* Eine auf einer anderen Annahme berahende Begriindung der Misesschen Theorie ist 1950 von 
D.C. Drucker gegeben worden. Vgl. Quart. Appl. Math. 8 (1950) S. 411. 
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Ermiidungserscheinungen sind dadurch gekennzeichnet, daB viele Werkstoffe bei lange andauernder, insbesondere periodisch 
mechanischer Beanspruchung in ihrem Widerstand dagegen nachlassen. Zahlreiche Beobachtungen liegen u. a. vor im Ma- 
schinenbau, im Verkehrswesen (Kraftfahrzeug, Eisenbahn, Schiffabrt, Flugzeug) und im Bauwesen (Briickenbau). 


Das vorliegende Werk enthalt eine Sammlung von Vortragen dieses Gebietes, die auf einer Tagung gehalten wurden, welche 


das Kénigliche Institut fiir Technologie in Stockholm vom 25. bis 27. Mai 1955 auf Veranlassung der Internationalen Union 
fiir Theoretische und Angewandte Mechanik veranstaltete. Diese 35 Vortrige international anerkannter Wissenschaftler 
reprasentieren den gegenwartigen Wissensstand auf dem Gebiete der: Ermiidungsfestigkeit. Der Inhalt der Darlegungen 
konzentriert sich im wesentlichen auf folgende Themen: Statistische Theorie der Ermitdungserscheinungen, ihre Schadens- 
hiufigkeit, die Grundlagen des Mechanismus, metallkundliche Betrachtungen, kinetische Uberlegungen, die Abhangigkeit 
von héheren Temperaturen und das gleichzeitige Auftreten mehrfacher Beanspruchung einzelner Teile von Maschinen und 
Bauwerken, . ; \ 

Die vorliegende erste abgeschlossene Darstellung dieses jungen Wissenszweiges laBt ein weitgehendes Interesse auf den ver- 
schiedensten Gebieten der Technik (Maschinenbau, Bauwesen, Materialkunde und -priifung) erwarten. 


Contents Sommaire Inhalt 


¥, BASTENAIRE: Etude critique de la notion de dommage appliquée 4 une classe étendue d’essais de fatigue. ; 

F. BASTENAIRE, R. CAZAUD, M. WEISZ: Essais de fatigue statistiques suivant la méthode de charge progressive, 

A. FERRO, U. ROSSETTI: Contribution a l’étude de la fatigue des matériaux avec essais A charge progressive. 

W. N. FINDLEY: Theories relating to fatigue of materials under combinations of stress. : 

A. FRANSSON: Effect of simultaneous cyclic variation of stress and temperature on a high temperature material. 

A.M. FREUDENTHAL: Physical and statistical aspects of cumulative damage. 

E. GASSNER: Ermiidungsfestigkeit bei statistisch verdnderlichen Spannungsamplituden. 

F. GATTO: New statistical methods applied to the analysis of fatique data. 

M. HEMPEL: Uber Verformungserscheinungen in Stahlen bei der Wechselbeanspruchung, 

R. B. HEYWOOD: The effect of high loads on fatigue. 

R. JACQUESSON: Modifications de texture cristalline produites par des efforts. alternés. 

A. JOHANSSON: Fatigue of steels constant strain amplitude and elevated’ temperature. 

A. KAMMERER: Une définition théorique de la limite de fatigue. : 

P. KUHN: Effect of geometric size on notch fatigue. SP ase : 

P. LAURENT: Les travaux récents de l'Institut de Recherches Métallurgiques de Sarrebruck dans le domaine de la fatigue. 

O. LISSNER: Einige Versuche iiber die Vorgiinge in der Oberflichenschicht yon Ermiidungsproben. 

L. LOCATI: Essais de fatigue par flexion avec fréquence’ superposées. 

L. MARTINAGLIA: Ziele der Ermiidungsforschung in der Schweiz, gezeigt am Beispiel von Dauerversuchen an Schrauben- 
verbindungen. 

F, A. McCLINTOCK: Variability in fatigue testing: sources and effect on notch sensitivity. 

I. A. ODING: Uber den Mechanismus der Zerstérung bei der zyklischen Belastung von Metallen. 

R. E. PETERSON: Torsion and tension relations for slip and fatigue. 

A. I. PETRUSSEWITSCH: Einige Besonderheiten der Kontaktermiidung. \ 

C, E, PHILLIPS: Some observations on the propagation of fatigue cracks. : 

F. J. PLANTEMA: Some. investigations on cumulative damage, 

J.T. RANSOM: A guide to statistical methods for use in fatigue testing. 

J. oe eee Die durch grofe Zug- und Druckermiidungsbelastungen hervorgebrachte mechanische Hysterese in 

tahlen, 

C. SCHAUB, W. LIEDTKE: Der Mechanismus des Dauerbruchs metallischer Werkstoffe. ‘ 

F.R. SHANLEY: A proposed mechanism of fatigue failure. y 

D. G. SOPWITH: Recent researches on fatigue at the Mechanical Engineering Research Laboratory, East Kilbride. 

R. J. TAYLOR: Experimental design and methods of analysis used in studying effects of metallurgical variation on fatigue. 

G. V. UZHIK: A contribution go the theory of the fatigue of metals. ‘ ; 

W. WEIBULL: Basic aspects of fatigue. 

F. WEVER: Ausscheidungsvorgiinge in Staihlen bei ruhender und wechselnder Beanspruchung. 

P. E. WIENE: Bending fatigue of large welded test pieces with regard to velocity of crack propagation. 

E, W. C. WILKINS: Cumulative damage in fatigue. : 


‘Subject index. Table des matiéres. Sachverzeichnis. 
SUEIEREnRIGRET URERARERRTET Eur mnencammencerepee es eae ease Ae NOS UE Sn or Mee MONEY Saat ANNL Se LS SU TEI ME I) tJ ey 
SPRINGER - VERLAG. BERLIN / GOTTINGEN / HEIDELBERG 


Springer-Verlag + Berlin / Géttingen / Heidelberg. Printed in Germany. V/12/6 0,8 K.B./Z. 040, 


ee ee 


ead 


